Synthetic Aperture Radar Using Non-Uniform Sampling

Total Page:16

File Type:pdf, Size:1020Kb

Synthetic Aperture Radar Using Non-Uniform Sampling -?. qS Synthetic Aperture Radar [Ising Non-Uniform Sampling by JoN¿IHEN ANOREW LECC Thesis submitted for the degree of Doctor of Philosophy Department of Electrical and Electronic Engineering Faculty of Engineering The University of Adelaide Adelaide, South Australia August 1997 Contents Abstract vu Declaration ix Acknowledgments xi List of Figures xrll List of Thbles xYll Glossary xtx Publications xxY I Introduction 1 l.l Motivation 1 I.2 Thesis Outline and Contributions 2 2 Background Information 7 2.I Synthetic Aperture Radar Fundamentals . 7 2.1.1 Range Information . 9 2.1.2 From Real to Synthetic Aperture l0 2.L3 Azimuthal Information 12 2.1.4 PRF Selection 13 2.2 Radar Moving Target Detection 15 2.3 The Effects of Moving Targets on SAR l6 2.4 Radar Waveforms l9 2.5 Existing SAR MTD Techniques 2t 2.6 Proposed SAR MTD Technique 29 ll CONTENTS 2.7 Non-Uniform Sampling 30 2.7.1 Heuristic Justification for Non-Uniform pRI SAR MTD 30 2.7.2 Introduction 31 2.7.3 Minimum SamplingRequiremeuts 55 2J.4 Pseudorandom Sampling JJ 2.7.5 SignalProcessing JJ 2.1.6 Timing Specifications 36 2.1.7 SAR Specifics 39 2.1.8 Non-UniformArrays 4t 2.7.9 Non-Uniform Sampling: Concluding Remarks 4t 2.8 Performance Requirements 42 2.9 Summary 42 3 Moving Thrget Ambiguity Function 43 3.1 Introduction 43 3.2 Ambiguity Function Derivation 44 3.2.1 Radar Scenario and Signals 44 3.2.2 RangeCompression 47 3.2.3 Slow Time Compression . 49 3.2.4 SAR Geometry 51 3.2.5 Stripmapping SAR Ambiguity Function . 52 3.3 Depth of Focus 54 3.4 Moving Target Focussing 55 '3.4.1 Linearity 56 3.4.2 Uniqueldentifiability 56 3.4.3 Squint Mode 58 3.5 Uniform Pzu . 59 3.5.1 Azimuthal Position/Range Velocity . 59 3.5.2 Azimuthal Position/Azimuthal Velocity 62 3.6 Non-UniformPRl 62 3.6.1 The Effect of Two Staggered PRIs . 65 3.7 The Ambiguity Function with Random Sampling 69 3.1.l Range Velocity/Azimuthal Position 1l 3.7.2 Azimuthal Velocity 76 CONTENTS lll 3.8 Conclusions 78 4 Moving Thrget Detection 79 4.1 Introduction 79 4.2 Data Model . 80 4.2.1 Azimuthal Signal 80 4.2.2 Hypotheses 80 4.2.3 Assumptions 8l 4.2.4 DisturbanceModelling 82 4.3 Optimal Detection 85 4.3.1 Numerical Examples . 89 4.3.2 Alternative Optimality Criterion 93 4.3.3 Signal-to-DisturbanceRatio Improvement 93 4.3.4 Receiver Operating Characteristic . 97 4.4 Moving Target Indication 102 4.5 ClutterCancellation 104 4.6 Practical Considerations 105 4.7 Conclusions 106 5 Moving Target Parameter Estimation 107 5.1 Introduction . 107 5.2 Signal Model . 108 5.3 Existing'Work . 111 5.4 Maximum Likelihood Estimation . 113 5.4.1 Bias . 115 5.4.2 Variance . rr7 5.5 Random Sampling . r2t 5.5.1 Example Sampling Distributions . r25 5.6 Simulations . t28 5.7 PracticalConsiderations . t29 5.8 Actual Unknown Parameters . 130 5.8.1 NumericalExamples . t32 5.9 Conclusions . t39 lv CONTENTS 6 Imaging r4t 6.1 Introduction l4t 6.2 Existing Work t42 6.3 Azimuthal Imaging Formulation r43 6.3.1 Focussing t46 6.4 ProcessingTechniques t47 6.4.1 Phase Corrections 148 6.4.2 Weighting 151 6.4.3 Maximum ISLR . t52 6.5 Effects of a Non-Uniform PRI 154 6.5.1 Resolution 156 6.6 The Azimuthal Response with Random Sampling 156 6.6.1 Case 1: Uniformly DistriburedTimes 157 6.6.2 Case 2: Uniformly Distributed Offsets 158 6.7 Conclusion 159 7 Application to Real Data t6l 7.1 Introduction 161 7.2 Moving Target Imaging r62 7.2.1 Non-Uniform PRI Simulation t6l 7.3 Moving Target Parameter Estimation 168 7.4 ClutterCancellation 170 7.5 Conclusions 17l I Summary 173 8.1 Summary of the SAR MTD Technique 173 8.2 Future Research 116 8.3 Conclusion 176 CONTENTS A Extended Ambiguity Function Properties 177 B SAR Clutter Modelling 181 8.1 Time Domain Clutter Modelling 181 8.2 Comparison between Doppler and Time Domain Models 184 C Biased Estimator Lower Variance Bounds 187 D The Complex Gaussian Fisher Information Matrix 189 E Notes Regarding the Optimal Tlansmission Times 193 E.1 Introduction . r93 8.2 Possible Schemes . t94 8.2.1 Deterministic: Heuristic . 194 8.2.2 Deterministic: Optimal . 194 8.2.3 Random Sampling Schemes . 195 E.3 Minimum Redundancy Sampling Scheme . r97 8.4 Conclusions . 198 VI Abstract Synthetic aperture radar is a well established technique for imaging the ground to one side of an airborne platform. Moving target detection is a very useful capability for a long range sensor. The combination of synthetic aperture radar and moving target detection has the potential to produce high resolution ground imagery with superimposed moving target information. Un- fortunately, using conventional imaging data for detecting moving targets leads to ambiguities in the targets' positions and velocities. By utilising a non-uniform pulse repetition interval, the proposed ground imaging / moving target detection radar overcomes this limitation and allows the azimuthal data to be focussed at any velocity of interest, whilst collecting data at the same average rate as a conventional synthetic aperture radar. This approach permits the flexible use of a multimode radar, relaxes the specifications of data acquisition systems, affords a degree of protection against electronic countermeasures and retains alarge unambiguous range swath, but with the added complexity of processing the non-uniform samples. This thesis investigates the technique in detail, incorporating optimal target detection strategies and azimuthal beampatterns, and demonstrates it using simulations of non-uniform transmissions with real data. vil vlll X Acknowledgments I wish to thank a number of people and organisations, without whom this thesis would not exist. Dr Alan Bolton is a colleague and friend who provided the original inspiration for this work and gave continual support and encouragement throughout its duration. Prof. Douglas Gray, my Ph.D. supervisor, changed my attitude from one of 'try it and see' to one of careful and thorough analysis. His expertise in sonar beamforming provided me with useful insights regarding the treatment of synthetic aperture radar as a phased array. Dr Don Sinnott, as Chief of Microwave Radar Division, approved my sponsorship through the Defence Science and Technology Organisation, Salisbury. My MRD supervisor, Mar- ian Viola, kindly decided he could do without me during the period of my candidature and eagerly awaited its completion. The Cooperative Research Centre for Sensor Signal and Information Processing provided excellent facilities and great people to work with. I will never forget the superlative morning teas. My family, and particularly my darling Jan, provided encouragement regarding the further- ing of my career. Jan was quite happy for me to work hard during the weekends! XI xll List of Figures 2.1 Simplified stripmapping SAR geometry 8 2.2 An overview of SAR processing. 8 2.3 The first order effects of moving targets on a SAR image 18 2.4 A moving target's return after being processed as if the target were stationary 19 2.5 The distances to two targets seen by the radar 30 2.6 Radar azimuthal phase returns. 3l 2.7 Radar azimuthal phase returns; non-uniform PR[. 32 2.8 Correlation filtering of non-uniformly sampled data. 35 2.9 Estimated spectrum; uniform sampling 37 2.ro Estimated spectrum; random offsets. 38 2.tt Non-uniform sampling notation 38 2.t2 Estimated spectrum; random intervals. 39 3.1 A block diagram of the radar model M 3.2 The SAR antenna pattern 45 3.3 SAR imaging scenario 49 3.4 The range velocitylazimuthal position ambiguity function . 60 3.5 The peak of the range velocitylazimuthal position AF . 60 3.6 The azimuthal velocity/azimuthal position AF . 62 3.1 The peak azimuthal velocitylazimuthal position AF . 63 3.8 The range velocitylazimuthal position AF; non-uniform PRI 63 3.9 Slices of the non-uniform PRI range velocitylazimuthal position AF 64 3.10 The non-uniform PRI azimuthal velocity/azimuthal position AF 64 3.1 l The peak non-uniform PRI azimuthal velocitylazimuthal position AF . ' 65 3.12 The normalised SAR azimuthal AF with one PRI plotted against range velocity. 67 3.r3 The normalised AF with superimposed sample times. 68 xiii XIV LIST OFFIGURES 3.t4 The normalised AF with two PRIs: the first is 0.3 x the nominal value. 68 3. r5 The normalised AF with two alternating PRIs: the first is 0.5 x the nominal value. 69 3.16 The expected range velocity/azimuthal position ambiguity function for uni- formly distributed sampling l2 3.t] The peak of the expected range velocitylazimuthal position ambiguity function for uniformly distributed sampling t3 3.18 The expected range velocity/azimuthal position ambiguity function for uni- formly distributed offsets 75 3.19 The normalised AF with randomly chosen PRIs . 16 4.1 Moving target detection spectra; high PRF 9l 4.2 MTD spectra; high PRF, non-uniform PRI 9l 4.3 MTD spectra;low PRF 92 4.4 MTD spectra; low PRF, non-uniform PRI 92 4.5 Signal to disturbance ratio improvement; high PRF 95 4.6 SDR improvement; high PRF, non-uniform PRI 96 4.7 SDR improvement; low PRF 96 4.8 SDR improvement; low PRF, non-uniform PRI 97 4.9 Receiver operating characteristics for an optimal filter, high PRF r0l 4.to ROC for an optimal filter,low PRF 101 4.tt ROC for a matched filter, low PRF .
Recommended publications
  • Ultrafast, Broadband and Multi-Pulse Transmissions for Ultrasonic Imaging Yanis Mehdi Benane
    Ultrafast, broadband and multi-pulse transmissions for ultrasonic imaging Yanis Mehdi Benane To cite this version: Yanis Mehdi Benane. Ultrafast, broadband and multi-pulse transmissions for ultrasonic imaging. Signal and Image processing. Université Lyon 1 - Claude Bernard, 2018. English. tel-02072384v1 HAL Id: tel-02072384 https://hal.archives-ouvertes.fr/tel-02072384v1 Submitted on 19 Mar 2019 (v1), last revised 11 Jul 2019 (v2) HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. N°d’ordre NNT : xxx THESE de DOCTORAT DE L’UNIVERSITE DE LYON opérée au sein de l’Université Claude Bernard Lyon 1 Ecole Doctorale N° 162 MECANIQUE, ENERGETIQUE, GENIE CIVIL, ACOUSTIQUE Spécialité de doctorat : ACOUSTIQUE ET TRAITEMENT DE SIGNAL Discipline : Imagerie ultrasonore Soutenue publiquement le 10/12/2018, par : Mehdi Yanis BENANE Ultrafast, broadband and multi-pulse transmissions for ultrasonic imaging Devant le jury composé de : Président du jury : …………… Rapporteur : Kouame, Denis Professeur des Universités, Université de Toulouse Petitjean, Caroline Maître de conférence,
    [Show full text]
  • Optimization of Short Probe Linear Frequency Modulated Signal Parameters
    ARCHIVES OF ACOUSTICS DOI: 10.2478/v10168-011-0057-1 36, 4, 861–871 (2011) Optimization of Short Probe Linear Frequency Modulated Signal Parameters Włodzimierz POGRIBNY(1),(2), Tadeusz LESZCZYŃSKI(1) (1)University of Technology and Life Sciences Kaliskiego 7, 85-796 Bydgoszcz, Poland e-mail: pohry, leszcz @utp.edu.pl { } (2)University of Economy in Bydgoszcz Garbary 2, 85-229 Bydgoszcz, Poland (received June 17, 2011; accepted October 20, 2011 ) In many physical experiments, linear frequency modulated (LFM) signals are widely used to probe objects in different environments, from outer-space to under- water. These signals allow a significant improvement in measurement resolution, even when the observation distance is great. For example, using LFM probe sig- nals in underwater investigations enables discovery of even small objects covered by bottom sediments. Recognition of LFM (chirp) signals depends on their compression based on matched filtering. This work presents two simple solutions to improve the reso- lution of the short chirp signals recognition. These methods are effective only if synchronization between the signal and matched filter (MF) is obtained. This work describes both the aforementioned methods and a method of minimizing the effects of the lack of synchronization. The proposed matched filtering method, with the use of n parallel MFs and other techniques, allows only one sample to be obtained in the main lobe and to accurately locate its position in the appropriate sampling period Ts with accuracy Ts/n. These approaches are appropriate for use in probe signal processing. Keywords: sonar, LFM signal, matched filtering. 1. Introduction For many physical investigations, short chirp signals are especially expedi- ent because they enable the probe rate to be increased (Grelowska, Koza- czka, 2010; Mahafza, Elsherbeni, 2004; Pogribny et al., 2002; Schock, LeBlanc, 1990).
    [Show full text]
  • Digital Beam-Forming for High Resolution Wide Swath Real and Synthetic Aperture Radar
    Forschungsberichte aus IHEIHE dem Institut für Höchstfrequenztechnik und Elektronik der Universität Karlsruhe (TH) Herausgeber: Prof. Dr.-Ing. W. Wiesbeck Marwan Younis Digital Beam-Forming for High Resolution Wide Swath Real and Synthetic Aperture Radar Band 42 Copyright: Institut für Höchstfrequenztechnik und Elektronik Universität Karlsruhe (TH) alle Rechte vorbehalten Druck: Druckerei Gunter Dünnbier, 02779 Großschönau, Tel. 035841-36757 ISSN: 0942-2935 Forschungsberichte aus dem Institut für Höchstfrequenztechnik und Elektronik der Universität Karlsruhe (TH) Vorwort des Herausgebers Als Anfang der 60-iger Jahre die vorteilhafte Nutzung der digitalen Technik gegenüber der analogen Technik in vielen Bereichen erkennbar wurde, waren die meisten der analogen Technik nahe stehenden Ingenieuren und Wissenschaftler verunsichert. Es wurden Fronten aufgebaut und in zahlreichen Publikationen nachgewiesen, warum die digitale Technik nicht leistungsfähig ist. Diese Einstellungen gehören zwischenzeitlich der Vergangenheit an. Die nachkommenden Generationen an Ingenieuren und Wissenschaftlern standen digitalen Problemlösungen offener gegenüber. Heute hat die Digitaltechnik in vielen Bereichen vorteilhaft die analoge Technik abgelöst. Um dies zu untermauern, muss man sich nicht auf die Errungenschaften in der Computertechnik beschränken. Herausragende Beispiele sind auch der digitale Mobilfunk, die digitale Fotografie und der digitale Rundfunk. Wie sehr gerade beim Letzteren die analoge Technik verwurzelt war und ist, zeigt sich daran, dass das Wort digital in den Bezeichnungen DAB, DVB-T explizit verwendet wird. So sind zwischenzeitlich die meisten analogen Basisstationen gefallen. Hartnäckig hält sich allerdings die analoge Technik in der Formung der Richtcharakteristiken von Antennen. Die Dissertation von Dr. Younis liefert in diesem Bereich grundlegende Beiträge zum sog. digital beam-forming speziell für die Radarsensorik. Hierfür wird in Mehrantennensystemen durch geeignete Raum-Zeit Signalverarbeitung eine äußerst flexible Strahlformung erreicht.
    [Show full text]
  • Asynchronous Chirp Slope Keying for Underwater Acoustic Communication
    Preprints (www.preprints.org) | NOT PEER-REVIEWED | Posted: 30 March 2021 doi:10.20944/preprints202103.0750.v1 Article Asynchronous Chirp Slope Keying for Underwater Acoustic Communication Dominik Jan Schott 1,†∗, Andrea Gabbrielli 2, Wenxin Xiong 2, Georg Fischer 3, Fabian Höflinger 1,3, Johannes Wendeberg2, Christian Schindelhauer 2, Stefan Johann Rupitsch 1 1 Dep. of Microsystems Engin. (IMTEK), University of Freiburg, Germany 2 Dep. of Computer Sci. (IIF), University of Freiburg, Germany 3 Fraunhofer EMI, Efringen-Kirchen, Germany * Correspondence: [email protected] † Current address: Georges-Koehler-Allee 106, 79110 Freiburg, Germany 1 Abstract: We propose an asynchronous acoustic chirp slope keying to map bit sequences on single 2 or multiple bands without preamble or error correction coding on the physical layer. Details of the 3 implementation are disclosed and discussed, the performance verified on laboratory scale in a 4 pool measurement, as well as simulated for a channel containing Rayleigh fading and Additive 5 White Gaussian Noise. For time-bandwidth products of 50 in single band mode, a raw data rate of 6 100 bit/s is simulated to achieve bit error rates below 0.001 for signal-to-noise ratios above -6 dB. 7 In dual-band mode and a data rate of 200 bit/s, this bit error level was achieved for signal-to-noise 8 ratios above 0 dB for time-bandwidth product of 25. The packet error rates follow this behavior 9 with an offset of 1 dB. 10 Keywords: Underwater Communication; Wireless Communication; Acoustic Communication; 11 Ultrasound Acoustics; Digital Signal Processing; Chirp Modulation; Chirp Slope Keying; Chirp 12 Spread Spectrum; Citation: Schott, D.
    [Show full text]
  • Asynchronous Chirp Slope Keying for Underwater Acoustic Communication
    sensors Article Asynchronous Chirp Slope Keying for Underwater Acoustic Communication Dominik Jan Schott 1,* , Andrea Gabbrielli 1 , Wenxin Xiong 2 , Georg Fischer 3 , Fabian Höflinger 1,3 , Johannes Wendeberg 2, Christian Schindelhauer 2 and Stefan Johann Rupitsch 1 1 Department of Microsystems Engineering (IMTEK), University of Freiburg, 79110 Freiburg, Germany; [email protected] (A.G.); fabian.hoefl[email protected] (F.H.); [email protected] (S.J.R.) 2 Department of Computer Science (IIF), University of Freiburg, 79110 Freiburg, Germany; [email protected] (W.X.); [email protected] (J.W.); [email protected] (C.S.) 3 Fraunhofer EMI, 79588 Efringen-Kirchen, Germany; [email protected] * Correspondence: [email protected] Abstract: We propose an asynchronous acoustic chirp slope keying to map short bit sequences on single or multiple bands without preamble or error correction coding on the physical layer. We introduce a symbol detection scheme in the demodulator that uses the superposed matched filter results of up and down chirp references to estimate the symbol timing, which removes the requirement of a preamble for symbol synchronization. Details of the implementation are disclosed and discussed, and the performance is verified in a pool measurement on laboratory scale, as well as the simulation for a channel containing Rayleigh fading and Additive White Gaussian Noise. For time-bandwidth products (TB) of 50 in single band mode, a raw data rate of 100 bit/s is simulated − Citation: Schott, D.J.; Gabbrielli, A.; to achieve bit error rates (BER) below 0.001 for signal-to-noise ratios above 6 dB.
    [Show full text]
  • A Multi-Path Compensation Method for Ranging in Wearable Ultrasonic Sensor Networks for Human Gait Analysis
    sensors Article A Multi-Path Compensation Method for Ranging in Wearable Ultrasonic Sensor Networks for Human Gait Analysis Karalikkadan Ashhar 1,* , Mohammad Omar Khyam 2 and Cheong Boon Soh 1 1 School of Electrical and Electronic Engineering, Nanyang Technological University, 50 Nanyang Ave, Singapore 639798, Singapore; [email protected] 2 Department of Mechanical Engineering, Virginia Tech, Blacksburg, VA 24060, USA; [email protected] * Correspondence: [email protected] Received: 15 February 2019; Accepted: 14 March 2019; Published: 18 March 2019 Abstract: Gait analysis in unrestrained environments can be done with a single wearable ultrasonic sensor node on the lower limb and four fixed anchor nodes. The accuracy demanded by such systems is very high. Chirp signals can provide better ranging and localization performance in ultrasonic systems. However, we cannot neglect the multi-path effect in typical indoor environments for ultrasonic signals. The multi-path components closer to the line of sight component cannot be identified during correlation reception which leads to errors in the estimated range and which in turn affects the localization and tracking performance. We propose a novel method to reduce the multi-path effect in ultrasonic sensor networks in typical indoor environments. A gait analysis system with one mobile node attached to the lower limb was designed to test the performance of the proposed system during an indoor treadmill walking experiment. An optical motion capture system was used as a benchmark for the experiments. The proposed method gave better tracking accuracy compared to conventional coherent receivers. The static measurements gave 2.45 mm standard deviation compared to 10.45 mm using the classical approach.
    [Show full text]
  • PERFORMANCE EVALUATION of CHIRP SPREAD SPECTRUM SYSTEM and LAND MOBILE SATELLITE SYSTEM by COMPUTER SIMULATION Bv Junghwan Kim
    PERFORMANCE EVALUATION OF CHIRP SPREAD SPECTRUM SYSTEM AND LAND MOBILE SATELLITE SYSTEM BY COMPUTER SIMULATION bv Junghwan Kim _ l Dissertation submitted to the Faculty of the Virginia Polytechnlc Institute and State University in partial fullillment of the requirements for the degree of Doctor of Philosophy in ° Electrical Engineering APPROVED: Timothy Prstt, Chairman - Richard O. Claus ‘ July, 1988 Blacksburg, Virginia PERFORMANCE EVALUATION OF CHIRP SPREAD SPECTRUM SYSTEM AND V LAND MOBILE SATELLITE SYSTEM BY COMPUTER SIMULATION bv Junghwan Kim Timothy Pratt, Chairman Electrical Engineering (ABSTRACT) 1— F The work presented in this dissertation examines the performance of two satellite radio j communication systems by computer simulation. Two simulations were separately performed for a spread spectrum chirp system as an analog communications system, and for the Land Mobile Satellite System (LMSS) channel as a digital communications system. For the simulation of analog communications, a spread spectrum system using chirp techniques called ’Coded Multiple Chirp Spread Spectrum’ was proposed as a simple, cost-effective alternative for conventional spread spectrum systems. Its application as a spread spectrum overlay service on analog FM-TV was examined through the mutual interference analysis and spectral analysis using soltware programming. For the simulation of digital communications, various digital modulation schemes as well as channel encoding, block interleaving/deinterleaving, and dltferential encoding techniques were used for a thorough performance evaluation of a Land Mobile Satellite System under fading conditions. For this purpose, an LMSS fading channel simulator capable of simulatlng diverse fading characteristics for a satellite channel was designed and tested to yield various performance measures such as symbol error rate and average bit error rate.
    [Show full text]
  • Effect of Photoacoustic Radar Chirp Parameters on Profilometric Information
    Effect of Photoacoustic Radar Chirp Parameters on Profilometric Information A thesis submitted to the Faculty of Engineering in partial fulfillment of the requirements for the degree of Master of Applied Science in Mechanical Engineering by Zuwen Sun Ottawa-Carleton Institute for Mechanical and Aerospace Engineering University of Ottawa Ottawa, Ontario, Canada, K1N 6N5 January 2018 © Zuwen Sun, Ottawa, Canada, 2018 Abstract Photoacoustic imaging for biomedical application has attracted much research in recent years. To date, most of the work has focused on pulsed photoacoustics. Recent developments have seen the implementation of a radar pulse compression methodology into continuous wave photoacoustic modality, however very little theory has been developed in support of this approach. In this thesis, the one-dimensional theory of radar photoacoustics for pulse compressed linear frequency modulated continuous sinusoidal laser photoacoustics is developed. The effect of the chirp parameters on the corresponding photoacoustic signal is investigated, and guidelines for choosing the chirp parameters for absorber profilometric detection are given based on the developed theory and simulations. Simulated results are also compared to available experimental results and show a good agreement. Key words: Photoacoustic, Pulse compression, Autocorrelation, photoacoustic radar. ii Acknowledgements My deep appreciation is presented to Dr. Natalie Baddour, who is my supervisor, for her patient help and professional advice during my study, as well as her financial support which made my life in Canada easier. I would also like to appreciate Dr. Liang Ming, who first gave me the opportunity to study at the University of Ottawa. At last, my grateful thanks to my parents, parents in law, and my uncle and aunt for their endless help and support.
    [Show full text]