SYNTHETIC-APERTURE RADAR IMAGING OF THE OCEAN SURFACE: THEORETICAL CONSIDERATIONS, AND EXPERIMENTS WITH SIMULATED AND ACTUAL SAR IMAGERY By PARIS W. VACHON B.A.Sc, The University of British Columbia, 1983 A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY in THE FACULTY OF GRADUATE STUDIES (Department of Oceanography) We accept this thesis as conforming to the required standard THE UNIVERSITY OF BRITISH COLUMBIA August 1987 ©Paris W. Vachon, 1987 In presenting degree this at the freely available copying of department publication of in partial fulfilment University of British Columbia, for reference and study. this or thesis thesis by this for scholarly his thesis or her the I agree requirements for may representatives. It be is Department The University of British Columbia 1956 Main Mall Vancouver, Canada V6T 1Y3 advanced that the Library shall make it I further agree that permission purposes an granted by the understood for extensive head that for financial gain shall not be allowed without permission. DE-6(3/81) of of my copying or my written Abstract Three key areas of controversy in synthetic-aperture radar (SAR) imaging of ocean surface waves are considered: first, the nature of Bragg scattering; second, the role, magnitude, and calculation of the scene coherence time; and third, the relevant ocean wave velocities for coherent Doppler modulations. This work begins with a re-derivation and extension of existing SAR imaging theory for point and diffuse targets. Generic, relatively simple, closed-form expressions for the impulse response, the resolution, and the image bandwidth summarize this unified treatment. Theoretical differences between the imagery of point and diffuse targets are pointed out. Based upon these fundamental differences, a statistical testing procedure is formulated to address the question of scene target density. Background ocean surface wave theory is outlined in preparation for discussions of SAR ocean imaging. Of central importance is the role of the phase velocity, which is the speed of translation of the mean pattern of reflectivity, and the orbital motion, which leads to coherent (phase) modulation, and hence to velocity bunching, acceleration defocus, and target decorrelation. Based upon this theoretical background, one- and two-dimensional simulation models are developed. The one-dimensional simulation addresses the effects of various parameters upon the mean image contrast in a velocity bunching model and guides the development of the two-dimensional simulation. The two-dimensional simulation is unique because each target which constitutes the scene is explicitly considered. This leads to a degree of control andflexibilitywhich is not available from actual SAR imagery. Qualitative and quantitative comparisons are drawn between the simulated and actual SAR imagery to address the key areas of controversy. The assertion that Bragg scattering is a coherent process is defended, despite inability to conclusively verify this ii using SEASAT data. Comparisons between simulation and C-SAR imagery of waves propagating into ice verify the roles of the scene coherence time and the wave phase velocity. iii Contents Abstract ii Contents iv L i s t of Tables vii L i s t of Figures viii L i s t of S y m b o l s ix Acknowledgements xii 1 Introduction 1 2 S A R P o i n t Target Imaging 4 2.1 2.2 2.3 2.4 3 5 28 Second Order Image Statistics 3.1.1 Static Diffuse Target 3.1.2 Dynamic Diffuse Target Point vs. Diffuse Scattering Summary 28 28 31 32 34 O c e a n Surface Waves 4.1 4.2 Subject Overview Summary •. S A R O c e a n W a v e Imaging 35 35 40 41 Scattering Mechanisms 5.1.1 Specular Reflection 5.1.2 Wedge Scattering 5.1.3 Bragg Resonant Scattering 5.1.4 Bragg Scattering Nature 5.2 Ocean Wave Velocities 5.2.1 The Phase Velocity 5.2.2 The Orbital Velocity 5.3 Scene Coherence 5.3.1 Intrinsic Scale 5.3.2 Decorrelation Scale 5.4 Waves in Ice 5.5 Narrow "V" Ship Wakes 5.6 Summary 41 41 42 42 43 44 44 46 49 50 50 53 53 54 D y n a m i c Reflectivity Patterns 56 5.1 6 4 5 14 15 20 25 S A R Diffuse Target Imaging 3.1 3.2 3.3 4 Static Point Target 2.1.1 System Description 2.1.2 Applications in SAR Ocean Imaging Dynamic Point Target Partially Coherent Point Target Summary 6.1 6.2 Development Summary 56 60 iv 7 8 9 Analysis Methodology 7.1 Simulation 7.2 Statistical Testing Methodology 7.2.1 Inter-Look Cross-Correlation 7.2.2 Peak-Height Scaling 7.2.3 Azimuthal Image Bandwidth 7.3 Summary S i m u l a t e d vs. A c t u a l S A R Imagery 62 62 63 64 64 65 66 68 8.1 Specular vs. Diffuse Scattering 8.1.1 SEASAT SAR Simulation 8.1.2 SEASAT SAR Control Scene 8.1.3 SEASAT SAR Test Scene 8.2 Linear Targets and Narrow "V" Ship Wakes 8.2.1 SEASAT SAR Simulations 8.2.2 SEASAT SAR Ship Wake Imagery 8.3 Scene Coherence Time 8.4 Translating Reflectivity Patterns 68 69 69 75 79 79 80 81 87 8.5 Summary 89 Conclusions . 91 Glossary 93 References 96 A Simulation Methodology 102 A.1 Introduction A.2 Method A.2.1 The Simulation Parameter File A.2.2 Complex Impulse Response Generation A.2.3 Target Parameter Generation A.2.4 Response Summation A.2.5 Complex Output A.2.6 Detection and Real Output A.2.7 Output to Other Analysis Systems A.3 Implementation A.4 Future Extensions A. 5 Conclusions B C 102 103 103 107 110 Ill Ill Ill 112 112 112 112 S i m u l a t i o n of S A R Image P h e n o m e n a 114 B. l B.2 B.3 B. 4 114 114 118 130 Introduction Varying Point Target Density Linear Targets Conclusions A c t u a l S A R Image D a t a Sets 133 Cl SEASAT SAR C.l.l SEASAT SAR Processing C.1.2 Selected SEASAT Scenes C. 2 C-SAR C.2.1 C-SAR Processing C.2.2 Selected C-SAR Scenes C.3 Conclusions 133 133 134 135 135 135 136 v D O n e - D i m e n s i o n a l V e l o c i t y B u n c h i n g Simulations D . l Introduction D.2 Method D.3 Results D . 4 Conclusions 137 137 138 139 143 E E s t i m a t i o n of the S y s t e m Transfer F u n c t i o n E . l Introduction E.2 Method E.3 Example E.4 Conclusions 144 144 144 145 145 vi List of Tables I II III IV V VI VII VIII IX X XI XII XIII Specific forms of the derived variables 15 Summary of generic impulse response forms 26 Summary of derived SAR temporal resolutions 27 Behaviour of key image properties for point versus diffuse targets. ... 33 The coherence time scale in msec for the various cutoff wavelengths considered 52 Peak-height scaling and inter-look cross-correlation for the Goldstone subscenes 73 Peak-height scaling and inter-look cross-correlation for the Duck-X subscenes 77 Contrast measures for one- and two-dimensional C-SAR simulations. . . 85 Radar parameters relevant to the simulations 106 Scene simulation parameters ftarget density study) 115 Scene simulation parameters (linear target study) 124 One-dimensional simulation parameters for C-SAR 139 One-dimensional simulation parameters for SEASAT 141 vii List of Figures 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 31 31 31 32 32 33 33 34 35 36 37 Block diagram for the S A R azimuthal channel S A R geometry A z i m u t h a l channel block diagram for partially coherent processing. . . . Predominant point target motion parameters A z i m u t h a l channel block diagrams for single-look of a multi-look set. . . Point target Doppler history T h e orbital velocities of points along an ocean wave Definition of the angle of propagation C o n c e p t u a l velocity bunching diagram S E A S A T S A R Goldstone California scene showing extracted subscenes. Goldstone subscene 1 with the five processor focus errors Goldstone subscene 2 w i t h the five processor focus errors Goldstone subscene 5 with the five processor focus errors Image bandwidth measure for the Goldstone subscenes S E A S A T S A R D u c k - X scene showing extracted subscenes D u c k - X subscene 1 w i t h the five processor focus errors Image bandwidth measure for the D u c k - X subscenes Examples of ship wakes in S E A S A T S A R imagery C o m p a r i s o n of C - S A R w i t h simulation (coherence time) C - S A R a n d simulation image Fourier transforms Range integrated two-dimensional simulations for the three coherence times considered C o m p a r i s o n of C - S A R w i t h simulation (envelope translation) Simulation sequence flow chart Representative examples of detected impulse responses Variable target density simulation with rj = 0.00 Variable target density simulation with rj — 0.03 Derived statistics for variable target density Variable focus error linear target simulations Variable relative cross-section linear target simulations Derived statistics for linear targets in zero background Derived statistics for linear targets in a diffuse background Continued Continued Continued Variable point target density linear target simulations Continued Derived statistics for linear target of varying point target density Continued Results of one-dimensional simulations for C - S A R Degradation in contrast as a function of coherence time for C - S A R . . . . Results of one-dimensional simulations for S E A S A T A z i m u t h a l spectrum corrected for S T F viii 5 7 13 16 18 21 38 45 48 70 71 72 73 74 76 77 78 80 83 84 86 88 104 109 116 117 119 121 122 123 125 126 127 128 129 130 131 132 140 141 142 146 List of Symbols a A a(t) b(t) B c C C C d E g p / fo f f(t) g g(t) h h(t) s H s i / k K K -Kemh KB L A l(t) m M MSVR N N n(t) p PRF R(t) Ro R (T) 0 u radial component of point target acceleration (page 15) ocean wave amplitude (page 35) real, zero-mean, wide-sense stationary Gaussian random variable (page 29) real, zero-mean, wide-sense stationary Gaussian random variable (page 29) Doppler bandwidth (page 6) subscript denoting the complex impulse response (page 107) contrast function (page 143) ocean wave group velocity (page 54) ocean wave phase velocity (page 35) subscript denoting a dynamic target (page 17) ratio of peak defocus contributions of along track orbital velocity to the radial acceleration (page 49) azimuthal Doppler frequency coordinate (page 8) Doppler frequency shift (page 19) sampling frequency (page 108) scene contrast function (page 5) gravitational constant output SAR image (page 5) platform height (page 6) or water depth (page 35) compression filter (page 5) ocean wave significant waveheight (2\/2A for a sinusoidal wave) (page 39) subscript denoting the i look of a multi-look set (page 9) intensity value (page 143) th radar wavenumber (page 6) the inverse of the square of the fractional reduction in bandwidth due to application of the processing and look-extraction windows (page 11) or ocean wavenumber (page 35) azimuth F M rate (page 8) electromagnetic-hydrodynamic separation wavenumber (page 42) Bragg resonance wavenumber (page 36) the inverse of the square of the fractional reduction in bandwidth due to application of the processing filter window (page 11) look extraction filter (page 10) subscript denoting multi-look processing (page 12) the inverse of the square of the fractional reduction in processor bandwidth due to the look-extraction window (page 17) or the number of looks summed of an AT-look set (page 65) mean-squared-to-variance ratio (page 29) number of looks processed (page 14) power spectral density of a white Gaussian process (page 29) white Gaussian noise process (page 5) area of the squared look extraction window (page 12) pulse repetition frequency (page 108) target range (page 6) target range of closest approach (page 6) autocorrelation function of the process v (page 20) ix s Sv(f) t to U T Ti T T u U v VE v V V w(t) subscript denoting a static target (page 10) power spectral density function of the process v (page 20) azimuthal temporal coordinate (page 8) azimuthal temporal Doppler shift (page 17) look extraction filter centre time (page 9) available coherent integration time (page 6) look extraction filter time extent (page 9) compression filter time extent (page 9) ocean wave period (page 35) azimuthal component of point target velocity (page 15) wind speed at 19.5 m (page 39) radial component of point target velocity (page 15) velocity of the envelope of mean reflectivity (page 57) root-mean-square radial velocity component (page 22) platform velocity (page 6) ship velocity (page 54) azimuth prefilter which includes antenna pattern weighting and Doppler encoding (page 5) w (t) azimuthal antenna pattern (page 6) W Weber number (page 36) Wd image bandwidth for a diffuse target (page 30) W image bandwidth for a point target (page 31) x azimuthal spatial coordinate (page 8) a velocity bunching parameter (page 47) /? azimuthal antenna pattern beamwidth (page 6) f3(t) multiplicative phase term (page 17) 7 angle of incidence (page 37) e processor decimation factor (page 105) e processor interpolation factor (page 105) f vertical displacement of the ocean surface (page 35) rj fractional azimuthal processor quadratic phase error (page 9) rj fractional range processor quadratic phase error (page 108) fj fractional scene-induced quadratic phase error (page 17) 9[t) wide sense stationary random phase process (page 20) A radar wavelength (page 6) A ocean wavelength (page 35) fi narrow "V" ship wake half-angle (page 54) £ fractional look overlap (page 105) p azimuthal temporal (or spatial) resolution (page 11) or water density (page 35) a received power (proportional to the radar cross-section) (page 10) or surface tension (page 35) cr i relative radar cross-section (page 120) T scene coherence time scale (page 22) <t> angle between ocean wavenumber vector and platform velocity vector (page 44) (bi angle between line of targets and platform velocity vector (page 54) x{t,f) ambiguity function (page 22) •ip point target intrinsic phase (page 6) \& point target phase history (page 6) p w r m s 3 m p d u r re C x n ocean wave r a d i a n frequency (page 35) xi Acknowledgements Many individuals and organizations must be thanked for help given in the production of this thesis. My supervisory committee, consisting of Dr. W.J. Emery, Prof. P.H. LeBlond, Mr. P.R. McConnell, Mr. D.J. Okerson, and Dr. R.K. Raney, offered guidance, expertise, and encouragement throughout. Dr. R.K. Raney is particularly thanked for the patience exhibited throughout our numerous discussions and in reviewing many working papers and early drafts of this thesis. More SAR data has been examined than actually appears in this thesis. I am grateful for the opportunity, and wish to thank: Mr. J. Princz, Mr. G. Choma and their colleagues at data processing (CCRS); Dr. C. Elachi and Dr. J. Curlander (JPL); Dr. I. Cumming, Mr. P. Hasan, and Ms. P. Kavanagh (MDA); Mr. A. Haskell and Mr. M. Elkington (RAE); and the LIMEX/LEWEX experimental team. This work began at the Satellite Oceanography Laboratory at UBC when it was under the direction of Dr. W.J. Emery, and concluded at the RADARSAT Project Office, thanks to the special interest shown by Dr. N.G. Freeman. Staff at both locations must be thanked for their support and the contribution of their expertise and resources.. Particular thanks goes to Dr. R. O'Neil. Personal funding came at various times from an NSERC post-graduate scholarship, teaching assistanceships while at UBC, a DSS contract, and through F.G. Bercha & Assoc. (Ont.) Ltd.. A special thanks for support and encouragement are directed towards my wife Cathy, and our son Philippe. This work is dedicated to them. xii 1 Introduction High resolution synthetic-aperture radar (SAR) systems, such as that carried aboard SEASAT or those planned for inclusion in the ERS-1 and RADARSAT payloads, potentially provide the opportunity to measure ocean surface features on spatial and temporal scales never before realized. Such systems are not limited to oceanic research, but can also supply operational measurements on a global scale. However, the ocean imagery produced by this active microwave instrument is not well understood at this time and is subject to controversy in the relevant oceanographic and engineering communities. Problems have arisen in several areas—conflicting hypotheses and interpretations of SAR ocean data sets are frequently encountered. The following questions summarize the key areas of controversy in current SAR research into the imaging of ocean surface gravity waves: 1. Given that Bragg resonant scattering is the dominant ocean scattering mechanism, is the nature of this scattering specular or diffuse? 2. What is the size and role of the scene coherence time in the degradation of SAR resolution in the presence of surface gravity waves? 3. Is the phase velocity or the orbital velocity the relevant velocity for coherent Doppler modulations? This thesis does not purport to fully answer each of these key questions. However, each of them is addressed from three points of view:first,via careful re-derivation and extension of the relevant theory for SAR imaging of point and diffuse targets; second, by one- and two-dimensional simulations of certain image phenomena; and third, by comparison of the simulation results with actual SAR imagery. As an aid in understanding the processes involved in SAR ocean surface imaging, a two-dimensional simulation package has been developed. In the simulation procedure, 1 careful account is taken of the scene scattering nature; there are no a priori assumptions as to the statistical behaviour of each resolution element, and each point target in the scene is explicitly considered. A corresponding degree of control and flexibility is not available when studying actual SAR imagery. Chapters 2, 3, 4, 5, and 6 provide the theoretical foundations necessary to understand the SAR, its interaction with a dynamic distributed scene, and the simulation model developed for this work. Chapter 2 discusses point target imaging theory starting from the simple case of a static point target, then incorporating target dynamics, and finally target partial coherence. Results are expressed in terms of closed form expressions based upon generic Gaussian weighting functions and fundamental SAR parameters. Many of the expressions which result are new. Chapter 3 discusses diffuse target imaging theory. The diffuse target is modelled as a dense array of statistically independent point targets. Certain key statistics are derived. Differences between point and diffuse targets which may be quantitatively measured are summarized. Chapter 4 is a brief summary chapter which outlines background oceanographic information relevant to radar imaging of the ocean. Chapter 5 provides the theory for SAR imaging of surface gravity waves. Scattering mechanisms, the role of the orbital and phase velocity, and calculation of the coherence time are discussed. Two special SAR image cases are considered: ocean waves in an ice-covered sea, a case in which the high frequency image disturbing components are suppressed; and narrow "V" ship wakes, a case in which a line of Bragg targets are explicitly observed. Both of these cases allow new insights into the phenomena of interest. Chapter 6 continues the discussion of Chapter 5 but focusses upon the imaging of translating patterns of reflectivity, specifically contrasting coherent and noncoherent velocity effects upon the SAR image. New results are offered at this level of discussion. The analysis methodology is outlined in Chapter 7. Included are the role of the 2 simulations and implementation of the distinctions between point and diffuse target imagery to test for the presence of each target type. Comparisons between simulated and actual SAR imagery are performed in Chapter 8. The actual SAR data sets consist of C-band airborne SAR imagery obtained during the Labrador Extreme Waves Experiment (LEWEX) (1987) and L-band satellite SAR imagery from SEASAT (1978). Comparisons are made to study target density, scene coherence time, and noncoherent scene motions. These comparisons directly address the three areas of controversy. Details of the simulation methodologies and analysis results along with a few supporting discussions are relegated to appendices. By its inherent nature, the study undertaken in this work is multi-disciplinary. This presents numerous problems with notation and terminology; the list of symbols included in this work should be helpful. Local usage of specific symbols is explicitly indicated. Note that the radar wave parameters are given lower case symbols, while ocean wave parameters are given upper case symbols. So as not to compromise content for readers from either an oceanographic or radar background, relevant terminology is underlined in the text and denned in the Glossary, and essential features from both disciplines are reviewed in the text. This Introduction offers only a very brief assessment of the state of the art in SAR observation of the ocean surface. It has been a challenging and active research area, particularly for the past ten years. Whereas it is beyond the scope of this Introduction to describe in more detail the situation, known and controversial issues together with commentary on the relevant literature are addressed in the chapters which follow. 3 2 S A R Point Target Imaging A b s t r a c t This Chapter discusses the SAR image of an isolated point target, which is proportional to the scene/SAR /processor impulse response. The scene itself has an impact on the impulse response through resolution degradation associated with target motion and partial coherence. The SAR, that is the actual radar mechanism, imposes weighting and Doppler encoding which dictate the best possible resolution. The processor has an impact on the resolution and may impose system partial coherence. Explicit generic forms for the impulse response based upon Gaussian weighting functions are presented. The first order effects of target motion, target partial coherence, processor focus, and processor partial coherence are formulated and discussed. Results are summarized in terms of the single-look of a multi-look set and multi-look impulse responses, and in terms of the impact upon the point target resolution. Contributions consist of a unified treatment of a variety of fundamental cases and the concept of coherence time-bandwidth product as a measure of potential coherent processing support. The output of a SAR is mathematically tractable in closed form for two distinct target types:first,an isolated point target; and second, a diffuse target consisting of a dense distribution of statistically independent point targets. The SAR image of an isolated point target is the subject of this Chapter. This image is proportional to the scene/SAR/processor impulse response (or point spread function) from which we may infer fundamental performance measures such as the system resolution. We treat static, dynamic, and partially coherent point targets in this Chapter. 2.1 S t a t i c P o i n t T a r g e t A SAR is an imaging radar system which achieves a high spatial resolution in both the range (across track) and the azimuth (along track) directions [11,19,78]. High range resolution is achieved by using pulse compression techniques [35]. High azimuthal resolution is achieved by scanning the radar via the platform vehicle motion, storing the amplitude and phase structure thus imposed upon the target returns, and subsequently focussing or compressing the stored data. Azimuthal processing is analogous to range compression, but is a function of range, and is scaled very differently in time. The range scanning is at the speed of light, while the azimuthal scanning is at the speed of the carrier vehicle. This leads to a substantial decoupling of the range and azimuthal 4 Figure 1: Block diagram for the S A R azimuthal channel. dimensions for processing purposes (ignoring range migration and curvature effects). In the ensuing discussion, the azimuthal channel is implied, if not explicitly stated. The azimuthal channel is of direct relevance to the dynamic problems encountered in oceanic imaging because the time scales between the system and the scene are of the same order. These are problems to which the range channel is largely immune. 2.1.1 System Description A compact analysis for the S A R imaging process is based upon a linear convolution model [11,19,88] in which an image g{t) is formed by taking the magnitude squared of (detecting) the output of the convolution of an appropriate prefilter w(t) and processing filter h(t) over an input scene function f(t) (Fig. 1). Then, we may write si*) = 1/(0 * w(0 * MO + MO * MO I 2 . (0 where n(t) is the receiver noise. If n{t) is zero-mean white Gaussian noise, which is statistically independent of the signal, then the optimum processing filter for azimuthal compression (from a signal to noise ratio point of view) is the matchedfilter[49]: MO = ^=0 . 5 (2) The prefilter w(t) consists of two parts: phase encoding due to the platform motion past the target, and amplitude weighting due to the azimuthal antenna pattern. The phase history of a given point target is = —2fejR(t) + xj) , where k = 2ir/\ (3) is the radar wavenumber, A is the radar wavelength, and tp is the target's intrinsic phase which is introduced upon reflection. From the S A R geometry of Fig. 2, if the azimuthal antenna beamwidth f3 (radians) is small, for the fiat earth case with broadside antenna pointing, the range for a static target is given approximately by VH 2 R(t) ~ ^ + ^ ' () 4 where RQ is the target range of closest approach, and V is the platform velocity. The appropriate prefilter (up to a constant phase term) may then be written w{t) = exp j - j T r ^ t } • w {t) , (5) 2 m where w (t) is the azimuthal antenna pattern. The quadratic phase modulation (its m length, rate, and possible perturbations) is the most important aspect of the prefilter and form the focus of the work of this Chapter. The available coherent integration time T and the Doppler bandwidth B are fundamental parameters. The upper limit of coherent integration time available to the processor is determined by the antenna beamwidth. In this model, the width of w (t) is 1 m V The corresponding Doppler bandwidth is (sec). (6) (hz). (7) 2 B = ^ If dealing with a satellite SAR, the roles of footprint velocity and spacecraft velocity must be treated carefully, and used appropriately in these expressions [65]. Replace V with the footprint velocity Vj in equation (6). 1 If dealing with a satellite SAR, replace V with the spacecraft velocity V / [65] in equation (7). Vf and V / could differ by 15% for a typical low Earth orbit SAR. 2 s s C 6 C Figure 2: S A R geometry. The quadratic phase of the azimuthal prefilter leads to a linear frequency modulated encoding of FM rate K a - T - \ R 0 - ( 8 ) For practical SAR systems, the time-bandwidth product TB *M- = A (9) satisfies TB ~> 1, a property which allows significant simplifications in analysis through use of the principle of stationary phase [36]. The time-bandwidth product represents the improvement in resolution which may be afforded through coherent processing. (Impulse response-type functions generally satisfy TB ~ 1 suggesting that such signals will not support coherent processing for improved resolution.) In the formulations in this work, the azimuthal spatial coordinate x is related to the azimuthal time coordinate t, and to the azimuthal Doppler frequency coordinate / (if TB » 1) via x = Vt = V^-, (10) for a static target. This property sets the static impulse response apart from the dynamic impulse response where the azimuthal coordinate may be modified by the target's motion parameters. This is of central interest in this work. Equation (10) shows that the azimuthal time coordinate is proportional to the Doppler frequency, a property specific to large time-bandwidth product signals. It turns out that this is a contributing factor to some of the controversy in the SAR ocean wave imaging problem. Appreciating this relationship helps to resolve an apparent paradox (see Chapter 6). It is convenient to introduce an explicit generic form for the azimuthal antenna pattern. A convenient choice which represents the antenna pattern's main lobe reasonably well, and allows simple closed form expressions to be found for the SAR/processor 8 impulse response and resolution, is Gaussian amplitude weighting : 3 (11) The time extent T of this form may be calculated using an equivalent rectangular width norm. 4 The processing filter, though nominally matched, may contain departures from its ideal form. There is a potential mismatch between the prefilter and processing filter azimuthal F M rates, referred to as a focus error or a quadratic phase error. Such a condition could arise intentionally or accidentally, and always occurs to some degree near segment edges in digital S A R processors which operate in a range-segment mode [9]. Another departure could be in the actual time (hence frequency) weighting of the processing filter. In keeping with our use of generic Gaussian weighting functions, the processing filter may be written (12) where T is the time extent of the processing filter and rj is the fractional error in the p processor F M rate . 5 The processing filter is frequently operated in a multi-look mode. That is, N portions of the Doppler spectrum, say those centred at times 3 and of time extent T/, It should be noted that some authors [4,77] use to define the azimuthal antenna pattern. This form and equation (11) are identical if T — y/2/n T. a The Doppler bandwidth should also be correspondingly scaled when interpreting results derived from these sources. The equivalent rectangle is the width of the rectangle of height equal to the maximum value and containing the same area. That is 4 For a Gaussian function, this resolution measure is related to the full-width-at-half-max ( F W H M ) by PFWHM = 0.94p, to the radius-of-gyration by p = 0.80/9, and to the standard deviation by a = 0.40/9. rg The fractional processor F M rate error is occasionally written as [l + rj') rather than (1 + rj). Then, in the event of a small focus error, these formulations are related by n = 2T;'. 5 2 9 are coherently processed, detected, and summed to produce the S A R image. Such a partially coherent processing scheme uses ensemble averaging to trade speckle smoothing (radiometric fidelity) for spatial resolution (see Chapter 3). Each portion of the Doppler spectrum so processed is called a look, and this procedure is termed multi-look processing. A n individual look extraction filter may be written M*)=expj-*fi-=^} - (13) Since we are dealing with a large time-bandwidth product signal, windowing in the frequency domain is equivalent to windowing in the time domain. Thus, the effective processing filter may be written h (t)=h(t)-l (t) i (14) i for the i look. This equality only applies prior to application of the processing filter th when the time-bandwidth product is large. Ignoring the receiver noise, the SAR/processor response to the static, impulsive target f(t) = y/a6(t) is g.{t,ti)=o\w{t)*hi{t)\ 2 , (15) where w(t) is defined in equations (5) and (11), and h{(t) is defined in equations (12), (13), and (14). The scaling constant o is the power received at the radar which is proportional to the impulsive target's radar cross-section, accounting for geometry effects in microwave propagation and reflection through the radar equation [81]. For convenience, o is set equal to 1 in the ensuing discussion. The required convolution is conveniently formed using Fourier transform techniques. The forward transforms of w(t) and hi[t) take advantage of these signal's large timebandwidth product and use the principle of stationary phase [19,63]. The transformed functions are then multiplied together and inverse Fourier transformed. For n <C 1, we 10 find that T g ° > [ t t i ] = B L tn 2 • \- *K f?) T w r r i exp x \ ~ *K+ 2 2 • e x p (16) where c = TBr). The fractional reduction in bandwidth due to the processing filter window is 1/VL where L = 1 + (T/T ) . 2 P The fractional reduction in bandwidth due to the processing filter window and the look extraction window is 1/y/K where K= 1 + [T/Tp) 2 + (T/Ti) . 2 This response (and all other generic forms derived) are listed in Table II on page 26. The term (T/Tj)£,-/K in the second exponentiation of equation (16) is the centre 2 of the effective data window consisting of the pre-filter, processing filter, and offset look-extraction windows. This term interacts with a focus error to govern the relative azimuthal locations of the impulse response associated with each look. That is, the response in the individual looks becomes misregistered in azimuth as well as broadened in response to the focus error. This effect turns out to be of importance in resolving the third key controversial question posed in the Introduction, which is discussed in Chapter 6. The misregistration has a further impact on the resolution after look summation. Some digital SAR processors utilize such azimuthal offsets to estimate the degree of processor focus error. Such a procedure has been termed "autofocus" [39]. The static single-look impulse response has an (equivalent rectangular norm) azimuthal resolution (or time duration) of K e K> ~ BXl 2~ 2 Pai 1+ (17) (The azimuthal spatial resolution is p V.) This resolution (and all other generic forms gi derived) are listed in Table III on page 27. A multi-look (or partially coherent) digital SAR processor is normally implemented as a discrete frequency-plane mixed-integrator; only looks from several discrete-look locations are extracted and summed to produce the final SAR image [9]. The block 11 diagram for this process is shown in Fig. 3 and may be represented by 9*m{t) =E () 18 i where the subscript m signifies multi-look summation. A partially coherent optical S A R processor can operate as a continuous or scanning frequency-plane mixed integrator in that looks extracted from a continuum of frequency locations are summed to produce the final S A R image [91]. This process may be represented by 9 (t) sm 1 r°° = - P 9,{t,ti)dti, (19) J-oo where p is a normalization constant which is equal to the area of the squared lookextraction window. It has been shown that discrete and continuous modes of operation produce substantially equivalent results as far as resolution degradation and speckle smoothing are concerned, although some improvement may be enjoyed by overlapping the looks [91]. Either approach is satisfactory for the work of this research, although closed form expressions are more readily obtained with equation (19). Consider now the multi-look case. Using equations (16) and (19) with p = Ti/y/2 we find that the multi-look impulse response is ' - w = v m ^ - ^ h ' i ^ i L ^ • (20) which has resolution e 2 - 1 1 - + • irr • Psm — ~B\^ 2 (21) Observe that in the event of a focus error the resolution of the multi-look response is degraded more than that of a single look of that same multi-look set. Care is required in using published results on focus sensitivity as the single-look of a multi-look set case or the multi-look case is seldom clearly stated. The susceptibility of SAR/processor resolution to a focus error may be judged by observing when the focus error term becomes of comparable size to the perfect focus 12 Figure 3: Azimuthal channel block diagram for partially coherent processing. 13 resolution. For a single-look of a multi-look set, this occurs when , , \ r K > \ > T B 2iV = 2 T W W For the multi-look case, this occurs when I "I TF = TB> (23) In each case, the second equality applies to matched filter processing. These results demonstrate that a focus error is more limiting in the multi-look case—the normal processor operating configuration. This is largely due to the responses within the individual looks becoming misregistered. A n alternative measure of the susceptibility to focus errors is the depth of focus [81], which is the nominal range times the right hand side of the appropriate expression in equation (22) or (23). This represents the error in target range which still allows reasonable focussing properties. For this alternate definition, the same sensitivities result for the single-look versus multi-look cases. 2.1.2 Applications in S A R Ocean Imaging Two specific processing models are generally found in the S A R ocean literature. The first model is the case of matched filter processing (T = T) which is what all practical p S A R processors attempt to approach. The second model is infinite rectangular filter processing (T — • oo) which is an unrealizable processing scheme, although somewhat p simpler to handle analytically. It is well known that a single-look of a multi-look set has resolution degraded in proportion to the number of looks extracted N for matched filter processing [19]. Using this result, the derived terms in the impulse response equations developed in this Chapter take the forms shown in Table I in terms of B, T, and N. matched filter processing (T = T, rj = 0), we find that p p intensity is T /2N. 2 s m = N/B and that the peak These are well known results in S A R analysis. 14 Note that for Table I: Specific forms of the derived variables. Cases of matched filter processing and infinite rectangular filter processing in terms of the fundamental parameters B, T, and N. Matched Filter Infinite Rectangular Generic Form rp oo T =T P N 2N 2 2 1 2 M = i + (a) ! 2N -1 2 1 2(N -1) 2 Pam peak intensity N B N 2N w It is helpful to have at hand Table I, since many results available in the literature are based upon one or the other assumptions for the filter weighting, although not always clearly. Since resolution of certain areas of controversy depends upon factors of two, the contents of the Table may be useful. 2.2 D y n a m i c P o i n t T a r g e t In the event of point target motion, the important target velocity parameters are the radial velocity v (the target velocity component away from the radar along the radial line of sight), the radial acceleration a, and the along-track velocity u [60] as shown in Fig. 4. (If the motion parameters change during observation, the average values should be used. The appropriate average is over the duration of the time domain lookextraction window.) To second order, the moving point target range equation may be written R(t) « Ro + vt + 2u 2Ro 1 ~T Ra 0 + T 2 (24) The target's dynamic properties may be represented by a multiplicative quadratic 15 Figure 4: Predominant point target motion paramet 16 phase term [12,13,51,61,63]. The resulting system block diagram is shown in Fig. 5b in which the multiplicative phase term is B (25) [t t + fjt ] , 2 D where tr, = RQV/V is proportional to the Doppler shift due to the radial velocity 2 component and where fj = Roa/V —2u/V is the contribution of the target motion to the 2 focus error . We find that the dynamic impulse response of this scene/SAR/processor 6 system is exp < — 27T 2TT exp 9d{t,U) = K + e /K 2 KT 2 t + t where Sn = n — fj and now e = TBSn. D 1 Lt\ + 2tit + D 6rj K) 2 D 6n K' (26) The fractional reduction in processor bandwidth due to the look-extraction window is l/y/M where M The +' \TJ Mt = 1+ (Ti/T ) . 2 p resolution (width of gd{t,ti)) is Pdi - K 77 ~2 1 + K 2 (27) the same as equation (17), the result previously derived for the static case. It is possible to compensate for resolution loss due to coherent target motion by matching rj, the processor focus parameter, to fj for a given target. This would be at the expense of all other targets not subject to the same motion parameters becoming correspondingly defocussed. Thus, an ensemble of point targets, each with different motion parameters, cannot all be matched simultaneously by the processor. Note also that attempts to It has become customary in the SAR ocean literature to express the focus error in terms of an equivalent azimuthal velocity Uf. Then, we may write 6 V , expresses the equivalent target velocity necessary to account for the net focus perturbation. This practice has lead to rather misleading interpretations of SAR image phenomena, since it suggests that the targets themselves are moving at such velocities. This is patently not the case from an oceanographic point of. view (see Chapter 6). 17 Figure 5: Azimuthal channel block diagrams for single-look of a multi-look set. a) Static point target, b) dynamic point target, c) and d) partially coherent point target. Representations c) and d) are equivalent if {6} is wide sense stationary in temporal variables. 18 estimate the focus error of moving targets by comparison of the look offsets produces incorrect results. Likewise, image improvement through focus adjustment ("autofo- cus") schemes should not be used in conjunction with moving target fields. These points are readdressed in subsequent sections, and cannot be overemphasized given the current state of confusion which is found in the literature [20,25,26,27,50,72,73]. The peak intensity of each look may be severely attenuated by the radial Doppler term io if the target Doppler shift is an appreciable fraction of the system Doppler bandwidth. This effect is called Doppler suppression and may be important for an airborne S A R observing the ocean surface. The Doppler suppression may be slightly reduced by increasing the width of the processing filter T . However, this is at the p expense of an increased signal-to-noise ratio, and could lead to azimuthal ambiguity problems [19]. In the event of multi-look summation, equation (19) gives >2 9 d m { t ) = 7 m r T e x p { - 2 " M ^ P , (28) - "KT?JL 2 which has resolution _ 1 This resolution is the same as equation (21), the result previously derived for the multi-look static case. The peak of the response is reduced by Doppler suppression (first exponentiation in equation (28)). The azimuthal location of a given target is governed by the time at which zero Doppler occurs. Thus, the target location is offset in azimuth due to the Doppler frequency shift fu = 2v/X in accordance with the to term in the second exponentiation in equation (28). The spatial offset in response to a radial Doppler shift term (in the absence of focus error) is A x = _ T V f D = „ 19 V t D = , (30) which is a well known result. Additional position modulations result from interaction with the focus error terms as indicated in equation (28) and as shown in Fig. 6. 2.3 P a r t i a l l y C o h e r e n t P o i n t T a r g e t If target motion is not adequately described by first and second order terms as in equation (24), or the target is subject to a multiplicative complex fading process, an additional phase perturbation parameter 6(t) should be introduced into the target's phase history (Fig. 5c). Let {0{t)} be a sample function of a random process, which for convenience may be taken to be zero-mean, wide-sense stationary. This parameter models random phase behaviour for the ensemble of targets which represents a partially coherent scene. The expected scene/SAR/processor impulse response for a single-look may be written E{g{t,ti)} = E{\w(t) • exp[j0(*)] * /i (i)| } , (31) 2 t where E{ } is the ensemble average operator, which is equivalent to a time average for an ergodic wide sense stationary random process. We may write E{g(t,U)} = fj 7 w{t )v4t2)h {t - t )hi[t - t )R (t - h) dh dt , l i x 2 a 2 2 (32) where (33) R (r) = E {exp[j0{t)]exp[-j6(t + r)]} e is the autocorrelation function of exp[j0[t)]. Then, we find that oo \w{t) exp[j2nvt] * h^t) | S {u) du , (34) 2 t / -oo If the look-extraction windows are explicitly included at this stage, the order of summation and integration may be exchanged to give 7 E{g{t)} = jj w(h)w(t )h{t - t )h(t - t )Re(t - t!)/(t - t ) dh dt 2 x 2 2 2 x 2 where ifo-h) = ^2l {t-t )l {t-t ) i 1 i 2 is a correlation function. This relates the development considered here to the quadratic filter theory formulation [61,63,64]. Unfortunately, the elegant quadratic filter formulation imposes increased mathematical complexity upon some of the issues considered in this work, and does not allow explicit exposure of certain parameter dependencies that are of central interest. 20 Figure 6: Point target Doppler history. The time of zero Doppler, hence the target location, is governed by the Doppler shift fn = 2v/X and the target motion induced focus error fj. Three cases are represented: a) }D = 0, fj = 0; b) f > 0, fj = 0; and c) f > 0, fj > 0. D D 21 where S$(f), the power spectrum of the random process exp[jd(t)], is the Fourier transform of R$(T). The expected impulse response may be written oo \Xv,h{t,v)\ S {v)du , (35) 2 / -oo e where Xwh{f,v) is the (non-normalized) cross-ambiguity function [70,88] for the SAR azimuthal channel. Thus, for a static target we may write oo 9d{t,U) S {v) du , / -oo t =Tv/B,fj=0 (36) e D and for a dynamic target oo / -oo 9d{t,ti) Sg (u) du . t ->t +Tu/B D (37) D It is convenient (and adequate for present purposes) to assume a normal correlation function for the random process exp[j'0(t)]. Then, we may write ^2\ ( R {r) 8 exp { - T r — \ , = (38) where r is the coherence time scale. This time may be denned on the basis of the c target's rms velocity as (39) 2Vr This formulation will be examined in greater detail in subsequent sections. The dynamic single-look impulse response in the face of a finite scene coherence time under the above assumptions is T E{g {t,t )} = d f 2TT exp I - - — 2 [LQt + 2Rtit " " l KT? 2 i 2 ^PK + PeyK + 2{T/T )* C ( •exp \^ —2?r-K B + e /K 2 + + MRt ] 2 D D 2 {2/P){T/r ) 2 c (40) 22 where 2M ( T \ ~ K \ T + 2 Q 1 2 ) t (41) p F ' T\ (I 2 KLP \TiJ (42) (T B) ' 2 C 2M / T \ R (TJ 1 (43) {T B) ' 2 C and the resolution is if Pi = 1+ B\ 2 K +' KP 2 Ty (44) \TJ We now introduce a new measure of fine-scale scene dynamics which is the concept of a coherence time-bandwidth product r B. The magnitude of this product seems to C differentiate between two distinct asymptotic operating regimes. If TB » 1 , (45) C we are in either a fully or partially coherent regime. In this case P,Q,R = 1, and if tr, = 0 we get the simplified and important result [61,63] B „ (B E{g{t,t )} = g(t,t )*-S (-t) i i (46) e In this case, the mean effects of a finite scene coherence time may be represented by a low-pass filter applied after detection. This case is shown in Fig. 5d. Thus, Fig. 5c and Fig. 5d are equivalent representations if inequality (45) is satisfied. Equation (46) directly applies to the expected response of an isolated point target only. In this case, the expected single-look static impulse response is E{g (t,U)} s ^ 2 + 2 e exp < — 27T 2 KT + (2/*r)(r/r ) c B Lt 2 T 2 2 exp < — 27r K + e /K + 2 ( r / r ) 2 c 2 [ VT,/ K \ (47) which has resolution K e 2 Pi = TT B\ 2 23 2 / T s 2 (48) After multi-look summation, the impulse response becomes T E{g {t)} sm 2 22 j- B t exp < —2?r jKL + e* + K + e*/L + (2/K)(T/T y c 2{T/r y (49) c which has resolution K B\ 2 1 + 2 (T +K KL 21 (50) x, , c These results become identical to those of the static impulse response in the limit of T/T — 0. C The condition of equation (45) is likely well satisfied for all satellite S A R platforms. In this case, B is O(10 ) Hz so that even a coherence time of 10 msec satisfies the 3 condition. However, in the case of an airborne S A R platform, B may be as small as O(50) Hz, and equation (45) may not hold. Then, the impulse response described by equation (46) is not valid. Thus, other operating regimes should be considered. A second asymptotic operating regime is one which is characterized by (51) TB < 1. C This represents the noncoherent limit. In this case, a S A R is expected to approach the behaviour of a real aperture radar. The expected single-look impulse response is E{g{t,ti)} = s/2K M exp < — 2TT exp < —27r Srj - 1 M 1 2' (52) which has resolution Pt = T K l 2M (53) After multi-look summation, we find that E{g {t)} = exp m (54) which has resolution Pm = IL T l- . 24 p] (55) As the target coherence time decreases, the resolution becomes limited by the breadth of either the processing or look-extraction window. Although not of direct relevance to this work, these results illustrate some of the subtleties in passing continuously from a coherent or partially coherent limit to a noncoherent limit [18]. Note that p m — • oo for infinite rectangular filter processing. The level of coherence being governed by the coherence time-bandwidth product is rather surprising. In fact, it might be expected that the coherence should be governed by the portion of the bandwidth relevant to the coherence time (i.e. K T ), a c in analogy to the case of a focus error. However, for both the single-look of a multi-look set case and the multi-look case, the impact upon the resolution by the coherence time term is the same and depends upon the entire available coherent integration time T. This is in contrast to the case of a focus error for which the resolution degradation due to the focus error depends upon the specific processing scheme. Thus, the coherence timebandwidth product is a measure of the degree of coherent, or of multi-look processing, which the signal will support. This restraint is a new result. 2.4 S u m m a r y This Chapter has presented reasonably simple, closed-form expressions for the scene/SAR/processor impulse response. The calculations were based upon generic Gaussian weighting functions and considered the effects of target motion, target partial coherence, and various processor configurations. It was suggested that the scene coherence time-bandwidth product represents a measure of the degree of coherent processing which the received signal from a given point target will support. The general operating regime is either fully or partially coherent. The derivations of this Chapter are summarized in Table II for the impulse responses, and in Table III for the corresponding resolution measures. 25 Table II: Summary of generic impulse response forms. Single-look of a multi-look set and multi-look processing for the cases of static, dynamic, and partially coherent targets are included. STATIC POINT T A R G E T single-look of a multi-look set *(*•*«) = 7^7? e X p{ 2 ^^-}- e X pJ 2 *K + (S/K) [* + ( £ Kl] ) J ( ) 16 multi-look DYNAMIC POINT T A R G E T single-look of a = vwh? multi-look set e x multi-look p|-2^F^w[ t • E X {-$rj i i Lt P + 2 i*D + Mt%]} +^( -7r) + (^) ^] 1 2 W = v / i f e ^ • P { *f$*} • P { *K&fL e x 2t e x [* + 2 | ( ) (l £ ) ] 26 } (28) PARTIALLY C O H E R E N T POINT T A R G E T single-look of a multi-look set • pJ e X multi-look g{ . (t)}g m ^ K V K L + e 3 + T = — * = ^ ^/K+2(T/r^ + [ 3 2 / i f ) ( r / r e ) 3 (sec) (H.) 26 [* + exp{ ( T , ) J ^ ] J ( ) 47 2 . ^ 1 ^ ^ } (49) Table III: Summary of derived S A R temporal resolutions. The generic resolution and resolution for the case of matched filter processing are included. single-look of multi-look multi-look set static or dynamic ft = ^ f [ l + p ] m „ partially coherent + * fe)'] * = 7^/f * = £\/ 1 + ^ ^ t e ) + T = ~Y~ * = ^ (48) 2 ^ec) (H.) 27 (21),(29) p = ^f[i+^r] (17),(27) — N + £ fe)'] *» = ^ / f = F^/ /l _i_ 1 + 477^ + TT" 5 fe) (50) 3 S A R Diffuse Target Imaging A b s t r a c t This Chapter discusses SAR imagery of diffuse targets. A diffuse target consists of many point targets per resolution element, and leads to the granular speckle image pattern associated with coherent systems. The second order image statistics of such patterns are calculated. It is shown that image resolution may be traded directly for radiometric fidelity, and that the speckle covariance function is of the same width as the static impulse response. Key issues are summarized which may be used to distinguish imagery of point targets from imagery of diffuse targets. These issues are the correlation between looks, peak-height scaling in response to multi-look processing, and the response of the image spectral width to a processor focus error. This Chapter ties the little known diffuse target analysis into the point target analysis framework developed in Chapter 2. A diffuse target represents the extreme case of many independent point targets per resolution cell. The SAR response to such a target is distinctly different from the case of the isolated point target which was treated in the previous Chapter. Here, we examine the S A R response to a diffuse target through the second order image statistics, and then, contrast the results with those of the isolated point target case. 3.1 S e c o n d O r d e r I m a g e S t a t i s t i c s A diffuse scene may be regarded as being composed of an ensemble of point targets, each characterized by a (statistically independent) radar cross-section and intrinsic phase [69,71]. Then, for each resolution element, the received complex amplitude is the vector summation of the contributions from each point target. If the number of point targets per resolution element is large (a classical guideline is > 5 [34]), the statistically independent intrinsic phases lead to a two-dimensional random walk and Rayleigh statistics for each image element [71]. The resulting grainy multiplicative noise, known as speckle, is characteristic of all coherent imaging systems. 3.1.1 Static Diffuse Target Consider the case of infinitely many independent, static, point targets per resolution element. For such a diffuse scene, the input to the azimuthal channel may be modelled 28 as f(t)=a(t)+jb{t) , (56) where a(t) and b(t) are real, zero-mean, statistically independent, Gaussian random variables with identical variances. If f(t) is input to a magnitude squared detector, that is 9(t) = \f(t)\ , 2 (57) then the second-order statistics of the output are [52] R {r) = 4R {0)+4R {r) (58) = 4R (0)6(f) + 4S(f) * S(f) , (59) 2 2 g and S (f) g 2 where R{T) = E{a{t)a{t + r)} = E{b{t)b{t + T)} = N 6{T) (60) 0 and S(f) = N 0 . (61) In a S A R system, the input f(t) passes through the linear filters w(t) and /i,-(i) prior to the detector. The input spectrum to the detector is modified to 5,(/) = | i y ( / ) 7 J ( / ) | 5 ( / ) , (62) 2 t where the subscript i refers to the i th look-extraction filter. Equation (62) may be substituted directly into equations (58) and (59) to give the image spectrum as S , ( / ) = 4 [ j V / |W(/)7J,(/)| df]\4N \W(f)H (f)\ 2 2 2 0 * |^(/)-ff,(/)| i 2 (63) . A norm for radar image radiometric resolution, the expected variation of the radar cross-section for each pixel, is the mean-squared-to-variance ratio (MSVR). Noting that the mean and variance add linearly after look summation, we find that [H\W(fmf)\ f f MSVR = 2 d dti M\W(\)H (\)\ \W(f-\)H (f-\)\ dfd\dt 2 i 2 i 29 i ' v ' Substituting the generic forms for W(f) and of Chapter 2 gives (65) MSVR = where the second equality applies to matched filter processing . The proportionality 8 to N, the number of looks processed, is a well known result [10,58] and indicates that resolution may be traded directly for radiometric fidelity through the ensemble averaging process of multi-look addition. This is possible since there is very little correlation between the speckle patterns in looks derived from non-overlapping portions of the Doppler spectrum. As the number of independent looks increases, the radiometric fluctuations decrease in proportion to l/\/~N. (The factor of A/2 in equation (65) arises because the scanning frequency-plane mixed integrator has an infinite scanning bandwidth, and yields some modest improvement over the discrete mixed integrator [91].) Note that radiometric fidelity is not impaired by multiplicative phase errors since the input signal phase is already a wide bandwidth random process. In other words, scene partial coherence cannot add more independence from look-to-look [64,65]. The azimuthal image spectrum for the diffuse input of equation (56), ignoring a real proportionality constant and a delta function at the origin, is (66) where the equality applies to matched filter processing. This image spectrum has (equivalent rectangular width norm) bandwidth W d = B B y/K (67) V2N ' where the second equality applies to matched filter processing. A n important consequence is that the image bandwidth of a diffuse scene is independent of the processor focus setting. This result may be contrasted with the azimuthal spectral response to the 8 This is the SAR version of a classic result, known since the 19 century [69]. th 30 image of a static point target (see Table II). Ignoring a real proportionality constant, the azimuthal image spectrum of an isolated, static point target is (68) where the equality applies to matched filter processing. This spectrum has bandwidth (69) where the second equality again applies to matched filter processing. In contrast to the diffuse case, the bandwidth measure for the point target is dependent upon the processor focus setting. We see also that the speckle covariance function is of the same shape as the zero-focus-error, squared impulse response . This result is used in 9 correcting for the system transfer function in ocean wave spectra derived from S A R imagery [5] as discussed in Appendix E . 3.1.2 Dynamic Diffuse Target In the event of a dynamic diffuse target, the situation becomes somewhat more complicated [23,61,63]. Each point target which constitutes the scene could be subject to different sets of motion parameters, resulting in each target having a different impulse response. The Doppler shift term t D causes each target to be shifted differentially in azimuth (and possibly leads to differing Doppler suppression for each such target). If the scene dynamics are random, the target redistribution is also random, and the resulting image will be indistinguishable from the case of a static diffuse target. However, if the scene dynamics are well structured, as can be the case for the ocean surface, otherwise uniformly distributed targets may have their images clustered producing a This assertion applies to the restraint of Gaussian weighting functions. equivalence applies more generally [63]. 9 31 An upper bound width nonuniform contrast in the S A R image [37]. This effect is termed velocity bunching and is an important mechanism for S A R imaging of the ocean surface (see Chapter 5). The scene-induced focus error fj causes the impulse response associated with each target to be broadened independently of the response of the other targets. 3.2 P o i n t v s . D i f f u s e S c a t t e r i n g It is useful to be able to differentiate between cases of point and diffuse targets in S A R imagery. As has been previously observed [19,63,81], three fundamental distinctions between S A R imagery of point and diffuse targets are apparent. It has been proposed [64] that these distinctions may be used to test control scenes derived from a given S A R and processor which include known examples of point and diffuse targets. These results may be applied to a test scene to establish the nature of the targets present. The first fundamental distinction is the correlation between looks in a multi-look set. For an isolated point target, the correlation between looks is large, particularly in the event of no processor focus error. As a focus error is introduced, the peak of the correlation function will remain large, but the responses in the individual looks will become misregistered. This fact is used in the auto-focus technique. For a diffuse target, however, the correlation between looks derived from nominally non-overlapping portions of the Doppler spectrum will be low. It is this property which allows improvement in radiometric resolution through multi-look processing. The second distinction is the peak scaling in the face of multi-look processing. From Table I on page 14, the peak intensity scales as 1/iV for point targets. If the processor is a discrete frequency-plane mixed-integrator, a large focus error could misregister the looks enough to cause this property to break down. The standard deviation of a diffuse scene is a measure of the peak heights, and in equation (65) this scaling is shown to be I/VN. 32 Table IV: Behaviour of key image properties for point versus diffuse targets property point diffuse correlation between looks high Peak intensity dependence on N focus sensitivity azimuthal spectral width a low 1 1 N s/N yes no (resolution) -1 bandwidth measure of: The third key difference is the response to a focus error. This response may be judged by measuring the image bandwidth. We have seen in equation (67) that in the event of a diffuse scene, the bandwidth will not respond to a processor focus perturbation. However, as shown in equation (69), the impulse response will broaden and the bandwidth will decrease as the processor focus error increases. These contrasting image properties are summarized in Table IV and are used as a basis for a statistical testing methodology which is found in Chapter 7, and is intended to address the issue of point versus diffuse scattering in S A R ocean surface imaging. Appropriate and tractable statistical description is awkward for the cases in which the average number of scatterers per resolution element is neither very small (purely isolated point targets) nor very large (purely diffuse target). A modified Bessel-function distribution has been proposed for such cases [29,30] and has been proven successful in describing many natural phenomena [32] including non-Rayleigh radar returns [31]. Unfortunately, a physical basis has not been found for the applicability of these Kdistributions . 10 T h e problem of detecting a specular scatterer in clutter is representative of the difficulties encountered [75,76], although this problem does not directly apply here. 1 0 33 3.3 S u m m a r y T h i s Chapter has cast the imaging of diffuse targets into the framework presented in Chapter 2. T h e response to a diffuse target is formulated in terms of the second order image statistics. These results are contrasted with those of an isolated point target. T h e key areas of departure, outlined in Table I V , are the inter-look cross-correlation, the peak-height scaling in response to multi-look processing, and the focus sensitivity which is manifested in the image spectral bandwidth. 34 4 Ocean Surface Waves A b s t r a c t This Chapter presents a succinct review of ocean surface waves which is based upon currently accepted oceanographic theory. The content is neither new nor controversial. However, many of the concepts presented become misconstrued when interpreted within the framework of SAR ocean surface imaging. Inclusion in this thesis is essential for an interdisciplinary understanding of the issues involved. The behaviour, analysis, and observation of ocean surface waves are well developed subjects which can be found in many introductory level [57] and advanced level oceanographic text books [36,38,48,53]. The content of this Chapter is drawn primarily from such sources. 4.1 S u b j e c t O v e r v i e w We begin by considering the ocean surface displacement in response to the monochromatic sinusoidal ocean surface wave given by f(x,r) = A cos(Kx - nt) , (70) where A is the wave amplitude, K is the wavenumber. and Cl is the wave frequency. This formulation represents a travelling train of waves of infinite extent with wavelength A = 2TV/K and period T w = 27r/Jl. A point of constant phase on the surface of the wave advances with a speed given by the magnitude of the phase velocity, which is n C = % • P ( ) 71 The phase velocity is in the direction of wave propagation. A n important characteristic of ocean surface waves is their dispersive nature. The dispersion relation, based upon a linearized analysis, is given by n 2 = [gK + j tanh(iC/i) , (72) where h is the water depth, a is the surface tension (~ 0.074 N / m ) , and p is the water density (~ 1000kg/m ). This result is independent of the wave amplitude, and holds 3 if A <C A. 35 In the open ocean, surface gravity waves often satisfy h ^> A. Then, Kh >• 1 and the tanh(i C/i) term approaches unity. This is the deep water or short wave limit. i The role of surface tension as a restoring force is quantified by the oK /p z term in the dispersion relation, while the role of gravity as a restoring force is quantified by the gK term. The importance of these two terms may be judged by their ratio. This number, called the Weber number, is given by W = oK 2 9P (73) The critical value occurs when the two restoring forces are of equal importance, in which case W = 1. For typical oceanic conditions this occurs for A = 1.73 cm. If A <C 1.73 cm, surface tension dominates as the restoring force. Such waves are called capillary waves. If A » 1.73 cm, gravity dominates as the restoring force. Such waves are called surface gravity waves. For surface gravity waves (74) illustrating that the longer waves will travel the fastest. As surface gravity waves propagate away from their point of origin, the various wave components separate from each other, with the longer waves moving in advance of the shorter waves. The wave field, at any time and some distance away from the generation area, becomes more purely of one wavenumber, assuming there are no locally generated waves. Such wave systems are called swell and form the basis of assumptions for SAR ocean observation. In contrast, the relatively confused seastates found at the areas of generation are often called wind- waves. It is well known that microwave scattering from rough surfaces, such as the ocean surface, is dominated by the scattering surface Fourier component whose wavenumber is given by [89,90] K B = 2k sin 7 , 36 (75) where k is the radar wavenumber and 7 is the angle of incidence (measured from vertical). This scattering mechanism is called Bragg scattering. The first-order Bragg wavelength is equal to the radar half-wavelength projected onto the ocean surface. For C-band and L-band radars, this wavelength falls into the short gravity wave regime, being typically 6 cm and 3 0 cm respectively. Thus, a microwave ocean imaging system, such as a S A R , images ocean surface features indirectly; the ocean is seen via Bragg scattering from short gravity waves which are subject to modulations by the longer scale waves (and features) which appear in the S A R imagery. The modulations include velocity changes, straining, and tilting. There is an important distinction between the velocity of a point of constant phase on the wave surface and the velocity of a packet of water on the ocean surface. While the wave surface advances at the phase speed C p in the direction of propagation, a packet of water at the ocean surface will be comparatively stationary, yet subject to vertical and horizontal oscillatory velocity components given by u(x,t) = VIA cos{Kx - Hi) , (76) and w(x, t) = ClA sin(Kx - fit) (77) in deep water. It is apparent that each such water packet is confined to a local domain and describes a circle rather than translating at a large constant rate such as the wave phase velocity. Thus, as a wave passes by a given point on the ocean surface advancing with phase speed C , the ocean surface heaves up-and-down and laterally, p causing the individual particles on the surface to describe a circle. This orbital motion is illustrated in Fig. 7. Any short waves riding upon the long waves will have their velocities modulated by the orbital velocity rather than the phase velocity. The straining of short gravity waves riding upon long gravity waves is formulated by appealing to conservation of wave action density and using a WKB-type of perturbation 37 38 approximation (so that the deep water dispersion relation applies locally to the short waves) [16,53,54]. The effect may be formulated as arising from the currents set up on the ocean surface due to the passage of the wave. As shown by Fig. 7, there is an area of convergence along the wave's leading face, and an area of divergence along the wave's trailing face. It has been proven that the short waves are shortest near the long wave crests and longest near the long wave troughs, that the short waves tend to align with the direction of long wave propagation near the long crests, and that the short waves become higher and steeper near the long wave crests. Thus, the wave crests will be rough relative to the wave troughs, and the small scale structure is dependent upon the large scale wave position. A sinusoidal ocean wave description can be quite inadequate. Instead, the ocean surface is often described in terms of its one-dimensional frequency spectrum ^(fi). This spectrum is defined by roo j o A% s (n)dn = {^) = i T , (78) where the second equality applies to the sinusoidal case discussed above. The significant waveheight is defined by B. = 4 <?> = 2v 2A . / (79) / V Useful generic forms for S ( are available [22,55]. The Pierson-Moskowitz wind-wave spectrum is given by ^)-(S-p{-"(^) }4 <»> where a = 0.0081 and j3 = 0.74 are dimensionless empirical constants, and U is the windspeed at 19.5 m above the surface. The windspeed is related to the significant waveheight (in M K S units) via H = 2.12(10)- t/ 2 3 39 2 (81) in the Pierson-Moskowitz model. This spectral form does not include fetch dependence or wind duration dependence. However, it is illustrative of the richness of the spectral description of a fully developed sea. 4.2 S u m m a r y In this Chapter, we have reviewed some basic oceanographic theory in preparation for the discussions of S A R ocean surface imaging which are found in Chapter 5. A n important point is the distinction between the phase and the orbital velocities. These velocities are subject to misleading interpretations within the framework of SAR ocean imaging because of the different ways in which they interact with the observing instrument. 40 S A R Ocean Wave 5 Imaging A b s t r a c t This Chapter discusses the SAR imaging of surface gravity waves, a topic which remains controversial. The physical scattering mechanisms are outlined. The dominant scattering mechanism is Bragg resonant scattering, which is suggested to be a coherent scattering process. This is contrary to the predominant point of view. Another controversy surrounds the role and calculation of the scene coherence time. Various decorrelation time estimates are presented and discussed based upon the rms orbital velocity and a relevant cutoff length scale. Yet another controversy surrounds the roles of the various ocean wave velocity components in the SAR imaging process. It is shown that the phase velocity cannot produce a coherent modulation. Discussion of the noncoherent role of the phase velocity is deferred to Chapter 6. The coherent modulation produced by the orbital velocity leads to velocity bunching. It is shown diagrammatically that the defocussing effects of the orbital acceleration may not be negligible in a velocity bunching model. Finally, two SAR image phenomena are pointed out, observations of which may help to settle aspects of the ocean imaging controversies: first, waves propagating through ice, a case in which the high frequency wave components are suppressed; and second, narrow "V" ship wakes, which are a result of scattering from Bragg scale waves produced by the ship. The study of SAR imaging of surface gravity waves is an interesting and challenging problem which has been the subject of many recent review articles [23,47,79,80,84,85]. There are at least three areas of potential controversy: first, the (nature of the) scattering mechanism; second, the importance and calculation of the scene coherence; and third, the appropriate wave velocity for coherent Doppler modulations. Each of these controversial issues are addressed in this Chapter. 5.1 S c a t t e r i n g M e c h a n i s m s It is generally accepted that three physical microwave backscattering mechanisms could be relevant for S A R imaging of the ocean surface. These are specular reflection, wedge scattering, and Bragg resonant scattering. 5.1.1 Specular Reflection Specular reflection is conventionally understood to be direct reflection by mirrorlike planes or facets for which the angle of incidence equals the angle of reflection. 41 Two conditions must be satisfied by the scattering surface [82]: first, the facet size must be at least of the order of the radar wavelength; and second, the facet must be oriented perpendicular to the radar line of sight. Wave steepness seldom reaches a typical angle of incidence before breaking occurs. Thus, even for a satellite S A R with an incidence angle of 22°, true specular reflection will seldom contribute to the energy backscattered to the radar. However, if the wave is breaking, or becomes concave curled while breaking, this mechanism could have a role. In such a case, the microwave returns are coherently modulated by the long wave phase velocity. 5.1.2 Wedge Scattering If the ocean wave is cusped, a wedge scattering component may be important in describing the backscattered microwave energy [44]. Such a mechanism applies to regions in which the scattering surface has a radius of curvature which is small compared with the radar wavelength and is most significant when the angle of incidence is relatively large (say > 45°). The microwave returns from this mechanism are specular in nature for a given wave crest and are coherently modulated by the velocity of the cusped region which is, again, the long wave phase velocity. 5.1.3 Bragg Resonant Scattering Bragg resonant scattering is considered to be the most important backscattering mechanism for S A R ocean surface imaging. The separation wavenumber between the long waves and the Bragg scale waves has been estimated to be [23] -Kemh = ~ - > (82) where the subscript emh stands for electromagnetic-hydrodynamic. Thus, the scattering may be regarded as arising from small facets, or "patches," of nominally five Bragg wavelengths across. The complex amplitude of the reflectivity from these patches is 42 assumed to vary slowly with time and be uncorrelated with neighboring or overlapping patches [23,79,80]. The use of JC mh as a parameter of interest is admittedly heuristic. e A n enhanced surface roughness and, hence, radar cross-section, is associated with the long wave crests due to straining of the Bragg scale waves. This radar cross-section modulation mechanism is termed the hydrodynamic modulation mechanism, and is relevant for all radars operating with a Bragg wavelength of the order of 10 cm, being most important for range travelling waves [3]. The tilting of short gravity waves riding upon long gravity waves results in a spatially variable local angle of incidence, and hence, relevant Bragg subset. It has been proven that HH polarization is more sensitive to wave slopes than VV polarization, that range travelling waves are more readily imaged by this mechanism, and that the maximum radar cross-section is associated with the trailing face of a wave moving away from the radar, and with the leading face of a wave moving towards the radar, being 90° out of phase with the hydrodynamic modulation [3]. This radar cross-section modulation mechanism is termed the tilt modulation mechanism, and is relevant for steeper incidence angle radars. 5.1.4 Bragg Scattering Nature A fundamental question is whether the nature of Bragg resonant scattering is specular or diffuse. Usually, it is dismissed as purely diffuse [23]. However, a single "patch" of Bragg resonant waves produces a coherent, hence, specular-like return [64]. Furthermore, a coherent geometry is required for this mechanism to produce any microwave reflection at all. Only when an ensemble of such "patches" are present in a given resolution cell can the return become diffuse in nature. 43 5.2 O c e a n W a v e V e l o c i t i e s Attention is now focussed upon the surface motions induced by the long waves. This is the third effect of the long waves upon the Bragg scale waves. Some investigators favour models which coherently respond to the long wave phase velocity, while others favour orbital velocity models. This issue is addressed in the balance of this Section and in Chapter 6. 5.2.1 T h e P h a s e Velocity The phase velocity cannot lead to coherent Doppler modulation of the S A R signal. This assertion is based upon an argument which is adapted from Raney and Lowry [66]. Assume that the radar return from an arbitrary point target on the ocean surface is coherently modulated by the long wave phase velocity C . p The radial component of this velocity is v = C sin <b sin 7 , p (83) where the angles are defined in Fig. 8. It is a fact that any coherently sensed radial velocity component leads to Doppler suppression of the response in accordance with the first exponential term in equation (28). For matched filter processing, the Doppler suppression is given by exp{—2TT(V//?V) }. Thus, the response to a point target will 2 be essentially obliterated if f I> fi• (84) For the target to be imaged at all, the phase velocity must satisfy C sin71 s'm6\ D Y n < fi • , , (85) For typical airborne S A R parameters (V = 130 m/s, fi = 1.9°, 7 = 45°) and a 150 m ocean wavelength (C p = 15.3 m/s), the wave must satisfy 4> < 23°. That is, the wave must be essentially azimuthally travelling if it is to be imaged at all. Rangeward propagating surface gravity waves, in fact, are regularly (and predominantly) observed 44 Figure 8: Definition of the angle of propagation. 45 in airborne SAR imagery. Thus, direct coherent sensing of the long wave phase velocity cannot be a relevant mechanism in the S A R imaging process. The fashion. long wave phase velocity may enter the imaging process in a noncoherent That is, patterns of enhanced reflectivity associated with the wave crests may translate at the phase velocity. This does not necessarily imply coherent Doppler modulation at the projected phase velocity. The implications of the noncoherent role of the phase velocity are discussed in Chapter 6. 5.2.2 The Orbital Velocity Vertical velocity components imparted upon the Bragg scattering patches by the orbital velocities associated with the passage of a long ocean wave lead to S A R imaging via the velocity bunching model [1,62,77]. Positive and negative velocity components along the SAR's radial line of sight lead to forward and rearward shifts in the apparent scatterer locations. This causes a redistribution of reflectivity density in the S A R image, appearing as a linear wave-like feature under some conditions. The radial orbital velocity component is 11 v(x,t) = CIA s'm(Kx - fit) cos7 . If the wave period T w (86) — 27r/fi is much longer than the available coherent integration time T , using equation (30), a target originally located at x is relocated to x' -x 1 1 + La s i n ( i i r x ) (87) In the event that the wave is not azimuthally travelling (<f> ^ 0) the radial velocity becomes [4] v(x, t) = OA (7(7, <f>) sin (it • x + 5j , where (7(7, <j>) = \Jsin 7 sin <f> + cos 7 , 2 2 2 and 6 = tan - 1 (tan 7 sin^) . This formulation accounts for the contribution of the horizontal orbital velocity component to the net radial velocity. 46 in the resulting image, where A = 27r/JC is the wavelength, and =v ^ ^ ? • a () 88 The parameter a is the maximum fractional amplitude of the azimuthal target shift. Fig. 9 conceptually shows the effect of the redistribution of target density. The vertical lines are displaced in azimuth in accordance with equation (87). The line density in the figure is proportional to the radar reflectivity which is taken to be uniform in azimuth before the motion effects are introduced . The differential velocities tend to bunch the 12 targets near the wave trough, regardless of the relative propagation direction (i.e. to the left or the right in the figure). This radar cross-section modulation mechanism is termed the velocity bunching mechanism, and is relevant for coherent radars such as SAR. The effect is strongest for azimuthally travelling waves. Raney [62] showed that the azimuth mapping becomes critical when a = 1/27T. If a < 1/2TT, the target redistribution is one-to-one and results in a single band of enhanced reflectivity near the wave trough in the resulting image. If a > l/2n, the mapping becomes many-to-one and results in a pair of bands of enhanced reflectivity for each wavelength in the resulting image. The radial acceleration component of the wave is a(x,t) 12 = - ft A cos(Kx 2 - fit) cos7 . (89) Assuming the initial uniform target density f (x) = ( ' \ 0, J xx y Kx < otherwise otherwi i f _ , r 7T is subject to the mapping x' = x + Aa sin(Kx) = g[x) the resulting target density is fx'(x') = fx{x) 2TT dx dx' a'=ff-M*) 1 1 + 2na cos (Kx) This mapping is one-to-one if 27ra < 1 and is linear (in a linear systems theory sense) if 2na < 1 . A guideline for linear mapping is 2wa <C 0.3. If 2ira > 1, the density function has one or more singularities. Such singularities are not of consequence since they are tempered by the system resolution. 47 Figure 9: Conceptual velocity bunching diagram. The abscissa is proportional to the azimuthal coordinate. The solid line represents the nominal ocean surface. The vertical lines are proportional to the radial acceleration component; their density is proportional to the mean reflectivity; and the lines have been relocated in azimuth by the radial velocity component. Cases shown are a) a = 1 / 3 7 T (one-to-one mapping), and b) a = l/ir (many-to-one mapping). This representation including velocity and acceleration effects is new. 48 The vertical lines in Fig. 9 are proportional to this acceleration component. The defocussing effects of the acceleration component are significant if >Z aR 0 (90) for a multi-look matched filter processor. Note that large radial accelerations occur near the regions of peak target density. Thus, the defocussing effect of the radial acceleration component opposes the effects of velocity bunching and may not be negligible. This is a new result and may be important in some applications. The azimuthal orbital velocity component u is 180° out of phase with the radial acceleration component a. The relative importance of a and u may be judged by \2u/V\ _ \R a/V>\ I2KV V^JZo' 0 1 the ratio of the along track velocity defocus term to the radial acceleration defocus term. For typical S A R parameters, and all practically observable ocean waves, £ < 1 , showing that acceleration defocus always dominates along track velocity components of the orbital motion. Appendix D discusses a one-dimensional simulation of monochromatic, azimuthally travelling ocean waves which has been utilized to further illustrate velocity bunching. It is shown that acceleration defocus may be important and will oppose the effects of velocity bunching. It is also shown that Doppler suppression may be important for an airborne S A R . 5.3 Scene Coherence We now turn to estimation of the relevant scene coherence time in the presence of ocean waves. In Chapter 2, it was shown that target motion may be modelled as a multiplicative phase given by B /3{t) = exp j - j [ - (t t D 49 + rjt ) - 0(r)] } , 2 (92) where the motion parameters tp = RQV/V and fj = R a/V 2 2 0 are the results of averaging over the duration of a given look. (The along track orbital velocity component has been neglected as per Section 5.2.2.) The random variable 6(t) represents the coherence of individual patches of Bragg waves. Of relevance here is the intrinsic coherence time of such a patch, and the decorrelation time of such a patch from other such patches. 5.3.1 Intrinsic Scale The intrinsic coherence time scale refers to the time over which a given patch of Bragg waves coherently retains it structure. There is no theory nor explicit estimates for such a time scale. However, a general assumption is that the complex reflectivity amplitude from such a patch varies slowly [23,79] so that the intrinsic coherence time scale is long compared with the available coherent integration time. 5.3.2 D e c o r r e l a t i o n Scale The decorrelation coherence time scale refers to the time over which a given patch of Bragg waves remains co-phase relative to itself, and relative to other such patches. This time scale is generally related to the rms surface velocity, v , as indicated in Tms equation (39). Thus, estimation of r requires estimation of v c Tms over the ensemble of wave motions. The magnitude of u r m s may be calculated from the integration [2,46,79,80] « £ n . = r n s {n)dn, 2 2 ( (93) where S (Q) is the one-dimensional frequency spectrum of surface elevation. A closed ( form expression may be found for u 2 m a by using the Pierson-Moskowitz spectrum. Then, we find that <. = /f(^)-e {^(4) }cosS[;. 2 rf It remains to specify the limits of integration. 50 M The high frequency limit Cl has not been controversial. This limit corresponds 2 roughly to the frequency of the Bragg waves and one can use nt — • oo with little 2 error. Then, we find that • "rms - ' ' ( ^ fix 4) (* ( fi a U 2 V \ 4 _r f (95) \unn I /^? • erf A 1 1 where the deep water dispersion relation is used to replace Jli by A i . If the argument of the error function is small compared with unity, we may write = ^A cos 7, (97) e— \ ~ r • y ag A i cos 7 (98) «r 2 2 m a 1 so that T The wavelength A represents a sharp cutoff between waves which are treated in x a statistical manner (A < A i ) and waves which are treated in a deterministic manner (A > A i ) . The interesting issue is reduced to consideration of the lower limit A . There x are several choices. Some investigators [2,79,80] choose A i = 2p: twice the matched filter processor, one-look resolution. Then T = 4.45— . sJP cos 7 C (99) Typical coherence time scales based upon this and other results derived in this section are summarized in Table V for the S E A S A T S A R and the C-SAR. Such a choice for the cutoff wavelength A i is questionable for several reasons. First, no account is given for the effect (if any) of multi-look processing. Second, the resolution scale is dependent upon the processor configuration (and for that matter, the resolution measurement norm) and hence, is arbitrary to within a factor of (say) two. Finally, the actual relevant one-look resolution may be substantially degraded through the presence of target radial accelerations, and may be orders of magnitude larger than the intended 51 Table V : The coherence time scale in msec for the various cutoff wavelengths considered. SEASAT C-SAR Ai 2p 445 315 J_rp2 550 355 oo 140 40 {H = 4 m) s static theoretical one-look resolution p. Given the variability of the effective resolution, this approach to estimating the cutoff wavelength scale might best be done through iteration. A second choice for the cutoff wavelength A i is based upon the available coherent integration time. Any wave with period less than or equal to T will pass through at least one complete cycle over the course of image formation. The corresponding cutoff scale is Hi = ^ T 2 , (100) which results in T C = 5.04 A T cos 7 . (101) A third choice for the cutoff wavelength Ai is to assume that the velocities of all spatial scales contribute to v ims (as a worst case estimate) [46]. Then A i — • oo so that (TT (aU \ 2 V "rms = 4 / \l p o cos 7 . (102) Thus, we find that r = 1.12c 'H cos 7 (103) s is a lower bound estimate of the coherence time. Results for this formulation based upon a 4 m significant waveheight are included in Table V . A n order of magnitude range of coherence time scales is available from these approaches as is illustrated in Table V . The estimates based upon inherent S A R length 52 scales seem independent of both sea-state and ambient environmental conditions—a situation which seems unlikely. The estimate based upon the entire wave spectrum seems more sensible from this point of view in that there is significant waveheight and hence, windspeed dependence. However, as we will see in Chapter 8, even the shortest estimated coherence time scale based upon the latter estimate is rather longer than that which is observed. This indicates the inadequacy of the decorrelation time scale approach, and might suggest that the intrinsic scale is somewhat shorter than indicated. 5.4 W a v e s i nI c e The first observations of ocean waves propagating into an ice-covered region in SAR imagery were made about 10 years ago by Dawe and Parashar [14]. Recently the analytical opportunities presented by such observations have been suggested by Raney [62]. A n important observation is that the intensity patterns of such waves appear cusped in S A R imagery when an azimuthal wavenumber component is present. This suggests that focus and coherence time limitations are not issues. Focus is shown to be negligible in Appendix D for multi-look SAR's in a one-to-one velocity bunching model. Furthermore, the ice-cover is a low-pass filter, effectively eliminating the high frequency surface waves which decrease the coherence time. Thus, the same gross dynamics are present in the ice-covered case as in the open ocean case, but there is no coherence time limitation. Recent L E W E X observations [15] (see Appendix C) successfully obtained S A R imagery of ocean waves with an azimuthal wavenumber component propagating into an ice covered region. Representative imagery of waves in ice is shown, discussed, and compared with simulations in Chapter 8. 5.5 N a r r o w The " V S h i p W a k e s S A R image of a ship wake seldom appears to be dominated by the classic Kelvin wake pattern. More usual, if the ambient conditions are suitable, is a narrow 53 "V" appearance, which is generally associated with azimuthally travelling ships. It has recently been established [41,42] that this S A R image phenomenon is a result of scattering from Bragg scale waves excited by the passage of the ship. This result is based upon the observation that the narrow "V" wake opening half-angle is consistent with the angle (104) where C g is the group velocity of the Bragg scale waves, <pi is the ship heading with respect to azimuth, and V is the ship speed. s Such imagery provides an opportunity to directly observe and assess the nature of Bragg scattering since the linear target pattern is a direct result of Bragg scattering. The coherence time associated with the point targets which constitute the wake is long because such image phenomena persist over large spatial scales. The imagery of a linear target (a line of point targets), a model for one arm of a narrow "V" wake, is discussed in Appendix B. It is shown that the width of the linear target modeled in such a fashion will respond to a focus error in the perpendicular image dimension—the range dimension for a narrow "V" ship wake. These points are discussed further in Chapter 8. 5.6 S u m m a r y This Chapter has discussed aspects of the S A R imaging of ocean surface waves. For example, it is asserted that the Bragg scattering associated with the ocean surface is a coherent, specular-like process. This is not the traditional interpretation of Bragg scattering, although it is consistent under logical analysis. Coherence time scales based upon the decorrelation of Bragg scattering cells were calculated for various cutoff wavelength scales and a Pierson-Moskowitz wave spectrum. The scales calculated are summarized in Table V and range from 40 msec to 500 msec. It was shown that the wave phase velocity cannot lead to coherent Doppler modu- 54 lations. The wave phase velocity enters the imaging process as the translation velocity of the mean reflectivity envelope. This concept is discussed further in Chapter 6 . Only the orbital dynamics may lead to coherent Doppler modulations. The radial component of the orbital velocity leads to velocity bunching while the radial component of the orbital acceleration may cause image defocus. This Chapter concludes by introducing two S A R image cases which may help to solve some of the controversy which surrounds S A R ocean imaging. These are waves in ice, a case in which coherence time limitations are not an issue, and the narrow "V" ship wake, a case in which Bragg scale waves are directly observed. 55 6 D y n a m i c Reflectivity Patterns A b s t r a c t In this Chapter, the effects of the dynamics of the mean pattern envelope of the ocean's reflectivity are distinguished from the effects of the dynamics of the individual point targets which constitute the reflectivity pattern. It is argued that the effect of frequency offset look-extraction windows combined with envelope motion causes the envelope image to become misregistered when combining looks. It is shown that the resolution, and hence the image contrast, is minimized when the fractional processor focus error is adjusted by n = VE/2COS^>. This result is typically interpreted as arising due to a coherent Doppler modulation at an azimuthal velocity of V E / 2 COS (f>. In fact it arises because one effect of a processor focus error is to adjust the registration of the individual looks. Thus, the processor focus adjustment is simply re-registering the individual looks which were misregistered by the noncoherent envelope motion. This re-registration is at the expense of actual misfocus of the processor. True optimal processing uses fixed correct focus and independent re-registration of the looks. These observations constitute the contributions of this Chapter. The case of an airborne multi-look S A R imaging an azimuthally travelling reflectivity pattern is considered in this Chapter. The objective is to clarify issues surrounding the controversy over the effect of a processor focus adjustment upon the image contrast of such scenes. Observations indicate that a processor focus perturbation proportional to the projected pattern velocity will improve the image contrast. We will show that this does not imply a coherent Doppler modulation at this velocity. The consequences have an important bearing upon the scattering mechanisms being investigated. 6.1 D e v e l o p m e n t Assume that a S A R is being used to observe a moving reflectivity pattern. By way of example, consider the imaging of a wind gust pattern or "cat's paw" excited upon an otherwise calm water surface. The wind roughened region may be modelled as consisting of many small coherent regions or patches of roughness, each of which produces a coherent return to the radar such as that derived from a point target. The mechanism for producing such returns is irrelevant to the present discussion. The ensemble of all such patches, or point targets, constitutes the noncoherent reflectivity pattern that we associate with the radar image. 56 There are two sets of dynamics which are relevant to this example. the dynamics of the individual point targets. The first is The motions of the individual targets can coherently (phase) modulate the returned radar signal. These effects (discussed in Section 2.2 for point targets and in Section 3.1.2 for diffuse targets) are the ones most frequently found in the S A R ocean literature. The second set of dynamics is the motion of the pattern itself, particularly the velocity of the envelope of the mean reflectivity pattern. The reflectivity envelope velocity cannot coherently modulate the returned radar signal. The reflectivity envelope velocity is not necessarily the same as the velocity of the individual point targets which constitute the pattern. Differences between these velocities are vital since the targets which constitute the pattern do not move at the wave phase speed (Chapter 5). A recognized effect of the reflectivity envelope is the scanning distortion which arises when the envelope velocity becomes an appreciable fraction of the platform velocity [66]. This effect is regularly observed in airborne S A R imagery of ocean waves and wave imagery derived from other airborne scanning sensors such as the surface contour radar [87]. We will now consider another role of the azimuthal velocity component of the reflectivity envelope VE, in the event of multi-look processing. Recall that in a multi-look processor, data windows are applied in the frequency domain to generate the individual looks. These data windows may be separated by a Doppler frequency of Af . From equation (10), this corresponds to a time separation between the looks of At = A fo/K . D a This implies that the image of the reflectivity envelope will be misregistered in the azimuth between the two looks by a distance Ax = vsAt. Then, in the output image plane, the corresponding temporal offset is 6t = ( U E / V ) A £ . Thus, the effect of the envelope velocity is to differentially shift the relative location of each look in azimuth. This same effect has been demonstrated in two-scale radar wave probe data [56]. The analysis of Chapter 2 may be extended to include the effect of the reflectivity 57 envelope velocity. The net effect is to shift each look in azimuth by an amount 2 (TV U_V_E 6t where (T/Ti) ti/K KV \Tj (105) ' is the centre of the effective look extraction window. This effect 2 may be introduced into the formulation after coherent processing is complete since it is a noncoherent effect. The single look dynamic impulse response is modified by: 9d{t,U) Qd (-( r\2 t c_ Tj K V y (106) ± M t Then, after look-summation we find that rp2 9dm{t) , = S • exp jKL + e 2 + i(f Y[el-2eeo} ^ l BK 2 exp t+ -27T- K 2 + \ ( g ) + Ke /L 2 {-27T^-\ 1- 6rj [e - 2ee ] 2 1 (T\ + LK\fJ 2 c v tD 0 (107) where e = TB(VE/V). 0 Pm This result has resolution = (108) \ It can be shown that the resolution p is minimum when m 1_ fT_\ K \T,) 2 VE _ V ~ N 2 N 2 - l V ' (109) where the second equality applies to both matched filter and infinite rectangular filter processing. For N > 1, we find that n ~ VE/V will minimize the resolution, and hence, maximize the image contrast . This result arises because one effect of a focus 13 adjustment is to change the azimuthal registration of the looks. Thus, a change in the coherent processor parameters decreases the impact of this noncoherent effect upon the output image. 13 See footnote 6 on page 17. 58 A better technique to provide such a correction is to noncoherently combine the sequence of coherently generated looks, shifting each look in accordance with the offset produced by the envelope velocity given in equation (105). Using this approach, the image contrast will be further enhanced because no broadening of the individual impulse responses will result, as is the case when the focus parameter is adjusted to align the misregistered looks. This degradation in contrast due to the envelope velocity also applies to the case of one-look processing. This is because a moving scene is being imaged over time by a moving sensor (analogous to a time-exposed photograph of a moving subject). The output is smeared in accordance with the velocity ratio VE/V over the course of the observation time. This effect applies within each look of a multi-look set as well. However, as N becomes large, the observation time is reduced, as is the degree of smearing. Thus, as the number of looks becomes large, the above analysis becomes more exact. A n identical focus perturbation will cause the S A R to track the moving reflectivity envelope and similarly improve the contrast of a one-look image. This also applies to the individual looks of a multi-look image [67]. If the processor focus error required to register the looks is related to an equivalent azimuthal velocity component v, we find that v = v /2. E It is tempting to interpret the effect of the envelope velocity as a coherent Doppler frequency modulation leading to a broadening of the response due to a focus error. However, to make such a leap is incorrect. It illustrates the danger of interpreting processor focus perturbations in terms of an equivalent coherently-sensed velocity component, which in turn implies oceanic targets moving with velocity vg/2. Of course, such is not the case. This issue is one of intense controversy in this area due to the experimental focus paradox which this section addresses. The effect of the velocity of an envelope of reflectivity applies to other notable 59 situations. One example is the imaging of ocean waves with an azimuthal wavenumber component. Through several mechanisms an imaged reflectivity pattern is produced which travels with the wave at the wave phase velocity C . p However, the individual point targets which constitute the pattern are coherently modulated by the wave's orbital velocity. Thus, processor focus perturbations proportional to the projected azimuthal component of the phase velocity will tend to register the looks and improve the contrast in both the wave image and wave spectra derived from the image. If the ocean wavenumber vector makes an angle <p with the platform velocity vector, the projected azimuthal component of the pattern velocity is C j cos <j). The processor focus v perturbation required to "optimally" focus S A R imagery of the waves is n = C /2 cos <f>. p This result is consistent with observations [28,72,74] that ocean wave image contrast is maximized by a focus adjustment proportional to the projected azimuthal phase velocity. (This effect is not observed in satellite S A R imagery [73] simply because C p / V is extremely small and the depth of focus is relatively large.) It is stressed that this result does not imply a coherent Doppler frequency modulation proportional to the projected phase velocity, and that true optimal processing uses fixed correct focus, and independent re-registration of the looks. 6.2 S u m m a r y This Chapter has addressed the apparent paradox that optimal S A R ocean wave image contrast is observed when the processor focus is perturbed in proportion to onehalf the projected phase velocity, yet the phase velocity cannot produce coherently sensible Doppler modulations. The paradox is resolved by noting that the mean reflectivity envelope translation velocity is the wave phase velocity, and that this translating pattern is imaged over a non-zero period of time by the SAR. In the event of multi-look processing, the envelope translation leads to misregistration of looks which may be corrected by the previously mentioned focus perturbation. This is because one effect of a 60 processor focus perturbation is to change the relative registration of the looks. True optimal processing should retain the correct focus and independently re-register the looks. 61 7 Analysis Methodology A b s t r a c t This Chapter presents the analysis methodology used in the remainder of this work. A new two-dimensional simulation model is motivated, with development contained in Appendix A. Outputs of the simulation may be compared directly with actual SAR image phenomena. This Chapter also contributes a methodology to extract three statistical parameters from multi-look SAR data which may be used to address the question of target density in an arbitrary SAR image. The parameters include the inter-look crosscorrelation, the peak-height scaling in response to look summation, and the azimuthal image bandwidth. These parameters may be extracted from simulated or actual SAR imagery in the same manner. The analysis tools used in this work fall into two categories: first, a two-dimensional S A R image simulation model is developed to produce controlled images that may be compared qualitatively with actual SAR imagery; and second, a series of statistical tests is developed specifically to address the nature of scattering in a S A R image. These tests may be applied to simulated imagery to quantify scene behaviour for known parameters. These tests also may be applied to actual S A R imagery, and the results compared with those from simulations to better relate S A R imagery to known parameter variations. 7.1 S i m u l a t i o n As an aid in understanding the processes involved in S A R ocean surface imaging, a two-dimensional simulation package has been developed. The simulation procedure takes careful account of the scattering nature; there are minimal a priori assumptions as to the statistical behaviour of each resolution element, and each point target in the scene is explicitly considered. Thus, the behaviour of the targets in each resolution element is explicitly known. This degree of control and flexibility is not available when studying real S A R imagery. Such a method presents the opportunity to produce a database for comparison to actual S A R imagery, and represents a distinct departure from other simulation based studies. SAR image simulation has a strong precedent in the design and performance evaluation of various aspects of S A R response and sensitivity. One thrust of such research 62 has been end-to-end S A R product simulation [21,24], intended primarily to assess geoscience applications, and to aid in S A R system design trade-offs. Such simulations generally require as input some sort of "radar reflectivity map," and in some manner generate a speckle pattern under the assumption of many independent targets, thus implying Rayleigh statistics and derivative forms for the multi-look case [68]. Although simulation should play a role in understanding S A R ocean surface imaging [72], the scene dynamics need to be appropriately represented, as they provide both the interest and difficulty essential to the problem. A one-dimensional (azimuthal) model has been presented by Alpers [l] in which velocity bunching allows image intensity modulations to arise from a randomly generated wave field. More recently, this type of work has been extended to two-dimensions [2,43,86]. These works generally contain one or several inherent assumptions which could prove to be restrictive: diffuse scattering, a large number of looks, or a long scene coherence time. Obviously, implementation of these assumptions into a simulation limits the utility of the simulation for testing the assumptions themselves. The details of the simulation methodology developed for this work are contained in Appendix A . The simulation is specifically designed to properly represent dynamic distributed targets. It has provisions to simulate varying number of targets per resolution cell, azimuthal wave imaging through the velocity bunching mechanism, linear or distributed targets, target acceleration and processor-induced focus errors, and finite scene coherence time. The simulation is based upon fundamental S A R imaging theory and gravity wave behaviour as developed in previous chapters. 7.2 S t a t i s t i c a l T e s t i n g M e t h o d o l o g y Table IV on page 33 outlines salient image features which may be used to differentiate between examples of SAR imagery derived from point and diffuse targets. In this section, these features are developed into a statistical testing methodology which may 63 be invoked to test for the presence of point targets in an arbitrary S A R image. 7.2.1 Inter-Look Cross-Correlation Observation of the correlation between looks should technically occur on a resolution cell-by-resolution cell basis. However, such a procedure is impractical for two reasons: first, a focus error shifts the relative azimuthal positions of the resolution cells from look to look, making it impossible to track the appropriate resolution cells for a correlation measure in the event of a dynamic ocean surface; and second, a low level of statistical significance is associated with resulting correlation estimates due to the small number of statistically independent samples (essentially equal to the number of looks N). These factors suggest that a spatial correlation measure would be more profitably employed. The adopted procedure is a straightforward calculation of the correlation coefficient between pairs of looks, applied at nominally perfect focus only. Such a measure is statistically rather reliable because the number of degrees of freedom in the estimate is twice the number of statistically independent resolution cells in each look. About 1000 degrees of freedom are expected for a four-look 128-by-128 pixel subscene from the S E A S A T S A R . In the event of target motion or look extraction filter overlap, the number of independent resolution cells is decreased. In this work, the derived statistic is the mean of the correlation from all possible pairs of looks. In the event of an isolated point target, this statistic is rather dependent upon the processor focus due to possible misregistration of looks in azimuth. Observation of the peak of a cross-correlation function may be more useful in such cases. However, in this work the correlation is measured only at nominally perfect focus. 7.2.2 P e a k - H e i g h t Scaling To observe the peak-height scaling as a function of the number of looks to which the S A R data is processed requires that the data be reprocessed with a variety of choices of 64 N. However, it is impractical and unnecessary that this level of processing be carried out when the individual looks of an N look set are available. In this case, an M look set may be simulated by summing M of the N available looks (where M < N) followed by a peak-height measurement. Using such a method, the results differ slightly from those indicated in Table IV, but there is no change in interpretation, as discussed presently. As a measure of peak height, we choose the standard deviation normalized by the mean for the M chosen looks. For either a point or diffuse target, the mean value is proportional to M. In the case of a diffuse target, the standard deviation is proportional to M / , and the peak-height measure will scale as M 1 2 - 1 / . However, for a point target, 2 the peak value will be proportional M so that the peak-height measure will scale as M ° , in contrast to N' 1 as expected for iV-look processing. This measurement is completed by calculating the mean and standard deviation of all possible combinations of M looks of the available N looks, averaging the results for the peak-height scaling, and performing a first-order least-squares fit to the logarithm of the peak-height scaling data. The coefficient of the linear term is the relevant measure. Again, this measure is dubious in the event of an isolated point target and a focus error. Point targets may become misregistered, and their sums will not scale as indicated if the misregistration is rather large. This is because the looks considered in this work are derived from a discrete frequency-plane mixed-integrator. Thus, in practice, this measure was performed only upon nominally perfectly focussed imagery. 7.2.3 Azimuthal Image Bandwidth For a point target, a processor focus error will broaden the image response, thus leading to a decreased image bandwidth. However, for a diffuse target, a processor focus error will cause no change in the second order statistical behaviour of the speckle pattern, and the image bandwidth will be unaffected. The azimuthal image bandwidth is considered to be the most important measured 65 parameter because it is explicitly dependent upon processor focus in a well defined fashion for the point and diffuse target cases. The inter-look cross-correlation and the peak-height scaling measures may be regarded as supporting the conclusions derived from this measure. The azimuthal image spectrum is calculated for each look by removing the mean, applying a data window and an F F T to each azimuthal line, and averaging in range to produce an azimuthal spectrum with a number of degrees of freedom equal to twice the number of statistically independent resolution cells in range. The bandwidth is measured on the basis of an equivalent rectangular width for a Gaussian shaped spectrum. This procedure is applied on a look-by-look basis rather than upon the multi-look image for several reasons: first, the variation in azimuthal bandwidth for a point target is adequate on a look-by-look basis for moderate focus errors in a S E A S A T SAR scene to indicate the presence of a point target; second, if the looks are sufficiently misregistered due to the processor focus error, the peaks in the individual looks will appear in the multi-look image, thus corrupting the spectral shape and corresponding bandwidth measure; and third, the bandwidth derived from each of the N looks may be averaged to provide a more statistically reliable measure of the azimuthal image bandwidth. 7.3 S u m m a r y This Chapter has presented the analysis tools which will be used in the remainder of this work. The two-dimensional simulation is unique in that it explicitly considers each of the point targets which are present in a given scene. Details of the simulation are found in Appendix A . The three image statistics which may be derived from multi-look S A R imagery are the inter-look cross-correlation, the peak-height scaling as looks are summed, and the azimuthal image bandwidth as a function of the azimuthal processor focus error. These measures are relatively easy to extract and subject to good statistical reliability. As well, they help to differentiate between S A R imagery of point and diffuse 66 targets. 67 8 S i m u l a t e d v s . A c t u a l S A R ImageryA b s t r a c t In this Chapter, we perform comparisons between simulated and actual SAR imagery. The comparisons are either quantitative on the basis of extracted observable statistics, or qualitative on the basis of one-to-one comparisons of actual SAR imagery with simulation products. Four sets of comparisons are performed to investigate various SAR ocean imaging characteristics. The first set of comparisons examines the specular versus diffuse nature of ocean surface scattering. It is shown that test statistics of the simulation and a SEAS AT SAR test scene demonstrate the appropriate statistical behaviour for point or diffuse scattering. When applied to an ocean scene the statistics suggest diffuse scattering only. This result is explained on the basis of random target motions and finite scene coherence times leading to large resolution cell sizes. The second set compares simulations of linear targets with SEASAT SAR imagery of narrow "V" ship wakes. The linear targets and narrow "V" wakes should behave in a specular manner. The "V" wake represents direct observation of ship generated Bragg waves, suggesting that Bragg scattering is a coherent process. The third set compares simulations of azimuthally travelling ocean waves having various mean scene coherence times with C-SAR imagery of waves in the open ocean propagating into an ice-covered region. The role and degree of coherence time limitation in airborne SAR imagery is demonstrated. A scene coherence time of 14 msec seems appropriate for the ocean case. The fourth set considers near and far range C-SAR imagery of waves propagating through an ice-covered region. It is shown that the range dependent azimuthal look misregistration due to the translating reflectivity envelope is consistent with theory. This reinforces the noncoherent motion effects discussed in Chapter 6. This Chapter draws quantitative and qualitative comparisons between simulated and actual S A R imagery. Comparisons are drawn in four areas: first, to assess the ocean surface target density in S E A S A T S A R imagery; second, to assess the nature of the scattering associated with narrow "V" ship wakes in S E A S A T S A R imagery; third, to assess the size of the coherence time in C - S A R wave imagery and wave in ice imagery; and fourth, to directly observe the effects of a moving envelope of reflectivity. 8.1 S p e c u l a r v s . D i f f u s e S c a t t e r i n g Recall that microwave reflectivity from the ocean surface is largely a Bragg scattering process. Also, recall the assertion that Bragg scattering is a coherent scattering process for each Bragg scattering region. We will now examine this assertion by applying the statistical tests developed in Chapter 7 to three data sets: first, to simulated 68 S E A S A T S A R imagery of variable target density and processor focus errors; second, to actual S A R data of a scene containing known examples of point and diffuse targets obtained by S E A S A T [33]; and third, to an actual S E A S A T S A R scene of the ocean surface. 8.1.1 S E A S A T S A R Simulation S E A S A T SAR scenes of variable target densities with various processor focus errors were simulated. The test statistics were extracted from the results. Details of this analysis are contained in Appendix B . In that Appendix, it is demonstrated that the simulations exhibit appropriate asymptotic behaviour for approach to the diffuse and specular scattering limits on the basis of the extracted test statistics. However, these asymptotes are approached rather slowly. It is also apparent that even cases which behave in a rather specular-like manner are not accompanied by noticeably different image appearance, perhaps suggesting why ocean reflectivity is traditionally dismissed as being diffuse in nature. Based upon these observations, subsequent simulations of ocean surface phenomena were performed with a rather small target density—generally 3 targets per resolution cell. This choice is less than the traditional guideline of 5 [34], yet does not noticeably impact qualitative comparison of actual S A R imagery with the simulation product, as shown presently. 8.1.2 S E A S A T S A R Control Scene The control scene chosen is Goldstone California, taken from S E A S A T orbit 882. This scene contains an array of immobile corner reflectors which are used as examples of discrete point targets, and regions of desert which are used as examples of diffuse targets. The scene is shown in Fig. 10, and the extracted 128-by-128 pixel subscenes used in the testing are indicated. Subscene 1 contains the corner reflector array which 69 Figure 10: S E A S A T S A R Goldstone California scene showing extracted subscenes. The azimuthal coordinate is horizontal while the (slant) range is vertical. 70 Figure 11: Goldstone subscene 1 with the five processor focus errors, a) 0%, b) 1%, c) 2%, d) 3%, and e) 5% corresponding to about 0, 1, 2, 3, and 5 depth of foci error in each look of this four look set. The azimuthal coordinate is horizontal. correspond to isolated point targets. Subscenes 2, 3, 6, and 7 contain expanses of desert which correspond to diffuse targets. Subscenes 5 and 9 contain imagery of mountains which are regions of high contrast. Subscenes 4, 8, and 10 contain regions with lower contrast. Fig. 11 shows subscene 1, the case of discrete point targets, with the five processor focus errors. As the focus error increases, the peaks associated with each corner reflector broaden in azimuth. Eventually, the individual peaks associated with the response of each look become apparent in the multi-look image. This illustrates why the peakheight scaling and cross-correlation measures do not apply to isolated point targets with a large processor focus error. (The look misregistration also illustrates the autofocus algorithm.) Fig. 12 corresponds to subscene 2, a diffuse case. As the processor focus error changes, there is very little change apparent in the speckle or the scene contrast. Note that the range of focus change is the same as that used in Fig. 11. 71 Figure 12: Goldstone subscene 2 with the five processor focus errors, a) 0%. b) 1%, c) 2%, d) 3%, and e) 5%. The azimuthal coordinate is horizontal. Fig. 13 corresponds to subscene 5, a high contrast mountain region. The image contrast is smoothed in this high contrast region as the processor focus error increases. The result of the statistical testing methodology are shown in Fig. 14 and in Table VI. The case of a point target (represented by subscene 1) exhibits a high inter-look correlation of 0.714 with a peak height scaling of M~° , and the image bandwidth is responding strongly, and as predicted, to the processor focus perturbations. The case of diffuse targets (represented by subscenes 2, 3, 6, and 7) exhibits an average inter-look correlation of 0.060 with a peak height scaling of M - 0 , 4 3 6 , and the image bandwidth is insensitive to the processor focus perturbations. The remaining cases of high and low contrast imagery exhibit intermediate results for each of the extracted statistics. The image bandwidth measure for the high contrast mountain scenes represented by subscenes 5 and 9 have a relatively low and constant image bandwidth for all processor focus settings. This is because the structure along the linear image features is being measured, rather than the bandwidth of the image speckle. Such effects are illustrated 72 Figure 13: Goldstone subscene 5 with the five processor focus errors, a) 0%. b) 1%, c) 2%, d) 3%, and e) 5%. The azimuthal coordinate is horizontal. Table VI: Peak-height scaling and inter-look cross-correlation for the Goldstone subscenes. subscene target type peak-height 1 2 point diffuse 3 4 diffuse low contrast -0.081 -0.464 -0.406 -0.366 0.714 0.028 0.086 5 high contrast 0.567 6 7 diffuse diffuse -0.139 -0.467 -0.405 8 9 low contrast high contrast -0.349 -0.134 10 low contrast -0.382 scaling 73 crosscorrelation 0.146 0.027 0.099 0.166 0.588 0.122 c o - CD • O | iiI 0.00 M I I li| I 0.01 M M I I ii|' ' 1 1 1 1 0.02 ' M 1 1 111111 0.03 M ' i M ii 0.04 I I I 0.05 Focus Error Figure 14: Image bandwidth measure for the Goldstone subscenes. 74 and discussed further in Appendix D in connection with imagery of linear targets and narrow "V" ship wakes. We see that nominally point target and diffuse target cases approach their predicted behaviour for each of the extracted statistics. Note, however, that even those cases chosen as being ideal are not behaving exactly in an ideal fashion. Having established the behaviour of S E A S A T S A R / G S A R imagery to the test statistics, we will now extract these same statistics from an ocean test scene. 8.1.3 S E A S A T S A R Test Scene The test scene chosen is a region of ocean in the vicinity of Duck Island, North Carolina, taken from S E A S A T orbit 1339. This data set, referred to as Duck-X, has been the subject of intense study [6,7]. Duck-X is known to contain a 200 m wave system propagating at about 45° to azimuth. The scene and extracted 128-by-128 pixel subscenes are shown in Fig. 15. The ten areas chosen for testing were extracted at random. The dominant wave system is apparent in the image. Fig. 16 shows subscene 1 with the five processor focus errors. As the processor focus error increases, the image contrast due to the wave system is lost, as discussed further in Appendix E . Results of the statistical testing methodology are shown in Fig. 17 and in Table VII. It is apparent that each of these regions is exhibiting the behaviour expected for nearly purely diffuse targets. This result is disappointing, but not surprising for several reasons. First, since the spatial resolution of S E A S A T is rather large, there is an ensemble of many point targets present in each subscene, and each target is subjected to the orbital dynamics of the underlying wave train. Thus, each target has different motion parameters and the observed result is an average over the ensemble of targets. In the case of the simulation and the control scene, no dynamics were involved and the focal 75 Figure 15: S E A S A T S A R Duck-X scene showing extracted subscenes. The azimuthal coordinate is horizontal while the (slant) range is vertical. 76 Figure 16: Duck-X subscene 1 with the five processor focus errors, a) 0%, b) 1%, c) 2%, d) 3%, and e) 5%. The azimuthal coordinate is horizontal. Table VII: Peak-height scaling and inter-look cross-correlation for the Duck-X subscenes. crosssubscene peak-height correlation scaling 1 5 6 7 -0.451 -0.436 -0.445 -0.444 -0.430 -0.439 -0.435 0.034 0.051 0.045 0.044 0.057 0.048 8 -0.433 9 10 -0.434 -0.392 0.059 0.053 0.108 2 3 4 77 0.053 O I i i i i I Ii 0.00 I I Ii 0.01 I Focus Error Figure 17: Image bandwidth measure for the Duck-X subscenes. 78 properties of the targets were systematic and identical. Second, but perhaps of less importance, the targets in the ocean scene could be subject to coherence time limitations, and hence be rather broad. This would soften any predominant peaks due to underlying specular processes. This makes the scene appear to be more diffuse in nature, even if the target density were rather low. The principle difficulty with this experiment has been the rather poor signal-tonoise ratio. A n individual target cannot be tracked as the focus is perturbed for the test scene, in contrast to the control Goldstone corner reflector array. In order to fully test the coherent Bragg cell hypothesis, we would like an ensemble of Bragg targets with a relatively long coherence time and a high signal-to-noise ratio. The image of a narrow "V" ship wake provides such an opportunity. 8.2 L i n e a r T a r g e t s a n d N a r r o w " V S h i p W a k e s As discussed previously, the narrow "V" ship wake presents an opportunity to explicitly observe an ensemble of Bragg scatterers. In this section, we discuss simulations of linear arrays of point targets, the resulting statistical measures, and draw qualitative comparisons with actual S A R imagery of narrow "V" ship wakes. 8.2.1 SEASAT SAR Simulations S E A S A T S A R scenes were simulated with various types of linear target features and with various processor focus errors. The test statistics were extracted from the results. Details of the analysis are contained in Appendix B . It is demonstrated that a large relative cross-section is required to differentiate the linear target from a nominally diffuse background, and that the best sensitivity to the processor focus test is achieved by a linear target oriented perpendicular to the S A R channel in which the processor focus is perturbed. This is the range channel for narrow "V" ship wakes—a S A R image phenomena which seems exclusively associated with near azimuthally travelling ships. 79 Figure 18: Examples of ship wakes in S E A S A T S A R imagery. 8.2.2 SEASAT S A R Ship W a k e Imagery Fig. 18 shows S E A S A T S A R imagery of several examples of ship wakes taken from a scene in the Irish Sea. Two types of wakes are apparent in this figure: first, the fishing fleet in the top image produces wakes with relatively large opening angles; and second, the ships in the middle and lower images produce narrow "V" wakes. The lines associated with these wakes are very similar to the linear targets simulated in Appendix B . Unfortunately, this imagery was not available in a form suitable for the 80 prescribed statistical testing. We may only assert that in the event of a change in the range channel processor focus error, the range impulse response will be broadened and the wake image will be smeared as was shown to be the case with the simulated linear target features. This, in itself, would suggest that the scattering from these Bragg scale waves is coherent in nature and that the proposed model for Bragg scattering is suitable. The fact that ship-generated Bragg scale waves, scatter microwaves in a coherent manner suggests that Bragg scattering, in general, is a coherent process. This may be a difficult assertion to defend because typical ocean surface Bragg scale structure is not as well organized as the waves produced by the passage of a ship. However, all ocean surface roughness or waves of the Bragg scale must obey the same dispersion relation, and so, must be of the same inherent nature within the narrow band observed by the radar. Differences between the general ocean roughness case and the ship wake case lie with the associated coherence time scale rather than in a different scattering process with a different inherent nature. This was demonstrated with the Duck-X ocean test scene. 8.3 S c e n e C o h e r e n c e T i m e We now turn to a very different dynamic phenomenon: scenes of ocean waves propagating from the open ocean into a region covered by floating ice floes as observed by an airborne S A R . The actual S A R data was obtained during the L E W E X experiment [15] by the C C R S C - S A R digital S A R [40]. The waves in ice imagery presents the same dynamics as in the open water case, but the coherence time has been effectively increased through the suppression of the high frequency wave components. This subtle but important point deserves emphasis: with floating ice, the targets are, in effect, permanently present—the issues of Bragg scattering and coherence time are removed. We are left with a pure experimental representation of a velocity bunching model. 81 Two actual S A R images were chosen for this comparison: first, a scene showing ocean waves propagating near and towards an ice edge; and second, a scene showing the waves propagating through the ice. These two scenes were obtained within 30 seconds of each other. The ice and ocean images are shown on the upper-left and lower-left of Fig. 19 respectively. The corresponding two-dimensional Fourier transforms are shown in similar positions in Fig. 20. The image transforms have been smoothed using a Gaussian kernel, and circles of constant wavelength are indicated. The image transforms of Fig. 20 indicate a rather interesting situation. We see that the same dominant wave component is present (labeled A l ) in both the ice and ocean cases (i.e. A « 150 m, 4> « 45°). In the ice case, the second harmonic is also apparent (A2), indicating that the wave image is rather cusped in this case. In the ocean case, the cusped nature of the response is smoothed by the effects of the finite scene coherence time, and the second harmonic is not readily apparent. Furthermore, it is apparent that in the ice case there is an additional wave component present which corresponds to a shorter, more azimuthally travelling wave (Bl). There is no evidence for the corresponding wave in the ocean case. Thus, based upon the ice image, the ocean has a bimodal wave system present and the corresponding images are cusped in nature due to velocity bunching. However, due to the smoothing effect of the finite scene coherence time, the ocean image appears relatively smooth and unimodal. Simulations were performed in order to observe the effect of the coherence time upon the principal wave component and to draw comparisons with the actual S A R imagery. The simulation parameters chosen correspond to a monochromatic azimuthally travelling wave (<f> = 0°) with a 150 m wavelength, a = 1/2.5TT (sharply cusped, but not double peaked), and coherence times of r ~> T, r = 47 msec, and T = 14 msec. The simulation results are presented on the right of Fig. 19 arranged from top-to-bottom with decreasing coherence time. The corresponding Fourier transforms are shown in Fig. 20. As the coherence time decreases, we move from a very cusped image with three 82 Figure 19: Comparison of C-SAR with simulation (coherence time). The C - S A R imagery is of waves in ice (upper-left) and of waves in the open ocean (lower-left). The simulated imagery corresponds to a 150 m wave with a = 1/2.57T having coherence times of r » T (upper-right), r = 47 msec (middle-right), and r = 14 msec (lower-right). 83 84 Table VIII: Contrast measures or one- and two-dimensional C - S A R simulations. T/NT r (msec) C i - d C*2-d ~ 0 > T 3.199 2.786 3 47 1.887 1.521 10 14 0.673 0.652 harmonics visible (upper-right A l , A2, A3) to a smooth image with only one harmonic visible (lower-right A l ) . It is apparent that the cases of r ^> T (no coherence time limitation) and r = 14 msec (coherence time limited) are similar in both appearance and in Fourier transform nature to the C - S A R images of waves in ice and waves in the ocean respectively. This suggests that a coherence time on the order of 14 msec is appropriate for C-band airborne S A R imagery of the ocean surface. This estimate is more than an order of magnitude smaller than some of the coherence times calculated based upon decorrelation time scales in Chapter 5. Simulation verifies that coherence time is not a limiting parameter for the images of the same waves moving through ice. The smoothing of the cusped image nature in the ocean case is enhanced by the looks becoming misregistered due to noncoherent reflectivity envelope translation. However, the scenes chosen were at near range so that the R/V parameter is as small as possible and the misregistration is minimized. This effect is discussed further in the following Section. Table VIII lists contrast measures for the simulations in Fig. 19 and the one-dimensional simulations of Appendix D . The two-dimensional simulation contrast measures were obtained by integrating the simulation in range to produce the results shown in Fig. 21. The contrast is measured as in Appendix D . In each case, the contrast is comparable between the one- and two-dimensional simulations for each simulated coherence time. Recalling that the target density used in performing the two-dimensional 85 Figure 21: Range integrated two-dimensional simulations for the three coherence times considered. 86 simulation is rather small, the approach to obtaining this comparison and relating the two-dimensional simulation to the one-dimensional simulation (and to most other SAR ocean simulations) shows how the assumption of diffuse scattering and a large degree of spatial averaging can be misleading. 8.4 T r a n s l a t i n g R e f l e c t i v i t y P a t t e r n s As discussed in earlier chapters, the S A R image of an ocean wave field constitutes the imaging of a translating pattern of microwave reflectivity which is made up of many independent point targets, each with its own dynamic properties. S A R imagery of waves in ice presents an ideal opportunity to observe the effect of the translating reflectivity pattern because coherence time limitations are not an issue for ice imagery. Each look is imaged at a different azimuth time. Also, the available coherent integration time is range dependent. Thus, we expect the relative misregistration between looks to increase for an ocean wave with an azimuthal wavenumber component. We expect to see a general loss in wave contrast at the far range compared with the near range. This effect should be particularly pronounced for the ice-covered case versus the open ocean case due to the absence of a coherence time limitation, and the resultant sharply peaked response functions for certain wave conditions [3,64,77]. C - S A R imagery of waves in ice with an azimuthal wavenumber component was compared with two-dimensional simulated imagery of waves with a large coherence time. The results are presented in Fig. 22. The figure shows the cases of near, mid-, and far range. Note that as the range increases, so does the degree of smear associated with the wave crests. The same effect is present in the ocean scene, but it is masked to a large degree by the finite scene coherence time. The upper-right hand image represents the multi-look image in the absence of noncoherent target motions. In principle, the images corresponding to the near, mid-, and far range cases presented in the figure can have their individual looks re-registered to produce the pattern with the sharpest 87 Figure 22: Comparison of C-SAR with simulation (envelope translation). The comparisons are at near range (top), mid-range, and far range (bottom). The parameters are A = 150 m and a = 1/2.5TT. 88 contrast shown in the upper-right. If the re-registration is made by adjusting the processor focus setting, the contrast, although improved, will be degraded from that which appears in the upper-right image. 8.5 S u m m a r y This Chapter has presented qualitative and quantitative comparisons between simulated and actual S A R imagery. Comparisons were made to address the Bragg scattering nature, the narrow "V" ship wake, the scene coherence time, and the role of the translating reflectivity pattern. The extraction of key statistics from simulated imagery of varying target density has demonstrated that convergence to diffuse scattering statistics is rather slow. The extraction of the same statistics from S E A S A T S A R imagery demonstrates the appropriate behaviour for point and diffuse targets in the test scene, but indicates that the ocean scene is diffuse in nature. This result probably arises due to the diversity of targets over the S E A S A T resolution cell in the ocean scene. On the basis of image appearance alone, relatively small target densities will produce simulated ocean scenes which compare favorably with actual S A R ocean surface imagery, as has been verified with high resolution airborne S A R imagery. One arm of a narrow " V ship wake was simulated as a linear target consisting of a line of discrete point targets in Appendix B . This model results in imagery which compares favorably with actual S A R ship wake imagery. The required testing procedure has been specified. The results indicate a specular-like scattering for this linear feature. As the processor focus is varied, the impulse response is broadened. Thus, it is expected that the linear features associated with the narrow "V" ship wake will be correspondingly broadened. This is the focus dependence expected for an isolated specular-like target. Scenes containing ocean waves with various mean target coherence times were simu- 89 lated and compared with actual S A R imagery of surface gravity waves in an ice-covered region and in the open ocean. Comparison between the ice imagery and simulations with long coherence times is favorable, as is comparison between the open ocean imagery and simulations with a coherence time of about 14 msec. This time is considerably shorter than the decorrelation time scales estimated in Chapter 5 and those conventionally assumed in the relevant literature. The translation of the reflectivity envelope at the wave phase velocity was demonstrated by comparing ice imagery at near, mid-, and far range with simulated imagery having appropriately misregistered looks. The misregistration is range dependent because the available coherent integration time is range dependent, as is the degree of look misregistration. Look re-registration confirms the results predicted in this work. These comparisons confirm the role of the translating reflectivity envelope. 90 9 Conclusions In this thesis, current S A R theory for imaging of the ocean surface has been re- derived and extended. A s well, simulated and actual S A R imagery has been compared to address three key issues in S A R ocean surface imaging. T h e first issue examines the coherent versus noncoherent nature of Bragg scattering. W h i l e the Bragg scattering phenomenon is traditionally dismissed as being a diffuse scattering process, this work defends the assertion that B r a g g scattering is a coherent, specular-like process, although experimental statistical tests designed to address this point were inconclusive for the available S E A S A T S A R test scene. A case w i t h a more advantageous signal-to-noise ratio is required by the analysis methodology. T h e narrow " V " ship wake often seen in S A R imagery of azimuthally travelling ships presents such a case. T h e required analysis has been described. T h e second issue is the magnitude, role, and calculation of the scene coherence time. Comparisons between simulated and actual S A R imagery has suggested a coherence time on the order of 14 msec is appropriate for the L E W E X ocean scene observed by an airborne C - b a n d S A R . T h i s result is more than an order of magnitude smaller than some estimates based u p o n conventional decorrelation time models and an idealized ocean surface wave spectrum. Furthermore, some of the theoretical estimates seem independent of b o t h seastate and ambient environmental conditions—a situation which seems unlikely. It was shown that considerably longer coherence times are required for imaging by a satellite S A R . O n l y very highly organized ocean surface structure, such as a swell system, should be imaged at all by a satellite S A R . T h e third issue is the role of the wave phase velocity and orbital velocity. It has been shown that coherent Doppler modulations by the long wave orbital velocity leads to velocity bunching, which is a well known phenomenon. Velocity bunching simulations in this work compare favorably with actual S A R imagery of waves in an ice-covered region—a situation in which coherence time limitations are not an issue. It was also 91 demonstrated that the phase velocity can only manifest itself as the translation velocity of the mean reflectivity envelope. This envelope will cause noncoherent effects in the S A R image, such as look misregistrations and scanning distortions, if the ratio C /V p is sufficiently large. The conclusions of this work rely heavily upon semi-quantitative comparisons drawn between actual and simulated S A R imagery. The simulation is based upon a S A R ocean scattering model which includes velocity bunching (coherent orbital velocity Doppler modulations) and a reflectivity envelope which translates at the wave phase velocity. The simulation is also unique in that explicit account is taken of each target which constitutes the scene. Simulation of ocean surface phenomena with relatively small target densities compare favorably with actual S A R ocean imagery. This fact strengthens the claim that Bragg scattering is a coherent phenomenon for each Bragg scattering region. The work presented in this thesis makes progress in understanding the underlying physical phenomena, and in defusing controversies surrounding three dominant issues in S A R ocean surface imagery. The work on these issues is by no means complete. However, the importance of target density, scene coherence time, and coherent/noncoherent scene velocities have been demonstrated in this work. The results provide new insights into the very complex issues involved, and should aid in improving the design and performance of imaging radar systems intended to faithfully observe global ocean resources. 92 Glossary along-track velocity The target velocity component parallel to the direction of platform motion. azimuth The S A R (image) dimension which is parallel to the platform velocity vector. Bragg scattering Diffracted beams of electromagnetic radiation resulting from constructive interference of coherent radiation homogeneously scattered from a regular lattice of targets. In the radar/ocean case, the radar wavelength and angle of incidence select the lattice, or ocean wavelength scale, which dominates the scattering. For the S E A S A T SAR, the Bragg wavelength is of order 30 cm, while for the C-SAR, the Bragg wavelength is of order 5 cm. C-Band Microwave frequency band with wavelength of order 5 cm. CCRS Canada Centre for Remote Sensing. C-SAR The C C R S C-band digital S A R and real-time processor on the C C R S Convair 580 aircraft. Capillary wave A wave whose restoring force is dominated by surface tension. Capillary waves generally satisfy A <C 1.73 cm. coherence time-bandwidth product A measure of the degree of coherent processing which a partially coherent signal will support. coherent A property implying that the phase is of relevance. A coherent signal has a deterministic phase function. coherent integration time The time over which a target is coherently observed by the radar. This is equal to the time spent within the antenna pattern. coherent modulation Modulation in which the signal phase is preserved or modified. depth of focus The error in target range which still allows reasonable focussing properties. Generally, one depth of focus corresponds to the focus error n = 2N/TB. dispersion relation The relationship between wave frequency and wavenumber. Doppler bandwidth The Doppler frequency range due to the target moving through the antenna pattern. Doppler frequency The frequency shift 2u/A resulting from reflection from a target moving with radial velocity v. Doppler suppression The loss of target signal due to shifts outside of the system bandwidth resulting from Doppler frequency shifts. FFT Fast Fourier transform. focus error The perturbation of the processor focus setting from nominally perfect focus. HH Horizontal transmit, horizontal receive polarization. 93 J P L Jet Propulsion Laboratory, Pasadena, California, U S A . G S A R Generalized S A R processor. The digital S A R processor developed by M D A . L - B a n d Microwave frequency band with wavelength of order 30 cm. L E W E X Labrador Extreme Waves Experiment. Took place in conjunction with L I M E X in March 1987 off the coast of Newfoundland in the Labrador Sea. L I M E X Labrador Ice Margin Experiment. See Lewex. look May be denned in two senses: first, a physical look is the S A R image resulting from applying a data window to extract a portion of the Doppler spectrum; and second, a statistical look is a statistical measure, usually the M S V R , which is equivalent to the number of effective statistically independent physical looks. MDA MacDonald Dettwiler and Associates, Richmond, B . C . , Canada. M S V R Mean-squared-to-variance ratio, a measure of radiometric resolution. For Gaussian signals, M S V R is a measure of the statistical independence of the data sources. multi-look A n image formed from more than one statistically independent look. noncoherent A property implying that the phase is irrelevant. A noncoherent signal has random or undefined phase. o r b i t a l velocity The speed and direction of travel of a particle or packet of water on the ocean surface as a result of the passage of a wave. p a r t i a l l y coherent A property implying the presence of both a coherent and noncoherent component. phase velocity The speed and direction of travel of a point of constant phase on the surface of a wave. R A E Royal Aircraft Establishment, Farnborough, Hampshire, U K . RMS Root-mean-square averaging operation. r a d i a l acceleration The target acceleration component away from the radar. r a d i a l velocity The target velocity component away from the radar. r a d i o m e t r i c resolution The expected variation of the radar cross-section for each pixel. This is often expressed as the M S V R . range The S A R (image) dimension which is perpendicular to the platform velocity vector. This coordinate is naturally the slant-range, along the radar line-of-sight, but may be converted to ground-range via interpolation. S E A S A T A US ocean remote sensing satellite which orbited the earth for approximately 100 days from July to October in 1978. S T F The system transfer function which is proportional to the Fourier transform of the static target, perfect focus impulse response. 94 significant waveheight May be denned in many ways, but here, we choose four time the rms waveheight. single-look of a multi-look set The result of compressing and detecting only a portion of the available Doppler spectrum. speckle The grainy multiplicative noise phenomenon found in all coherently produced imagery of diffuse targets. surface gravity wave A wave whose restoring force is dominated by gravity. Surface gravity waves generally satisfy A > 1.73 cm. swell A highly organized (generally long) surface gravity wave which has propagated away from its point of generation. time-bandwidth product A measure of signal complexity which indicates the degree of coherent processing which a fully coherent signal will support, and the degree of improvement in resolution as a result of such processing. VV Vertical transmit, vertical receive polarization. velocity bunching The concentration of dispersed targets in a S A R image due to spatially variable radial velocity components. wave frequency n = 27r/T . w wavenumber K = 27r/A. wind-waves A system of ocean waves resulting from local wind generation. 95 References Alpers, W . R . , 1983; "Monte Carlo Simulations for Studying the Relationship Between Ocean Wave and Synthetic Aperture Radar Image Spectra," J . Geophys. Res., 88 (CS), pp 1745-1759. Alpers, W . R . , and C . Bruening, 1986; "On the Relative Importance of MotionRelated Contributions to the S A R Imaging Mechanism of Ocean Surface Waves," IEEE Trans. Geoscience Rem. Sens., GE-24 (6), pp 873-885. Alpers, W . R . , D . B . Ross, and C L . Rufenach, 1981; "On the Detectibility of Ocean Surface Waves by Real and Synthetic Aperture Radars," J. Geophys. Res., 86 (C7), pp 6481-6498. Alpers, W . R . , and C L . Rufenach, 1979; "The Effect of Orbital Motions on Synthetic Aperture Radar Imagery of Ocean Waves," IEEE Trans. Antennas Propagat., AP-27 (5), pp 685-690. Beal, R . C , P.S. DeLeonibus, and I. Katz E d . , 1981; Spaceborne Synthetic Aperture Radar for Oceanography, The Johns Hopkins University Press, Baltimore, 215p. Beal, R . C , D . G . Tilley, and F . M . Monaldo, 1983; "Large- and Small-Scale Spatial Evolution of Digitally Processed Ocean Wave Spectra From S E A S A T Synthetic Aperture Radar," J. Geophys. Res., 88 (CS), pp 1761-1778. Beal, R . C , T . W . Gerling, D . E . Irvine, F . M . Monaldo, and D . G . Tilley, 1986; "Spatial Variations of Ocean Wave Directional Spectra From the Seasat Synthetic Aperture Radar," J. Geophys. Res., 91 (C2), pp 2433-2449. Beal, R . C , F . M . Monaldo, D . G . Tilley, D . E . Irvine, E . J . Walsh, F . C Jackson, D . W . Hancock III, D . E . Hines, R . N . Swift, F.I. Gonzalez, D.R. Lyzenga, and L . F . Zambresky, 1986; "A Comparison of SIR-B Directional Ocean Wave Spectra with Aircraft Scanning Radar Spectra," Science, 232, pp 1531-1535. Bennett, J.R., and I.G. Cumming, 1978; "Digital Techniques for the MultiLook Processing of S A R Data with Application to S E A S A T - A , " 5 Canadian th Symposium on Remote Sensing, Victoria, pp 506-516. Bennett, J.R., and P.R. McConnell, 1980; "Considerations in the Design of Multilook Processors for Image Quality," Proc. 3 SEASAT-SAR Workshop on 'SAR Image Quality,' Held at Frascati, Italy, 11-12 D e c , 1980, pp 5-12. r d Brown, W . M . , 1967; "Synthetic Aperture Radar," IEEE tron. Syst., AES-S Trans. Aerospace Elec- (2), pp 217-229. Brown, W . M . , and C J . Palermo, 1962; "Theory of Coherent Systems," IRE Trans. Military Electronics, pp 187-196. Brown, W . M . , and C J . Palermo, 1962; "System Performance in the Presence of Stochastic Delays," IRE Trans. Inform. Theory, IT-8 (5), pp s206-s214. Dawe, B.R., and S.K. Parashar, 1978; "SAR Imaging of Waves in Ice," Oceans '78 4 Ann. Combined Conf., Washington D C , Sept 6-8, 1978, pp 379-384. th 96 [15] Freeman, N . G . , P.W. Vachon, and C E . Livingstone, 1986; "Airborne S A R imaging of Azimuthally Travelling Ocean Surface Waves: The L E W E X Experimental Plan," Proc. IGARSS '87 Symposium, Ann Arbor, May 19-21, 1987, pp 13931398. Garrett, C . , and J . Smith, 1976; "On the Interaction between Long and Short Surface Waves," J. Physical Oceanography, 6, pp 925-980. Goldfinger, A . D . ; "Estimation of Spectra from Speckled Images," IEEE Aerospace Electron. Syst, AES-18 Trans. (5), pp 675-681. Goodman, J . W . , 1968; Introduction to Fourier Optics, McGraw-Hill, New York, 287p. Harger, R . O . , 1970; Synthetic Aperture Radar Systems; Theory and Design, Academic Press, New York, 240p. Harger, R . O . , 1981; "SAR Ocean Imaging Mechanisms," in Beal et al, pp 41-52. Hartl, Ph., A . Popella, and A . J . Sieber, 1985; "SAR Product Simulation," Proc. Workshop on Thematic Applications of SAR Data, Frascati, Italy, Sept. 1985 (ESA SP-257), pp 65-72. Hasselmann, K , T . P . Barnett, E . Bouws, H . Carlson, D . E . Cartwright, K . Enke, J . A . Ewing, H . Glennop, D . E . Hasselmann, P. Kruseman, A . Meerburg, P. Miiller, D . J . Olbers, K . Richter, W . Sell, and H . Walden, 1973; "Measurements of Wind-Wave Growth and Swell Decay During the Joint North Sea Wave Project (JONSWAP)," Deut. Hydrogr. Z., Suppl. A. Nr. 12.. Hasselmann, K . , R . K . Raney, W . J . Plant, W . Alpers, R . A . Shuchman, D.R. Lyzenga, C . L . Rufenach, and M . J . Tucker, 1985; "Theory of S A R Ocean Wave Imaging: A M A R S E N View," J. Geophys. Res., 90 (CS), pp 4659-4686. Holtzman, J . C , V . S . Frost, J . L . Abott, and V . H . Kaupp, 1978; "Radar Image Simulation," IEEE Trans. Geoscience Elect., GE-16 (4), pp 296-303. Ivanov, A . V . , 1982; "On the Synthetic Aperture Radar Imagery of Ocean Surface Waves," IEEE J. Oceanic Eng., OE-7 (2), pp 96-103. Jain, A . , 1978; "Focussing Effects in Synthetic Aperture Radar Imagery of Ocean Waves," Appl. Phys., 15, pp 323-333. Jain, A . , 1981; "SAR Imaging of Ocean Waves: Theory," IEEE OE-6 (4), pp 130-139. J. Oceanic Eng., Jain, A . and O . H . Shemdin, 1983; "L Band S A R Ocean Wave Observations During Marsen," J. Geophys. Res., 88 (C14), pp 9792-9808. Jakeman, E . , 1980; "Speckle Statistics with a Small Number of Scatterers," SPIE Vol. 243 Applications of Speckle Phenomena, pp 9-19. Jakeman, E . , and P.N. Pusey, 1976; "A Model for Non-Rayleigh Sea Echo," IEEE Trans. Antennas Propagat., AP-24 (6), pp 806-814. Jakeman, E . , and P.N. Pusey, 1977; "Statistics of Non-Rayleigh Microwave Sea Echo," Proc. RADAR '77 Int. Conf., pp 105-109. 97 Jakeman, E . , and P.N. Pusey, 1978; "Significance of K Distributions in Scattering Experiments," Phys. Rev. Let., 40 (9), pp 546-550. Jordan, R . L . , 1980; "The Seasat-A Synthetic Aperture Radar System," J. Oceanic Eng., OE-5 (2), pp 154-164. IEEE Kerr, D . E . , 1951; Propagation of Short Radio Waves, McGraw Hill, New York. Klauder, J.R., A . C . Price, S. Darlington, and W . J . Albersheim, 1960; "The Theory and Design of Chirp Radars," Bell Syst. Tech. J., 39, pp 745-808. Lamb, H . , 1932; Hydrodynamics 6 th ed., Dover Publications, New York, 738 pp. Larson, T . R . , L.I. Moskowitz, and J . W . Wright, 1976; "A Note on S A R Imagery of the Ocean," IEEE Trans. Antennas Propagat., LeBlond, P.H.,and L . A . Mysak, 1978; dam, 602p. AP-24, pp 393-394. Waves in the Ocean, Elsevier, Amster- L i , F . - K . , D . N . Held, J . C . Curlander, and C . Wu, 1985; "Doppler Parameter Estimation for Spaceborne Synthetic-Aperture Radars," IEEE Trans. Geoscience Rem. Sens., GE-23 (1), pp 47-56. Livingstone, C . E . , A . L . Gray, R . K . Hawkins, J . G . Halbertsma, R . A . Deane, and R . B . Olsen, 1987; "CCRS C-band Airborne Radar—System Description and Test Results", submitted to Canadian Symposium on Remote Sensing, Waterloo, Ontario, June 22-25. Lyden, J . D . , 1985; Analysis of SEASAT Revolution 407 Ship Wake Data, EREVI Topic Report 155900-32-T, 17p. Lyden, J . D . , R.R. Hammond, and D.R. Lyzenga, 1987; "Synthetic Aperture Radar Detection of Surface Ship Wakes," submitted to J. Geophys. Res.. Lyzenga, D.R., 1986; "Numerical Simulation of Synthetic Aperture Radar Image Spectra for Ocean Waves," IEEE Trans. Geoscience Rem. Sens., GE-24 (6),pp 863-872. Lyzenga, D.R., A . L . Maffett, and R . A . Shuchman, 1983; "The Contribution of Wedge Scattering to the Radar Cross Section of the Ocean Surface," IEEE Trans. Geoscience Rem. Sens., GE-21 (4), pp 502-505. Monaldo, F . M . , 1984; "Improvement in the Estimation of Dominate Wavenumber and Direction from Spaceborne S A R Image Spectra when Corrected for Ocean Surface Movement," IEEE Trans. Geoscience Rem. Sens., GE-22 (6), pp 603-608. Monaldo, F . M . and D.R. Lyzenga, 1986; "On the Estimation of Wave Slopeand Height-Variance Spectra from S A R Imagery," IEEE Trans. Geoscience Rem. Sens., GE-24 (4), pp 543-551. Moore, R . K . , 1985; "Radar Sensing of the Ocean," IEEE 10 (2), pp 84-112. Neuman, G . and W . J . Pierson Jr., 1966; Principles Prentice-Hall, London, 545 p. 98 J. Oceanic Eng., OE- of Physical Oceanography, North, D . O . , 1943; "Analysis of the Factors Which Determine Signal/Noise Discrimination in Pulsed-Carrier Systems," RCA Lab. Rept. PTR-6C; reprinted in Proc. IEEE, 51 (7), pp 1015-1027. Ouchi, K . , 1986; "Multi-look Processing of Synthetic Aperture Radar Data from Dynamic Ocean Surfaces," Pattern Recognition Letters, 4, PP 305-314. Palermo, C . J . , 1963; Theory of Stochastic Delays, Rep. 2900-383-T Institute of Science and Technology, University of Michigan, A n n Arbor (Also available from University Microfilm Inc., Ann Arbor, Michigan, Document 64-868) 83p. Papoulis, A . , 1965; Probability, Random Variables, and Stochastic Processes, McGraw-Hill Book Co., New York, 583p. Phillips, O . M . , 1980; The Dynamics of the Upper Ocean 2 nd ed., Cambridge University Press, 326p. Phillips, O . M . , 1981; "The Structure of Short Gravity Waves on the Ocean Surface," in Beal et al, pp 24-31. Pierson, W . J . , and L.I. Moskowitz, 1964; "A Proposed Spectral Form for Fully Developed Wind Seas Based on the Similarity Theory of S.A. Kitaigorodskii," J. Geophys. Res., 69 (24), pp 5181-5190. Plant, W . J . and W . C . Keller, 1983; "The Two-Scale Radar Wave Probe and S A R Imagery of the Ocean," J. Geophys. Res., 88 (C14), pp 9776-9784. Pond, S. and G . L . Pickard, 1983; Introductory Dynamical Oceanography, 2 nd ed., Pergamon Press, Oxford, 329 p. Porcello, L . J . , N . G . Masey, R . B . Innes, and J . M . Marks, 1976; "Speckle Reduction in Synthetic Aperture Radars," J. Opt. Soc. Amer., 66 (11), pp 1305-1311. Pravdo, S.H., B . Huneycutt, B . M . Holt, and D . N . Held, 1983; Seasat Aperture Radar Data User's Manual, Synthetic- J P L Publication 82-90. Raney, R . K . , 1971; "Synthetic Aperture Radar Imaging and Moving Targets," IEEE Trans. Aerospace Electron. Syst., AES-7 (3), pp 499-505. Raney, R . K . , 1980; "SAR Response to Partially Coherent Phenomena," Trans. Antennas Propagat., AP-28 IEEE (6), pp 777-787. Raney, R . K . , 1981; "Wave Orbital Velocity, Fade, and S A R Response to Azimuth Waves," IEEE J. Oceanic Eng., OE-6 (4), pp 140-146. Raney, R . K . , 1983; "Transfer Functions for Partially Coherent S A R Systems," IEEE Trans. Aerospace Electron. Syst., AES-19 (5), pp 740-750. Raney, R . K . , 1985; "Theory and Measure of Certain Image Norms in SAR," IEEE Trans. Geoscience Rem. Sens., GE-23 (3), pp 343-348. Raney, R . K . , 1986; "Doppler Properties of Radars in Circular Orbits," Int. J. Rem. Sens., 7 (9), pp 1153-1162. Raney, R . K . , and R . T . Lowry, 1978; "Oceanic Wave Imagery and Wave Spectra Distortions by Synthetic Aperture Radar Imagery," 12 International Symposium on Remote Sensing of Environment, Manila, Phillipines, 20-29 April, 1978, pp 683-702. th 99 Raney, R . K . , and P.W. Vachon, 1987; "SAR Imaging of Ocean Waves from an Airborne Platform: Focus and Tracking Issues (Part I)," submitted to IEEE Trans. Geoscience Rem. Sens.. Raney, R . K . , and G . J . Wessels, 1987; "Spatial Considerations in Speckle Simulation," Proc. IGARSS '87 Symposium, A n n Arbor, 18-21 May, 1987, pp 1341- 1348. Rayleigh, J.W.S., 1896; The Theory of Sound 3 rd ed., MacMillan, London. Rihaczek, A . W . , 1969; Principles of High-Resolution Co., New York, 498p. Radar , McGraw-Hill Book Rice, S.O., 1944; 1945; "Mathematical Analysis of Random Noise," Bell Syst. Tech. J., 23, pp 282-352; 24, pp 46-156. Shemdin, O . H . et al, 1986; TOWARD Field Experiment Interim Report, Vol. I: Data Summary and Early Results, Vol. II: Contributions of Individual Investigators, J P L . Shuchman, R . A . , 1981; "Processing Synthetic Aperture Radar Data of Ocean Waves," in Oceanography from Space, J . F . R . Gower E d . , Plenum, New York, pp 477-496. Shuchman. R . A . , and J.S. Zelenka, 1978; "Processing of Ocean Wave Data from a Synthetic Aperture Radar," Boundary Layer Met., IS, pp 181-191. Swerling, P., 1950; "Probability of Detection for Fluctuating Targets," Rand Report, RM-1217; reprinted in IRE Trans., IT-6 (2), pp 269-308. Swerling, P., 1957; "Detection of Fluctuating Pulsed Signals in the Presence of Noise," IRE Trans. Information Theory, IT-3 (3), pp 175-178. Swift, C T . and L.R. Wilson, 1979; "Synthetic Aperture Radar Imaging of Moving Ocean Waves," IEEE Trans. Antennas Propagat., AP-27 (6), pp 725-729. Tomiyasu, K . , 1978; "Tutorial Review of Synthetic Aperture Radar (SAR) With Applications to Imaging of the Ocean Surface," Proc. IEEE, 66 (5), pp 563-583. Tucker, M . J . , 1985a; "The Imaging of Waves by Satelliteborne Synthetic Aperture Radar: The Effects of Sea-Surface Motion," Int. J. Rem. Sens., 6 (7), pp 1059-1074. Tucker, M . J . , 1985b; "The Decorrelation Time of Microwave Radar Echoes from the Sea Surface," Int. J. Rem. Sens., 6 (7), pp 1075-1089. Ulaby, F . T . , R . K . Moore, and A . K . Fung, 1982; Microwave Remote Sensing Vols. I, II, and III, Addison Wesley Publishing Co., Reading Mass., 2162p. Valenzuela, G . R . , 1978; "Theories for the Interaction of Electromagnetic and Oceanic Waves—A Review," Boundary Layer Met., IS, pp 61-85. Vant, M . R . , 1982; "A Spatially Variant Autofocus Technique for Synthetic Aperture Radar," Proc. Radar '82, IEE Conf. Pub. No. 216, pp 159-163. 100 Vesecky, J . F . , S.L. Durden, M . P . Smith, and D . A . Napolitano, 1984; "Synthetic Aperture Radar Images of Ocean Waves, Theories of Imaging Physics and Experimental Tests," Proceedings of the URSI Commission F Symposium and Workshop, Shoresh, Israel, May 14-23, 1984, pp 137-148. Vesecky, J . F . , and R . H . Stewart, 1982; "The Observation of Ocean Phenomena Using Imagery from the S E A S A T Synthetic Aperture Radar: A n Assessment," J. Geophys. Res., 87 (C5), pp 3397-3430. Wakasugi, K . , N . - F . Kishi, and M . Matsuo, 1986; "A Three-Dimensional Formulation for Synthetic Aperture Radar Images of Ocean Waves in Orbital Motions," IEEE Trans. Geoscience Rem. Sens., GE-24 (5), pp 732-737. Walsh, E . J . , D . W . Hancock III, D . E . Hines, R . N . Swift, and J . F . Scott, 1985; "Directional Wave Spectra Measured with the Surface Contour Radar," J. Physical Oceanography, 15 (5), pp 566-592. Woodward, P . M . , 1953; Probability and Information Theory, With Applications to Radar, The MacMillan Company, New York, 128p. Wright, J . W . , 1966; "Backscattering from Capillary Waves with Application to Sea Clutter," IEEE Trans. Antennas Propagat, AP-14 (6), pp 749-754. Wright, J . W . , 1968; "A New Model for Sea Clutter," IEEE Propagat., AP-16 (2), pp 217-223. Trans. Antennas Zelenka, J.S., 1976; "A Comparison of Continuous and Discrete Mixed Integrator Processors," J. Opt. Soc. Amer., 66 (11), pp 1295-1304. 101 A Simulation Methodology A b s t r a c t A two-dimensional SAR image simulation model developed for this research is discussed in this Appendix. The model is unique in that it explicitly considers the contribution of each target to the final image. System linearity is used to produce the complex image in a reasonably efficient manner. The simulation can vary the target density, models velocity bunching for an azimuthally travelling monochromatic wave, includes a processor focus error, and includes either point target accelerations or scene partial coherence. A.1 I n t r o d u c t i o n There are four main objectives for the S A R simulation model developed for this research. The first objective is to study the role of a variable target density upon SAR image statistics. Second is to study the detectability of a linear image feature in a nominally speckled background. Third is to study the velocity bunching mechanism for the imaging of azimuthally travelling surface gravity waves. The fourth main objective is to study the role of a finite scene coherence time in a velocity bunching model. The SAR simulation considered in this work is seen to be unique because the effect of each point target which constitutes the scene is explicitly considered. Such a simulation philosophy is in contrast with all other S A R simulation models [1,21,24,68,86]. Most simulations contain inherent assumptions of diffuse scattering (dense targets) and heavy post-detection image smoothing (no speckle-like processes). The simulation is two-dimensional, although a one-dimensional simulation could be adequate to address many aspects of the stated objectives. However, a two-dimensional simulation yields imagery for explicit comparison with actual S A R imagery, allows evaluation on an image analysis system, and allows direct application of the statistical testing procedures developed elsewhere in this work. 102 A.2 Method The simulation is designed to be as efficient as possible by taking advantage of the SAR system's linearity up to the point of detection. Thus, the complex "image" or complex look is the superposition of the appropriate complex impulse responses each weighted by the target's reflectivity 14 and intrinsic phase. Subsequently, the complex looks may be detected and summed to produce the real S A R image product. This impulse response superposition procedure nicely circumvents the need to generate new range-Doppler responses for each and every target and then perform the azimuthal compression. (The latter approach, if implemented, would be a very computationally expensive procedure.) All targets subject to the same motion parameters will have the same complex impulse response for each look. Thus, a set of responses which adequately spans the range of motion parameters presented by the scene is pre-calculated and treated as a set of complex look-up tables. A field of static, fully coherent point targets requires only one such impulse response. A dynamic scene may require many. Fig. 23 is a flow chart showing the major steps and components in the simulation sequence. These steps are now each briefly discussed. A.2.1 T h e Simulation Parameter File The simulation parameters are initially read from a disk file in three groups: 1. Scene Parameters • output grid size (generally 128-by-128 pixels) • target density (or set of densities) • target grid type • wavelength of a monochromatic azimuthally travelling wave A 1 4 The reflectivity amplitude is proportional to the square-root of the radar cross-section. 103 slmulaCIon parameter f i l e 3 target parameter generat ion 2 complex response generation \ target parameter f i l e target response summation 5 6 complex output detection and real output complex image real image f i l e f i l e Figure 23: Simulation sequence flow chart. 104 • velocity bunching parameter a which implies the wave amplitude A • mean scene coherence time r c 2. S A R Parameters • platform velocity V • platform height h • radar wavelength A • angle of incidence 7 • available coherent integration time T • Doppler bandwidth B • pulse repetition frequency P R F 3. Processor Parameters • azimuthal processor focus error n • range processor focus error r\ r • number of looks processed N • look overlap £ • interpolation factor e • decimation factor e u d Two S A R systems were simulated in this work: the satellite S E A S A T L-band SAR processed on the G S A R ; and the C C R S C-Band S A R processed on the C-SAR realtime processor. Nominal parameters used in the simulations for these SAR/processor combinations are listed in Table IX. 105 Table IX: Radar parameters relevant to the simulations. S E A S A T refers to the L-band S A R on the S E A S A T satellite and the G S A R processor while, C - S A R refers to the C-band S A R on the C C R S Convair 580 aircraft and the real-time processor. SEASAT C-SAR V (m/s) h (m) 1 Ro (m) A(m) T(s) T •B PRF (Hz) JV e U f, (Hz) A x (m) Pm (m) 7200. 800,000. 20° 850,000. 0.235 2.3 2760. 1600. 4 0.42 130. 3600. ~ 45° 5090. 0.050 1.0 130. 2.57V = 334.1 7 0.12 V2 4. 566. 12.7 24. 1. 10. 33.4 3.89 7. 106 A.2.2 C o m p l e x Impulse Response G e n e r a t i o n A series of complex impulse responses are pre-calculated and stored for each look generated. The number of responses generated depends upon the range of target accelerations and target coherence times presented by the scene. If scene coherence time is incorporated in a simulation run, the situation is assumed to be coherence time-limited so that processor and acceleration defocus are ignored. Up to eleven responses may be generated as required. The complex azimuthal response used in this work is ^*-* ' , V53^"" {" * p = S ( J T , " 1 ) P} •^{-^Jw^h^^^TF ''] }' 2 1 2 where the subscript c denotes the complex impulse response and e = TB8n. <> 110 The focus perturbation in the case of a long coherence time is Sn = r)-ri = r i - — (Hi) or, in the case of coherence time limited imaging, is S "-¥B(B- (112) That is, the mean first-order effects of a finite scene coherence time are modelled as a (possibly large) focus error. Such a focus error will likely violate 8n <C 1—an assumption used in the derivation of the impulse response. However, the expected behaviour of an individual response is consistent with the expected broadening in the face of a finite scene coherence time. Furthermore, the phase structure of this response preserves the system bandwidth and hence, the scene behaviour in the event of a large target density . Various other computationally efficient methods of modelling the finite scene coherence time were investigated including variation of the decimation rate. However, the bandwidth is generally not preserved as required. Short of explicit compression of each target's Doppler response, incorporation of a finite scene coherence time as a large focus error is a reasonable approach. The multi-look image represents the approximation of the finite scene coherence time. The individual looks of the multi-look set behave in a different manner for a coherence time error than for a focus error. 15 15 107 Several nonrestrictive assumptions are used in deriving the above form of the complex impulse response. These assumptions are matched filter processing, a large platform velocity implying no Doppler suppression, perpendicular antenna pointing, and a flat earth. The look-centre-times are calculated by specifying the overlap of the equivalent rectangular width of the look-extraction windows. Then, we find *.•= ( ^ y ^ ) {2i-N-l)Ti l < i < N where £ is the fractional window overlap and Tj = T/y2[N 2 (113) — l) is the look-extraction window width. The impulse response is generated at the sampling frequency / . = PRF — , (114) where P R F is the radar pulse repetition frequency, e is the processor interpolation facu tor, and ed is the processor decimation factor. The interpolation factor is occasionally applied prior to detection to avoid aliasing during the spectrum spreading detection operation, while the decimation factor is used to discard redundant information after look-extraction. A n additional over-sampling factor arises (and is subsequently removed) by stepping through the look-up table. This permits freedom for accurate sub-pixel target placement. The number of points generated for each impulse response is based upon twice the equivalent rectangular width of the detected multi-look impulse response. This could cause many zero contributions in the event of a large focus error. However, the relative look registration is preserved—a critical property to accurately maintain the multi-look system resolution. Representative examples of sets of impulse responses are shown in Fig. 24. The complex range impulse response is generated using the above complex impulse response with N = 1, = 0, T = \[2, B —»• B/N, 108 and 8r) = n (the range processor r o Azimuth Coordinate (m) Figure 24: Representative examples of detected impulse responses, a) The case of the S E A S A T S A R , and b) the case of the C C R S C-band S A R each with n = 0.05. 109 focus error). Sampling is again at the rate f . s The range impulse response has the same resolution as the multi-look azimuthal response and has unity amplitude. It is assumed that the scene does not contribute any range channel focus errors, and that there is no coherence time impact upon the range channel. The two-dimensional complex impulse response is formed as the matrix product of the azimuthal and range responses. A.2.3 Target Parameter Generation The three target distributions currently available in the simulation are a regular grid, a random distribution, and a line of targets at an angle fa to the azimuth. The regular grid is generally used as a test pattern. For each target location generated, the target's radial velocity v, radial acceleration a, and coherence time r are calculated. The velocity and acceleration are based upon the orbital motion of the azimuthally travelling wave. The coherence time is generated independently for each target based upon a Rayleigh distribution of coherence times with a mean of r . The corresponding c focus errors then follow the distribution (115) Focus errors are divided into eleven bins of equal probability ranging from 0.05 to 0.95. The azimuthal target location is shifted by the radial velocity component. The acceleration and coherence time are used to specify the appropriate response to select from the look-up table. The target location and response parameters are stored in a target parameter disk file. This reduces the amount of core computer memory required since the simulation proceeds on a look-by-look basis requiring this information be retained for each look. 110 A.2.4 Response S u m m a t i o n Each target is weighted by its reflectivity and intrinsic phase prior to summation into the complex grid. To date, the reflectivity used has always been unity since we are primarily interested in the dynamic effects. The intrinsic phase is assumed to be uniformly distributed on 0 to 2ir from target-to-target and from look-to-look (even in the event of overlapping look-extraction windows). A.2.5 Complex Output The complex looks may be optionally output to the disk to permit later use for summation with other complex outputs. Such a procedure is possible because of the SAR's linear nature. Also, the complex looks may be output at any intermediate stage in the simulation. This permits, for example, a single run of ten targets per resolution element to provide any smaller target density case. A line of point targets in a nominally diffuse background is generated by making two separate simulation runs and combining the complex outputs of each. The line and the background are each independently simulated. The line is weighted by its relative reflectivity -y/ovei (based upon the mean value of the detected multi-look diffuse background) and each of the complex looks associated with the line are summed with the corresponding complex looks of the diffuse background. A.2.6 Detection and Real Output As an option, the detected real looks may be output to the disk along with the complex looks. The real outputs are generally scaled to between 0 and 255 in preparation for interpretation on an image analysis system. Ill A.2.7 O u t p u t to O t h e r A n a l y s i s Systems The simulation products may be output to any other systems and the detected products treated in the same manner as actual S A R imagery. A.3 I m p l e m e n t a t i o n The simulation was written in Fortran-77 and currently resides on the V A X 11750 at the R A D A R S A T Project Office (RPO) in Ottawa, Ontario. The R P O system runs the V A X / V M S operating system and has an interactive Dipix Aries III image analysis work station. Original simulation development was on the V A X 11-750 at the Satellite Oceanography Lab at The University of British Columbia (Department of Oceanography) in Vancouver, British Columbia. A.4 F u t u r e E x t e n s i o n s Many enhancements could be incorporated relatively simply into the framework of the simulation developed thus far. For example, it is straightforward to include Doppler suppression for the case of an airborne S A R , range dependent radar crosssection modulation mechanisms so that off-azimuthally travelling waves may be imaged, a full spectrum of waves or at least several wave components, and variable/selectable target reflectivities. Such extensions would increase the simulation's flexibility and applications. A.5 C o n c l u s i o n s The simulation model presented in this Appendix is unique in its explicit treatment of each target which constitutes the scene being imaged. This approach allows point and diffuse targets (and all intermediate cases) to be simulated in exactly the same manner—that is, via explicit summation of the appropriately weighted and shifted complex impulse responses. 112 This model is applied to the simulation of several S A R image phenomena in Appendix B , and provides some rather insightful results. There are potential future applications as well. These could include, the a priori assessment of arbitrary S A R systems for ocean surface imaging, and the investigation of corrections for the system transfer function in the event of non-diffuse scattering. 113 B Simulation of S A R Image Phenomena A b s t r a c t Varying target densities and a linear target, consisting of a line of point targets, are investigated using the two-dimensional simulation model of Appendix A and the statistical analysis methodology of Chapter 7. It is shown that significantly different image statistics due to varying target densities are not necessarily accompanied by a noticeably different image appearance—a fact which suggests why ocean scattering is traditionally dismissed as being diffuse in nature. A linear target composed of a line of discrete point targets is proposed as a model for a narrow "V" ship wake. The results of the statistical analysis are consistent with predictions if the linear target is oriented perpendicular to the channel in which the processor focus is perturbed. B . l I n t r o d u c t i o n As an aid in understanding how a S A R images the ocean surface, a simulation model has been developed. This model is unique in that careful account is taken of the scattering nature; there are no a priori assumptions as to the statistical behaviour of each resolution element since each scattering centre is explicitly considered. In this Appendix, two experiments with this simulation model are outlined. The first addresses the scattering nature associated with Bragg resonant scattering while the second considers the imaging of narrow "V" ship wakes. B.2 V a r y i n g P o i n t T a r g e t D e n s i t y It is relatively well established that a dominant ocean surface scattering mechanism is Bragg resonant scattering [89]. It is also true that if a single coherent scattering centre dominates the return from a resolution cell, then the scattering is rather more coherent than diffuse. A single coherent patch of Bragg waves for a given resolution cell likewise results in a coherent specular-like return [64]. Thus, Bragg scattering is not necessarily equivalent to diffuse scattering, a position which is generally not taken [23]. Such specular occurrences impact both the way in which S A R imagery is formed, and quantitative measures derived from such S A R imagery. This experiment with the simulation model is intended to help establish the relative 114 Table X : Scene simulation parameters (target density study). Target densities and focus perturbations for a 128-by-128 pixel subscene with nominal S E A S A T S A R parameters. number of focus depth of focus case log (density) density targets error per look multi-look 1 2 -1.0 -0.8 0.10 0.16 460 729 0 0.86 0 3.45 1155 0% 1% 2% 3 -0.6 0.25 1.73 6.90 4 -0.4 0.40 1831 3% -0.2 0.0 0.63 1.00 2902 4600 5% 2.59 4.31 10.35 17.25 5 6 7 0.2 1.58 7290 8 9 10 11 0.4 0.6 0.8 1.0 2.51 3.98 6.31 11,555 18,313 29,024 10.00 46,000 density of ensembles of point targets for which specular or diffuse scattering properties dominate S A R imagery. The number of targets is varied in a controlled manner based upon the target density per iV-look, zero focus-error resolution cell size. Each target is randomly positioned across the simulation region. The appropriate complex impulse response is weighted by its cross-section and its intrinsic phase, and is summed to produce N complex looks. In this experiment, the target cross-sections are all chosen to be unity. This means that the results are best-case estimates for convergence towards diffuse scattering. (A distribution of target cross-sections would cause the scattering to retain aspects of a specular-like nature for larger scatterer densities.) The intrinsic phase is assumed to be distributed uniformly on 0 to 2TT. The target densities considered range from 0.1 to 10 targets per resolution cell and the azimuthal processor focus is varied through roughly 0, 1, 2, 3, and 5 depth of foci in each individual look for nominal S E A S A T S A R parameters. Specific simulation values considered are listed in Table X and some representative outputs from the simulation model are shown in Fig. 25 and in Fig. 26. Plots of the image statistics considered, 115 Figure 25: Variable target density simulation with n = 0.00. The densities range from 0.1 (a) to 10 (k) targets per resolution element according to Table X . 116 Figure 26: Variable target density simulation with n — 0.03. The densities range from 0.1 (a) to 10 (k) targets per resolution element according to Table X . 117 specifically the inter-look cross-correlation coefficient, the scaling of peak-heights as looks are summed, and the image spectral width as a function of processor focus error are given in Fig. 27. The ideal behaviors are illustrated by the diffuse (dashed line) and specular (dotted line) cases included in Fig. 27b. It is apparent that the asymptotic behaviour of these cases is correct. However, it is equally apparent that these asymptotes are approached rather slowly. That is, even the extreme cases represented by 0.1 (nominally specular) and 10 (nominally diffuse) targets per resolution cell have not completely reached their limiting behaviors. As may be judged from Fig. 25 and from Fig. 26, there is not a significant difference visually between the appearance of say, 0.63 targets per resolution cell (image e) and the denser cases considered in those figures. It is evident that significantly different image statistics are not necessarily accompanied by noticeably different image appearance. This suggests why even rather specular-like cases could be construed as being purely diffuse in nature on the basis of image appearance alone. Perhaps ocean reflectivity, dominated by Bragg resonant scattering, is often dismissed as being diffuse in nature since the imagery "looks" like it is derived from a diffuse scattering mechanism. Bearing this assertion in mind, actual S A R ocean surface imagery may be subjected to the same analysis as has been done to the controlled simulated images to assess the true nature of Bragg resonant scattering. B.3 L i n e a r T a r g e t s The imaging of ship wakes by S A R could prove to be a valuable tool in assessing S A R ocean imaging mechanisms. It is well established that ship wakes often manifest themselves as a bright narrow "V" pattern in S A R imagery, rather than as the classic Kelvin wake pattern, which is, in fact, rarely observed in S A R imagery. It has been established on the basis of group velocity versus ship velocity arguments and the opening angle of the "V" that this feature is a result of reflectivity from Bragg scale waves 118 10 ' Target Density a.oo 10 0.01 0.02 0.03 Focus Error o.o* 0.05 Figure 27: Derived statistics for variable target density, a) Azimuthal image bandwidth as a function of target density, b) azimuthal image bandwidth as a function of processor focus error, c) peak-height scaling as a function of target density, and d) correlation as a function of target density. 119 produced by the passage of the ship [41]. This phenomena seems exclusively associated with nominally azimuthally travelling ships. In this experiment with the simulation model, one arm of the narrow "V" ship wake is modeled as a linear target consisting of a line of coherent point targets. There are two objectives: first, to assess the detectability of such a well defined structure in a nominally speckled background subject to the chosen analysis methodology; and second, to assess the applicability of such a model and analysis to the narrow "V" wake problem. In this experiment the linear target is composed of individual targets separated by a single resolution element. The complex linear target in zero background is first generated with unity cross-section and randomly distributed intrinsic phase at angles ranging from <p t — 0° (along azimuth) to <pi = 90° (along range) in 10° increments. Subsequently, the linear target is weighted and summed into the nominally diffuse complex background composed of 100 unity cross-section targets per resolution element. The weighting is based upon the single iV-look resolution element being cr i dB above re the mean background level, that is, the mean value of the nominally diffuse, multi-look background scene. The cross-sections considered range from a i = 0 dB to a \ = 10 dB re re in 1 dB increments. The simulation is again subjected to the processor focus errors of Table X . Examples of these simulation outputs are given in Fig. 28 and in Fig. 29. Plots of the statistics for the linear target in zero background, that is, the infinite signal-to-noise ratio case, are shown in Fig. 30. The plot of azimuthal image bandwidth versus linear target angle in Fig. 30b appears to be tending towards zero bandwidth at <j>i = 0° {i.e. azimuthally oriented). However, for small angles this trend disappears since the bandwidth has become a measure of the structure along the linear target. Note that the separation between the curves (Fig. 30a) and the sensitivity to focus (Fig. 30b) are reduced as the linear target becomes azimuthally oriented. The sensitivity is best, and in fact compares with the 0.1 target per resolution cell case of Fig. 27b, for the 120 Figure 28: Variable focus error linear target simulations. The case of <f> = 4 0 ° , a = 8 dB, and the five focus errors considered ranging from n = 0.00 (a) ton = 0.05 (e) according to Table X . t r e l range oriented linear target. This lack of sensitivity for the azimuthally oriented linear target is unfortunate since all observed cases of narrow "V" wakes are azimuthally oriented. However, in the event of a static target, as is very nearly the case with the Bragg scale waves, the range and azimuthal S A R channels are logically interchangeable. Thus, focus sensitivity tests may be performed using the range channel with no loss in generality. In this case, the sensitivity to focus will correspond to the most sensitive case presented (given an equivalent time-bandwidth product). The peak-height scaling and correlation measures are not affected by such considerations since they are based upon nominally zero focuserror cases. Fig. 31 shows the analysis results for the linear target of varying orientations, crosssections, and azimuthal focus-errors in the case of a nominally diffuse background. It is apparent that although the results are tending towards those of Fig. 30 as the crosssection (T i re increases, they do so rather slowly. In the case of the azimuthally oriented 121 Figure 29: Variable relative cross-section linear target simulations. The case of fa = 4 0 ° , rj = 0.00, and c i ranging from 0 dB (a) to 8 dB (i) in 1 dB increments. re 122 20 40 Angle (°) 60 80 0.60 0.01 C.02 0.03 0.04 0.05 Focus Error Angle (°) Figure 30: Derived statistics for linear targets in zero background, a) Azimuthal image bandwidth as a function of linear target angle, b) azimuthal image bandwidth as a function of processor focus error, c) peak-height scaling as a function of linear target angle, and d) correlation as a function of linear target angle. 123 Table XI: Scene simulation parameters (linear target study). Target densities along the azimuthally oriented linear target and range channel focus perturbations for a 128-by-128 pixel subscene with nominal S E A S A T S A R parameters. case log (density) density error depth focus 1 -0.6 0.25 0% 0 2 -0.3 0.0 0.3 0.50 1.06 0.6 4.00 0.5% 1.0% 1.5% 2.5% 3 4 5 1.00 2.00 2.12 3.18 5.30 linear target, Fig. 31a(i) shows the image bandwidth measure decreasing with increasing cross-section. This is due to the structure along the linear target dominating the width measure as the linear target emerges above the background speckle pattern. In the case of the range oriented linear target, Fig. 31a(iv) shows that the width measure is responding to the linear target broadening as the processor focus is perturbed. Fig. 31b and Fig. 31c further illustrate the slow tendency towards the behaviour of the linear target in zero background as shown in Fig. 30a and in Fig. 30b. The peak-height scaling and correlation measures of Fig. 31d and of Fig. 31e are essentially independent of linear target orientation, but again, are only slowly approaching their anticipated specular-like behaviour as the linear target cross-section increases. With these results in mind, two further cases were simulated to observe the sensitivity to target density within the linear target. The linear target is fixed at 4>i — 0° (azimuthally oriented) and the cross-section set to o~ \ = 8 dB. The target densities ie considered range from 0.25 to 4 targets per resolution element and the range and azimuthal processor focus is varied through roughly 0, 1, 2, 3, and 5 depth of foci. Specific values considered are listed in Table X I . Note that the smaller time-bandwidth product associated with the range channel, coupled with one-look range processing, imposes a stricter tolerance on the focus error than in the range channel. Some representative outputs from the simulation model are shown in Fig. 32. Fig. 33 shows the analysis results for the linear target of varying point target den- 124 0 2 4 6 Cross-section 8 10 0 (dB) 2 4 6 8 Cross-section 10 (dB) Figure 31: Derived statistics for linear targets in a diffuse background, a) Azimuthal image bandwidth as a function of linear target cross-section for various target angles at various azimuthal processor focus errors, (continued) 125 Angle (°) Angle ( ° ) Figure 31: Continued, b) Azimuthal image bandwidth as a function of target angle for various azimuthal processor focus errors at various target cross-sections, (continued) 126 Figure 31: Continued, c) Azimuthal image bandwidth as a function of azimuthal processor focus error for various target angles at various target cross-sections, (continued) 127 Figure 31: Continued, d) Peak-height scaling as a function of linear target angle, and e) correlation as a function of linear target angle. sity with azimuthal and. range processor focus errors. It is of note that there is almost no focus error sensitivity along the linear target (i.e. in azimuth) while there is high sensitivity in range. However, the sensitivity to a range focus error is essentially independent of the actual point target density within the linear target. This is expected since in the case of a focus error, the linear target is being subjected to a broader point spread function. Then, as long as the individual target's cross-sections are large relative to the mean speckle-noise level, the measured spectral width in range will decrease with increasing focus error. Bearing these results in mind, the S A R imagery of narrow "V" ship wakes may be subjected to the same analysis methodology, that is, introduction of a focus error in range followed by a rangeward image bandwidth measurement. A similar pattern of responses would suggest that the Bragg resonant scattering associated with narrow " V ship wakes is the result of specular-like scattering. This result could be generalized to any ocean surface Bragg-scale structure since an identical dispersion relationship is 128 a .v b c Figure 32: Variable point target density linear target simulations. The case of varying point target density with I) no processor focus errors, II) azimuthal processor focus errors ranging from n = 0.00 (a) to 77 = 0.05 (e) according to Table X . (continued) 129 Figure 32: Continued. Ill) Range processor focus errors ranging from n = 0.00 (a) to n = 0.025 (e) according to Table XI. r r obeyed for any ocean surface roughness of this scale. B.4 Conclusions This Appendix has presented simulation results and corresponding statistics for two SAR image phenomena. The first phenomenon is the case of a varying target density. Simulated target densities ranging from 0.1 to 10. targets per resolution element have demonstrated that: 1. A target density of 0.1 targets per resolution element (nominally a point target) does not behave in a completely specular manner. 2. A target density of 10. targets per resolution element (nominally a diffuse target and well above the classical guideline of 5 [34]) does not behave in a completely diffuse manner. 130 Figure 33: Derived statistics for linear target of varying point target density, a) Azimuthal image bandwidth as a function of point target density for the various azimuthal processor focus errors, b) azimuthal image bandwidth as a function of azimuthal processor focus error for the various target densities, c) peak-height scaling as a function of point target density, and d) correlation as a function of point target density in a nominally diffuse background, (continued) 131 o CD ao • o- 10 ' Target Density O 10 | I I I 0.000 I 0.005 0.010 0.015 0.D20 0.025 Focus Error Figure 33: Continued, e) Rangeward image bandwidth as a function of point target density for the various range processor focus errors, and f) rangeward image bandwidth as a function of range processor focus error for the various target densities in a nominally diffuse background. 3. Significantly different image statistics due to different target densities are not necessarily accompanied by noticeably different image appearance. The second phenomenon is the case of a linear target consisting of a line of discrete point targets. This target is proposed as a model for one arm of a narrow "V" ship wake. The methodology to perform a comparison between this model and actual S A R imagery of a ship wake exists, but thus far, appropriate S A R data has been unavailable. It is also shown that the range channel focus, rather than the azimuthal channel focus, should be perturbed to maximize the sensitivity of the statistical measures. 132 C A c t u a l S A R Image D a t a Sets A b s t r a c t The two types of actual SAR imagery used in the comparisons in this work are discussed in this Appendix: the first, is satelliteborne SAR imagery from the SEASAT SAR processed to four complex looks on the GSAR processor; and the second, is airborne SAR imagery from the CCRS C-band SAR processed to seven looks on the real-time processor. Difficulties in processing the data are discussed, particularly subtleties in processing the complex looks. C l S E A S A T S A R The L-band S A R on the S E A S A T satellite provides the most extensive satellite S A R data set available to date (and available for the next several years). The S E A S A T satellite system became operational in July 1978 and ceased operation after about 100 days due to an electrical failure. Despite such a short operational time span, and requirements of direct transmission of raw signal data to one of five ground receiving stations, the available S E A S A T S A R data set is rather extensive [59]. C.1.1 S E A S A T S A R Processing All of the S E A S A T S A R imagery used in this work was processed on an M D A G S A R processor at either M D A or C C R S . The data is required as individual detected looks for the prescribed statistical tests. Unfortunately, the G S A R cannot directly provide such a data product. Instead, the processor must be halted prior to the detection and look-sum module, and each of the looks extracted in their complex form. The looks may be subsequently registered, interpolated, and detected to provide the individual looks suitable for the testing procedure. The look-detection procedure was performed upon 128-by-128 (post interpolation) pixel sized subscenes. The interpolation requires some care due to the GSAR's rangesegment mode of operation [9]. The data is sectioned into segments of range data, and a single compression filter matched to the centre of the range segment is used to 133 compress all of the azimuthal lines in the segment. Focus errors resulting from range changes towards the edges of the segments lead to several problems in the resulting imagery. The G S A R tolerates broadening of the individual responses due to the focus error, and it corrects azimuthal shifts due to the focus errors by physically shifting the azimuthal lines back into alignment. The G S A R accounts for large phase rotations across range segment seams due to a large net focus error by using two interpolation procedures: first, a complex cubic convolution resampling is used within the segments; and second, a nearest neighbor resampling is used near the segment seams. This procedure is used because the phase rotation across the seam, in conjunction with the complex cubic convolution resampling, results in a loss of signal energy along the seam. This interpolation procedure was duplicated prior to detecting the extracted subscenes. These phase rotations must be considered when dealing with complex imagery produced by a processor operating in such a mode. C.1.2 Selected S E A S A T Scenes Three S E A S A T S A R scenes were chosen for this work: 1. The Irish Sea off Ballyquinton Point (orbit 633, provided by R A E ) contains examples of various ships and their wakes. The data was processed through to a multi-look image, and is only used for qualitative comparisons with linear target simulations. 2. Goldstone, California (orbit 882, provided by CCRS) contains many features of interest. A corner reflector array provides examples of point targets while regions of desert provide examples of diffuse targets. The scene was chosen to supply control data for the statistical analysis work. 3. Duck Island, North Carolina (orbit 1339, provided by C C R S ) , known as Duck-X, has been subjected to intense study because of its ocean wave imagery [6,7]. This 134 scene was chosen as the test scene for the analysis methodology. Scenes 2 and 3 were each processedfivetimes with focus errors of 0%, 1%, 2%, 3%, and 5%. C.2 C-SAR The C-band SAR on the CCRS Convair 580 marked itsfirstoperational use during the recent LIMEX/LEWEX experiments which took place off the east coast of Newfoundland over the Labrador Sea [15]. One objective of this experiment was to image ocean waves as they propagate from the open ocean into ice covered regions. C.2.1 C-SAR Processing The C-SAR has a real-time processor which produces a seven look, square-root mapped image. The C-SAR system can also store the raw signal data on high density digital tape for subsequent processing on the GSAR. Unfortunately, this system was not enabled for the scenes of interest for this work. Thus, the individual looks are not available from this data set. C.2.2 Selected C - S A R Scenes Two scenes were selected from the LIMEX/LEWEX data set: 1. 87-03-13 (Line 10 Pass 8) contains examples of open ocean waves propagating into an ice covered region. The waves are at about 45° to the azimuth. 2. 87-03-14 (Line 10 Pass 1) contains an example of nearly azimuthally oriented waves propagating through an ice covered region. These scenes were processed in narrow-swath mode, meaning that the near scene edge is at an incidence angle of 45°. A 16 km ground range swath is imaged. 135 C.3 C o n c l u s i o n s This Appendix has discussed the preprocessing of and indicated the scenes chosen from the actual SAR image data sets considered in this work. Some care is required in dealing with the complex SEASAT imagery due to the range-segment GSAR processing. Three SEASAT scenes were chosen: one contains examples of narrow "V" ship wakes; another contains examples of static point and diffuse targets; and the last contains an arbitrary ocean scene showing ocean waves. The C-SAR imagery is more straightforward to deal with. Two scene were chosen which show ocean surface waves propagating through an ice-covered region. 136 D One-Dimensional Velocity Bunching Simulations A b s t r a c t This Appendix presents an idealized one-dimensional velocity bunching simulation model designed to study the contrast of a scene produced by velocity bunching. The simulation is over one cycle of an azimuthally travelling monochromatic wave and includes acceleration defocus, Doppler suppression, and finite scene coherence time. It is demonstrated that acceleration defocus and Doppler suppression may be important in some cases, and oppose the effects of velocity bunching. Furthermore, a finite scene coherence time can smooth the cusped wave image and can obliterate the possible doublepeaked nature of some images. D.l I n t r o d u c t i o n The objective of this simulation work is to systematically study the roles of target acceleration, Doppler suppression, and a finite scene coherence time upon the velocity bunching wave imaging mechanism. The simulation performs a pseudo-convolution of the detected impulse response over a dense, uniform array of point targets. We term the procedure a pseudo-convolution because the appropriate impulse response varies from target-to-target dependent upon the motion parameters. Thus, fast convolution techniques may not be used. This simulation procedure, while straightforward to implement, has several inherent assumptions which could prove to be limiting. Only the mean contrast is considered so that a large target density is coupled with heavy averaging. The one-dimensional procedure is related to the two-dimensional point-by-point simulation developed in Appendix A by integrating the two-dimensional simulation in range (or more generally in the direction of the wave crests). This removes the effects of the target density leaving only the mean contrast function (providing the area integrated over is large). The one-dimensional simulation procedure developed has similarities to the one- dimensional Monte-Carlo simulations of Alpers [l]. However, there are departures from Alpers' work in several key areas. One departure is that we treat only a single cycle of a monochromatic wave rather than a full wave spectrum. While this simulation is 137 not restricted to monochromatic waves, this approach permits the explicit assessment of the roles of specific parameters to velocity bunching. On the other hand, Alpers was interested in studying the modulation transfer function which relates the radar image to the ocean wave field. Another departure from Alpers' work is that we include a finite scene coherence time while Alpers did not. Alpers has since recognized the role of the scene coherence time and it has been included in his subsequent two-dimensional simulation work [2]. D.2 Method The one-dimensional simulation proceeds by generating a monochromatic azimuthally travelling surface gravity wave consisting of a large number of point targets uniformly spaced at locations tj and each having unity cross-section. Associated with each point target is a specific radial velocity component Vj and radial acceleration component aj based upon the wave's orbital motion, and an average scene coherence time r . The c appropriate (expected) detected impulse response g {t) is generated for each target 3 md along the wave. Included in the impulse response are the azimuthal shift, acceleration defocus, Doppler suppression, and finite scene coherence time. Each impulse response contribution is summed to produce the final S A R image. The output image is 0(*) = £*iLi(*i-*). (lie) 3 where each target which contributes energy to the output at time t is included in the sum (including contributions from targets on adjacent wavelengths). The impulse response of the j target is given by th B 2 • exp { —27T N 2 + {TBTJ /2) j r 2 + {T/T ) 138 C 2 / . x -.2' r+ 1 - | W (117) Table XII: One-dimensional simulation parameters for C - S A R . a A a D T/NT T (msec) Fig. 34 C a (i) (ii) (iii) b (i) (ii) (iii) (iv) c (i) (ii) (iii) (iv) d (i) (ii) (iii) (iv) (v) (vi) (vii) (viii) (ix) where rjj = ajRo/V 2 0.90 1/27T 1.34 1/TT 2.69 1/TT 2.69 1/TT 2.69 1/TT 2.69 1/TT 2.69 1/2.57T 1.08 1/2.57T 1.08 1/2.57T 1.08 1/2.57T 1.08 1/2.57T 1.08 1/2.57T 1.08 1/2.57T 1.08 1/2.57T 1.08 1/2.57T 1.08 1/2.57T 1.08 1/2.57T 1.08 1/2.57T 1.08 1/2.57T 1.08 1/37T and tuj = VJRQ/V . 2 no no no no no no no yes no no no yes yes yes no no yes no no yes yes yes no no no no no no no no no no no no no no no no no no ~ 0 ~ 0 ~ 0 ~ 0 ~ 0 ~ 0 ~0 ~ 0 ~ 0 ~ 0 ~0 ~ 0 1 2 3 4 5 10 15 20 > T > T > T > T > T > T > T > r » T > r > T > T 140 70 47 35 28 14 9 7 1.858 7.070 3.326 3.326 2.729 2.947 2.431 3.199 2.825 3.321 2.929 3.199 2.845 2.308 1.887 1.575 1.335 0.673 0.364 0.190 The impulse response is augmented by selec- tively enabling the acceleration defocus contribution rjj and the Doppler suppression contribution (the first exponentiation) in order to directly assess the effects of these terms. D.3 Results Simulations using the radar parameters of Appendix A were performed. The results are summarized in Table's XII and XIII and in Fig.'s 34, 35, and 36. For each configuration, simulations were performed to demonstrate the role of four parameters: first, the role of a showing one-to-one, critical, and many-to-one mapping regimes; second, the role of acceleration defocus and Doppler suppression in the many-to-one mapping case; third, the role of acceleration defocus and Doppler suppression in a near 139 Figure 34: Results of one-dimensional simulations for C-SAR. 140 10 ( T / N r ) Figure 35: Degradation in contrast as a function of coherence time for C-SAR. Table XIII: One-dimensional simu ation parameters for S E A S A T . T/Nr T (msec) a C a A Fig. 36 1.271 0.23 no ~ 0 1/3* » T a (i) 2.519 ~ 0 0.35 no » T 1/27T (ii) 1.673 ~ 0 1/TT 0.70 no > T (iii) 1.673 ~ 0 0.70 no » T 1/TT b (i) 1.532 ~ 0 » T 1/TT 0.70 yes (ii) 1.742 ~ 0 > T c (i) 1/2.5* 0.28 no 1.567 ~ 0 » T (ii) 1/2.57T 0.28 yes d (i) (ii) (iii) (iv) 1/2.5* 1/2.5* 1/2.5* 1/2.5* 0.28 0.28 0.28 0.28 yes yes yes yes 141 ~ 0 0.94 1.88 9.40 » T 575 288 1.567 1.294 0.898 58 0.001 Figure 36: Results of one-dimensional simulations for S E A S A T . 142 critical one-to-one mapping case; and fourth, the role of a finite scene coherence time. For each simulation, a wavelength of 150 m was chosen. As well, a = 1/2.5* was chosen for the one-to-one mapping case. The contrast measure used in this work is C = J p e a k ~ J m e a n , (118) -^mean where the J's are the indicated intensity values. Note that it is not relevant to directly compare the contrast of the double peaked cases to the single peaked cases. D.4 Conclusions Using the one-dimensional velocity bunching simulation we have demonstrated that: 1. Target acceleration-induced defocus will smooth a cusped peak in a velocity bunching model. 2. Doppler suppression will reduce the maximum intensity, but may not reduce the contrast. Doppler suppression is only important if the platform velocity is relatively slow. Thus, its impact will be more apparent (and critical) in airborne S A R ocean imagery than in similar imagery from a satellite S A R . 3. A finite scene coherence time will smooth a cusped peak and may completely obliterate a double peak. We have also judged the parameters which are relevant and should be included in the two-dimensional simulation. A coherence time of 14 msec is acceptable for C SAR, but 60 msec may be too short for the imaging of azimuthally travelling waves by SEASAT. 143 E E s t i m a t i o n of the System Transfer Function Abstract The system transfer function (STF) is required to correct wave spectra derived from SAR ocean imagery. A new approach to estimating the STF is to defocus a portion of the wave image to simulate a large, diffuse target rather than to use a different nominally smooth diffuse scene. The concept is shown to work well with a 5% focus error for a ~ 200 m azimuthal wave component in SEASAT SAR imagery. E.l Introduction Ocean wave spectra derived from S A R imagery require correction for the system transfer function (STF) to allow reliable extraction of ocean wave properties [6]. The procedure involves division of the image spectrum by a replica of the spectrum of an impulse response. The impulse response spectrum, the spectrum of the perfect focus processor response to a stationary point target, has been shown to be proportional to the speckle spectrum, the spectrum of the response to a nominally diffuse target [63]. This normalization procedure "whitens" the effect of the speckle noise which may then be removed by subtraction of an easily calculated bias value [17]. The usual procedure for estimating the S T F involves calculation of the spectrum of a nominally diffuse scene [6,7,8,45,46]. Unfortunately, this requirement could lead to complications in that a known, large, diffuse target is required. That target should be from the same S A R , recorded at the same ground station, and the imagery produced on the same processor. Furthermore, the example of a diffuse target should be gathered at about the same instant in time to ensure reliable estimation of the S T F . In the past, the diffuse target constraint has been satisfied by a fetch-limited body of water. Such a scene is difficult to find for open ocean cases. E.2 Method A new approach to the calculation of the S T F hinges upon the fact that as the processor focus is perturbed, the scene contrast is lost at the processor output. The 144 output becomes diffuse in nature—that is, a uniform speckle pattern is created. Thus, the scene in use for calculation of the ocean wave spectrum could be reprocessed with a focus error to produce the uniform speckle pattern, the output of a nominally diffuse target. This approach to estimating the S T F has the advantage that it is not necessary to find an example of a diffuse target which satisfies the previously mentioned constraints. As long as the focus perturbation is large enough, all the contrast in the processor output will be obliterated, and the S T F may be reliably calculated. The necessary focus error can be estimated from the equation for focus limited resolution, which is V [TB)r] P - = B 2 - . ( 1 , 1 9 ) If the azimuthal resolution cell dimension is of the order of the azimuthal wavelength component imaged we find, for S E A S A T S A R parameters and a wavelength of 200 m, that n ~ 0.024 or a 2.4% error in the processor's azimuthal F M rate will suffice. E.3 Example This procedure was performed upon a sample of the azimuthal spectrum of S E A S A T SAR imagery from orbit 1339. The imagery was available with azimuthal focus errors of 0% and 5%. Fig. 37 shows the results of using a smoothed version of the azimuthal spectrum with a 5% focus error to normalize the spectrum calculated from the 0% focus error case (solid line). This normalization compares well with the normalization based upon a nominally diffuse section of desert scene from the Goldstone California orbit 882 scene (dashed line). The differences between the two spectra in Fig. 37 are small, and may not be statistically significant. E.4 Conclusions This approach to the calculation of the S T F would seem to be feasible. The disadvantage is that the scene must be processed twice. However, only a relatively small 145 o CL ' CO I 1 0.00 I I I I I I I I I I I I I I I I 50.00 I I I I I I I I I I 100.00 I M I I I I I I I 150.00 Frequency I I I I I I I I I I 200.00 250.00 (Hz) Figure 37: Azimuthal spectrum corrected for S T F . The solid line is the result of using the 5% focus error case for normalization while the dashed line is the result of using a nominally diffuse region of the Goldstone scene. The azimuthal wavenumber coordinate in cycles per meter is found through division by V. The peak corresponds to about a 220 m wavelength. 146 region need be reprocessed with a focus error, or perhaps a post-compression refocussing scheme could be used if the complex image data is available [83]. Note that the reprocessed version must have a focus perturbation in both the azimuth and the range to remove the S T F effects in both image dimensions. 147
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Synthetic-aperture radar imaging of the ocean surface...
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Synthetic-aperture radar imaging of the ocean surface : theoretical considerations, and experiments with.. Vachon, Paris W. 1987
pdf
Page Metadata
Item Metadata
Title | Synthetic-aperture radar imaging of the ocean surface : theoretical considerations, and experiments with simulated and actual SAR imagery |
Creator |
Vachon, Paris W. |
Publisher | University of British Columbia |
Date | 1987 |
Date Issued | 2010-08-20T17:41:26Z |
Description | Three key areas of controversy in synthetic-aperture radar (SAR) imaging of ocean surface waves are considered: first, the nature of Bragg scattering; second, the role, magnitude, and calculation of the scene coherence time; and third, the relevant ocean wave velocities for coherent Doppler modulations. This work begins with a re-derivation and extension of existing SAR imaging theory for point and diffuse targets. Generic, relatively simple, closed-form expressions for the impulse response, the resolution, and the image bandwidth summarize this unified treatment. Theoretical differences between the imagery of point and diffuse targets are pointed out. Based upon these fundamental differences, a statistical testing procedure is formulated to address the question of scene target density. Background ocean surface wave theory is outlined in preparation for discussions of SAR ocean imaging. Of central importance is the role of the phase velocity, which is the speed of translation of the mean pattern of reflectivity, and the orbital motion, which leads to coherent (phase) modulation, and hence to velocity bunching, acceleration defocus, and target decorrelation. Based upon this theoretical background, one- and two-dimensional simulation models are developed. The one-dimensional simulation addresses the effects of various parameters upon the mean image contrast in a velocity bunching model and guides the development of the two-dimensional simulation. The two-dimensional simulation is unique because each target which constitutes the scene is explicitly considered. This leads to a degree of control and flexibility which is not available from actual SAR imagery. Qualitative and quantitative comparisons are drawn between the simulated and actual SAR imagery to address the key areas of controversy. The assertion that Bragg scattering is a coherent process is defended, despite inability to conclusively verify this using SEASAT data. Comparisons between simulation and C-SAR imagery of waves propagating into ice verify the roles of the scene coherence time and the wave phase velocity. |
Subject |
Radar meteorology Synthetic aperture radar |
Genre |
Thesis/Dissertation |
Type |
Text |
Language | eng |
Collection |
Retrospective Theses and Dissertations, 1919-2007 |
Series | UBC Retrospective Theses Digitization Project |
Date Available | 2010-08-20 |
Provider | Vancouver : University of British Columbia Library |
Rights | For non-commercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https://open.library.ubc.ca/terms_of_use. |
DOI | 10.14288/1.0053106 |
URI | http://hdl.handle.net/2429/27556 |
Degree |
Doctor of Philosophy - PhD |
Program |
Oceanography |
Affiliation |
Science, Faculty of Earth, Ocean and Atmospheric Sciences, Department of |
Degree Grantor | University of British Columbia |
Campus |
UBCV |
Scholarly Level | Graduate |
Aggregated Source Repository | DSpace |
Download
- Media
- UBC_1987_A1 V32.pdf [ 13.69MB ]
- Metadata
- JSON: 1.0053106.json
- JSON-LD: 1.0053106+ld.json
- RDF/XML (Pretty): 1.0053106.xml
- RDF/JSON: 1.0053106+rdf.json
- Turtle: 1.0053106+rdf-turtle.txt
- N-Triples: 1.0053106+rdf-ntriples.txt
- Original Record: 1.0053106 +original-record.json
- Full Text
- 1.0053106.txt
- Citation
- 1.0053106.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Country | Views | Downloads |
---|---|---|
China | 11 | 38 |
United States | 7 | 0 |
France | 3 | 0 |
Japan | 3 | 0 |
India | 2 | 0 |
Switzerland | 2 | 1 |
United Kingdom | 1 | 0 |
Russia | 1 | 0 |
City | Views | Downloads |
---|---|---|
Shenzhen | 5 | 35 |
Ashburn | 5 | 0 |
Guangzhou | 4 | 0 |
Unknown | 4 | 0 |
Tokyo | 3 | 0 |
Zurich | 2 | 1 |
Harbin | 2 | 0 |
Boardman | 1 | 0 |
Tomsk | 1 | 0 |
Mumbai | 1 | 0 |
Ogden | 1 | 0 |
Jammu | 1 | 0 |
{[{ mDataHeader[type] }]} | {[{ month[type] }]} | {[{ tData[type] }]} |
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.831.1-0053106/manifest