RADAR CROSS SECTION ENHANCEMENTFOR RADAR NAVIGATION AND REMOTE SENSINGByDavid George MichelsonB. A. Sc., The University of British Columbia, 1982M. A. Sc., The University of British Columbia, 1986A THESIS SUBMITTED IN PARTIAL FULFILLMENT OFTHE REQUIREMENTS FOR THE DEGREE OFDOCTOR OF PHILOSOPHYinTHE FACULTY OF GRADUATE STUDIESDEPARTMENT OF ELECTRICAL ENGINEERINGWe accept this thesis as conformingto the required standardTHE UNIVERSITY OF BRITISH COLUMBIADecember 1993®David George Michelson, 1993In presenting this thesis in partial fulfilment of the requirements for an advanced degree at theUniversity of British Columbia, I agree that the Library shall make it freely available for reference and study. I further agree that permission for extensive copying of this thesis for scholarlypurposes may be granted by the head of my department or by his or her representatives. Itis understood that copying or publication of this thesis for financial gain shall not be allowedwithout my written permission.David G. MichelsonThe University of British Columbia,Department of Electrical Engineering,2356 Main Mall,Vancouver, B.C.CanadaV6T 1Z4Date: 29 December 1993AbstractResearch Supervisor: Prof. E.V. JullRecent developments in radar navigation and remote sensing have led to a requirement forrugged yet inexpensive location markers and calibration targets which present both a very largescattering cross section and a specified polarization response over a wide angular range. Thisstudy considers several problems related to the analysis and design of passive radar targetsderived from corner reflectors.Transformation of the polarization response of a target between global and local coordinateframes is shown to correspond to rotation of the polarization basis by a prescribed angle which isa function of both the coordinate transformation matrix and the direction of propagation. Oncethe angle of rotation has been determined using either spherical trigonometry or vector algebra,any polarization descriptor can be transformed between coordinate frames by application of asuitable rotation operator.The scattering cross section and angular coverage of a conventional trihedral corner reflectorcan be altered by modifying the size and shape of its reflecting panels. A numerical algorithmbased on physical optics is used to predict the contribution of triple-bounce reflections to theresponse of a reflector with polygonal panels of arbitrary shape. If three-fold symmetry isbroken and the reflector is simply required to present bilateral symmetry, it is found that thescattering cross section, elevation beamwidth, and azimuthal beamwidth of the reflector can bechosen independently of each other.A method for altering the polarization response of a conventional trihedral corner reflectorby adding conducting fins or corrugations to one its interior surfaces is proposed. In calculatingdesign curves for twist-polarizing or circularly polarizing reflectors by mode-matching, optimum accuracy and efficiency are obtained by setting the ratio of free space to groove modesequal to the ratio of groove width to the period. Methods for obtaining linear and circularIIpolarization selective responses are considered. The contribution oftriple-bouhce reflections tothe response of such reflectors is a function of the direction ofincidence, the orientation of thereflector, the dimensions of the corrugations, and the size andshape of the reflecting panels.Experimental results show that prototype twist-polarizing and circularly polarizingreflectorsrespond essentially as predicted.InTable of ContentsAbstractjjTable of Contents ivList of TablesviiList of Figures viiiAcknowledgements xviii1 Introduction1.1 Background and Motivation11.2 Outline 7References 92 Transformation of Polarization Descriptors Between CoordinateFrames2.1 Introduction 112.2 Wave Polarization132.3 Transformation of Coordinates162.4 Evaluation of the Coordinate Transformation Matrix232.5 Rotation of the Basis of Common Polarization Descriptors262.6 Conclusions44References453 Truncation and Compensation of TrihedralCorner Reflectors3.1 Introduction473.2 Analysis49iv3.3 Reflectors with Three-Fold Symmetry3.4 Reflectors with Bilateral Symmetry3.5 Effect of Errors in Construction on Reflector Performance3.6 ConclusionsReferences4 Depolarizing Trihedral Corner Reflectors4.1 Introduction4.2 Scattering by a Conducting Grating with Rectangular Grooves4.3 Scattering by a Depolarizing Trihedral Corner Reflector4.4 Numerical and Experimental Results4.4.1 Design and Construction of the Prototype Trihedral Corner4.4.2 Polarization Response4.4.3 Azimuthal Response4.5 ConclusionsReferences5 Summary, Conclusions, and Recommendations5.1 Summary and Conclusions5.2 Recommendations for Further WorkAppendicesA Design Curves for Top Hat ReflectorsA.1 IntroductionA.2 AnalysisA.3 Design CurvesA.4 Design ExampleReferences63718485879092104112Reflectors . . 112117120126128131136140141146152156V• 157• 159• 159• 165169174187C Circular Polarization Selective ReflectorsC.1 IntroductionC.2 ConceptC.3 Proposed ImplementationsC.4 DiscussionReferencesD Experimental ArrangementD.1 IntroductionD.2 OverviewD.3 CW Radar ApparatusD.4 Digital Pattern Recorder.D.5 Facility EvaluationReferences189189192198199B Scattering by a Conducting Grating withRectangular GroovesB.1 IntroductionB.2 AnalysisB.2.1 TM PolarizationB.2.2 TE PolarizationB.3 Verification of Numerical ResultsB.4 ImplementationReferences200200202208214221viList of Tables3.1 Response Characteristics of Selected Trihedral Corner Reflectorswith Three-foldSymmetry 703.2 Response Characteristics of Selected Trihedral CornerReflectors with BilateralSymmetry 834.1 Tolerances on the Corner Angles of a Prototype Trihedral Corner ReflectorwithTriangular Panels1164.2 Dimensions of the Prototype Reflection Polarizers116A.1 Response Characteristics of Selected Top Hat Reflectorsatf= 10 GHz 154B.1 Definition of Symbols158D.1 Digital Pattern Recorder Program Modules209D.2 Response of Trihedral Corner Reflectors with TriangularPanels at 9.445 GHz . . 216ViiList of Figures1.1 Detection of a point targetin ground clutter 31.2 Probability distribution functions ofclutter and a signal embedded in clutter. . 31.3 NASA/JPL synthetic aperture radarcalibration site at Goldstone, California.. 51.4 Radar-assisted positioningwith respect to cooperative shore-based targets. . . 51.5 Use of range and azimuth gatesto isolate shore station reflectors from surrounding clutter61.6 Relative size ofcorner reflectors which present the same maximum radarcrosssection (4500 m2) atf= 10 GHz62.1 Coordinate frames used to definepolarization state in radar scattering problems 112.2 Coordinate system for a plane wavepropagating in the direction A 142.3 A polarization ellipse with semi-majoraxis OA, semi-minor axis OB, and tiltangle r152.4 Mapping of polarization states onto aPoincar sphere 152.5 A polarization ellipse showingthe relationship between the tilt angles Tand Tin the xyz and x’y’z’ coordinate frames172.6 A perspective projectionof the parallels of two coordinate frameswhich arerelated by pure rotation192.7 A Mercator projection ofthe parallels of two coordinate frames whichare relatedby pure rotation192.8 Poincaré sphere representationof the polarization state W of aplane wave. . 272.9 Definition of the phasereference for orthogonal circular components28viii2.10 Coordinate systems and scattering geometry for the forward scatteringalignment(FSA) convention412.11 Coordinate systems and scattering geometry for the backscatter alignment (BSA)convention413.1 Relative sizes of trihedral corner reflectors used as location markers and calibration targets in radar navigation, radar-assisted positioning, and radar remotesensing 473.2 Scattering by a trihedral corner reflector with triangular panels 503.3 Problem geometry and coordinate system forscattering by a trihedral cornerreflector composed of triangular, elliptical, or rectangular panels with corners ofarbitrary length 503.4 Spencer’s model for the equivalent flat plate area of a trihedra)corner reflector 513.5 An implementation of Robertson’s model for thecontribution of triple-bouncereflections to the response of a trihedral corner reflector 543.6 Examples of trihedral corner reflectors for which Spencer’smodel gives accuratepredictions of the equivalent flat plate area 553.7 Examples of trihedral corner reflectors for which Spencer’smodel gives erroneouspredictions of the equivalent flat plate area 563.8 Transformation of coordinates by pure rotation583.9 Execution of the Weiler-Atherton polygon clipping algorithm623.10 Reflector coordinate frame and global coordinate frame 633.11 Angular coverage of a trihedral corner reflector with triangular panels 663.12 Angular coverage of a trihedral corner reflector with circular panels 663.13 Angular coverage of a trihedral corner reflector with square panels 673.14 Boresight view and panel geometry of trihedral corner reflectors with Robertsonpanels 67ix3.15 Angular coverage of atrihedral corner reflector with Robertsonpanels whered/d = 0.25683.16 Angular coverage ofa trihedral corner reflectorwith Robertson panels whered/d = 0.50683.17 Angular coverage of atrihedral corner reflector with Robertsonpanels whered/d = 0.75693.18 Cumulative probabilitydistribution of the response oftrihedral corner reflectorswith triangular, circular, square,and Robertson panels(where d/d = 0.25)over the quadrant defined by theaxes of the trihedral693.19 Angle of maximumresponse6mof bilaterally symmetrictrihedral corner reflectors with triangular, elliptical,and rectangular panels vs. the reflectoraspectratio c/a743.20 Angle of maximumresponse8mand angle of the normal to theaperture Oof bilaterally symmetric trihedralcorner reflectors with triangular panelsvs. thereflector aspect ratio c/a743.21 Maximum response of a bilaterallysymmetric trihedral corner reflectorswithtriangular, elliptical, and rectangularpanels vs. the reflector aspect ratioc/a. . 753.22 Azimuthal and elevationbeamwidths of a bilaterally symmetrictrihedral cornerreflector with triangularpanels vs. the reflector aspect ratio c/a753.23 Azimuthal andelevation beamwidths of a bilaterallysymmetric trihedral cornerreflector with elliptical panels vs.the reflector aspect ratio c/a763.24 Azimuthal and elevationbeamwidths of a bilaterallysymmetric trihedral cornerreflector with rectangular panels vs.the reflector aspect ratio c/a763.25 Angular coverage of abilaterally symmetric trihedral cornerreflector with triangular panels and reflector aspectratio c/a = 0.25773.26 Angular coverage of a bilaterallysymmetric trihedral corner reflectorwith triangular panels and reflector aspect ratioc/a = 4.077x3.27 Angular coverage of a bilaterally symmetric trihedralcorner reflector with triangular side panels and a circular center panel 783.28 Angular coverage of a bilaterally symmetric trihedral corner reflector withtriangular side panels and a square center panel783.29 Angular coverage of a bilaterally symmetric trihedralcorner reflector with circular side panels and a triangular center panel 793.30 Angular coverage of a bilaterally symmetric trihedral corner reflectorwith circular side panels and a square center panel 793.31 Angular coverage of a bilaterally symmetric trihedral corner reflectorwith squareside panels and a triangular center panel 803.32 Angular coverage of a bilaterally symmetric trihedral corner reflector with squareside panels and a circular center panel 803.33 Evolution of Lanziner’s bilaterally symmetrictrihedral corner reflector 813.34 Angular coverage of a bilaterally symmetric trihedral corner reflectorwith truncated triangular side panels and a triangular center panel813.35 Angular coverage of a bilaterally symmetric trihedralcorner reflector with truncated and compensated triangular side panels and a triangular centerpanel. . . 823.36 Angular coverage of a bilaterally symmetric trihedralcorner reflector with truncated, compensated, and extended side panels and a triangular center panel.. . 823.37 Effect of errors in all three corner angles on theresponse of a trihedral cornerreflector with triangular panels for incidence along the the symmetry axis844.1 Methods for altering the polarization response of aconventional trihedral cornerreflector904.2 Problem geometry for scattering by a conductinggrating with rectangular grooves 934.3 Visible diffracted orders as a function of the gratingperiod (in wavelengths) andthe angle of incidence where the plane of incidence is normal to the grating axis93xi4.4 Normalized response of a regular reflector as a functionof the polarization stateof the incident wave964.5 Normalized response of a twist-polarizing reflector asa function of the polarization state of the incident wave964.6 Normalized response of a circularly polarizing reflector asa function of the polarization state of the incident wave974.7 Normalized response of a vertical polarization selective reflectoras a function ofthe polarization state of the incident wave 974.8 Twist polarizer design curves for normal incidence1004.9 Twist polarizer design curves for 45 degree incidence1004.10 Circular polarizer design curves for normal incidence1014.11 Circular polarizer design curves for 45 degree incidence1014.12 Linear polarization selective reflectors derived fromcorrugated surfaces 1024.13 Replacement of one panel of a trihedral cornerreflector by a reflection polarizerderived from a corrugated surface1054.14 Twist polarizer design curves for incidence along thesymmetry axis of a trihedralcorner reflector1064.15 Circular polarizer design curves forincidence along the symmetry axis of a tnhedral corner reflector1064.16 Reflector coordinate frame and global coordinateframe 1084.17 Angle of rotation of the projection of the gratingaxis onto the view plane forincidence along the symmetry axis and 30 degrees off thesymmetry axis 1084.18 Profile view of the RCS measurement range1134.19 Measurement of the response of a prototype trihedral corner reflector1134.20 Construction details of the prototype depolarizingtrihedral corner reflector. . . 1154.21 Photograph of the prototype twist-polarizing trihedral corner reflector mountedon the antenna range model tower115xii4.22 A single segmentof a prototype reflection polarizer1164.23 Evaluationof the polarization response ofa radar target1184.24 Polarization responseof a prototype regulartrihedral corner reflector as a function of rotation about the boresight1184.25 Polarization responseof a prototype twist-polarizingtrihedral corner reflectorasa function of rotationabout the boresight1194.26 Polarizationresponse of a prototype circularly-polarizingtrihedral corner reflector as a function of rotation aboutthe boresight1194.27 Co-polar azimuthalresponse pattern of a prototyperegular trihedral corner reflector for rotation anglea= 00and vertically polarized transmission1214.28 Cross-polar azimuthalresponse patterns of a prototyperegular trihedral cornerreflector for rotation anglea00and vertically polarized transmission1214.29 Co-polar azimuthalresponse pattern of a prototypetwist-polarizing trihedralcorner reflector for rotationangle a= 00and vertically polarizedtransmission. . 1224.30 Cross-polarazimuthal response patterns of aprototype twist-polarizingtrihedralcorner reflector for rotation anglea= 00and verticaily polarized transmission.. 1224.31 Co-polar azimuthalresponse pattern of a prototypetwist-polarizing trihedralcorner reflector for rotationangle a = 45° and verticallypolarized transmission. 1234.32 Cross-polar azimuthalresponse patterns of a prototypetwist-polarizing trihedralcorner reflector for rotationangle a = 45° and verticallypolarized transmission. 1234.33 Co-polar azimuthalresponse pattern of a prototype circularly-polarizingtrihedralcorner reflector for rotationangle a = 0° and verticallypolarized transmission. . 1244.34 Cross-polar azimuthalresponse patterns of a prototypecircularly-polarizing tnhedral corner reflector forrotation angle a =00and vertically polarized transmission1244.35 Co-polar azimuthalresponse pattern of a prototypecircularly-polarizing trihedralcorner reflector for rotation angle a =45° and vertically polarizedtransmission. 125XIII4.36 Cross-polar azimuthalresponse patterns of a prototype circularly-polarizingtnhedral corner reflector forrotation angle o — 45° and vertically polarizedtransmission125A.1 Problem geometry forscattering by a top hat reflector141A.2 A simplified ray-opticalmodel for scattering by a top hat reflector142A.3 Problem geometry forforward scattering by a cylinder143A.4 Angle of maximum responseof a top hat reflector vs. , the ratioof the annuluswidth to the cylinder height146A.5 Maximum scatteringcross section of a top hat reflectorvs. , the ratio of theannulus width to the cylinder height forfixed values of a and c148A.6 Scale factors for the cylinderradius a and height c of a top hat reflectorvs. theratio of the annulus width to thecylinder height148A.7 Half-power elevationbeamwidth of a top hat reflector vs.the ratio of the annuluswidth to the cylinder height150A.8 1 dB elevation beamwidth of atop hat reflector vs. the ratio of theannulus widthto the cylinder height151A.9 Angle of maximumresponse and angles of median responsefor 1 and 3 dBelevation beamwidths vs. the ratio of theannulus width to the cylinderheight. . 152A.10 Relative size of selectedtop hat reflectors which present thesame maximumscattering cross section154A.11 Elevation response patternsof selected top hat reflectors whichpresent the samemaximum scattering cross section155B.1 A unit cell of a conducting gratingwith rectangular grooves158xivB.2 Convergence of the phasedifference between the TE and TMspecular reflectioncoefficients with the number of groove modes for areflection twist polarizer withd = 0.3333A, aid = 0.5000, and h= 0.2302A illuminated by a planewave incidentat qS 45 degrees172B.3 Convergence of the phase differencebetween the TE and TM specular reflectioncoefficients with the number of groove modes fora reflection twist polarizer withd = 0.3333), aid = 0.9999, and h = 0.3172A illuminatedby a plane wave incidentat = 45 degrees172B.4 Convergence of the phase differencebetween the TE and TM specular reflectioncoefficients with the number of groove modes for areflection circular polarizerwith d = 0.3333A, a/d = 0.5000, and h = 0.1466Ailluminated by a plane waveincident at = 45 degrees173B.5 Convergence of the phase differencebetween the TE and TM specular reflectioncoefficients with the number of groovemodes for a reflection circular polarizerwith d = 0.3333), a/d = 0.9999, andh = 0.1641A illuminated by planewaveincident at = 45 degrees173B.6 Hierarchy of subprogramscalled by subroutines TMREFL and TEREFL175B.7 Combined execution time ofsubroutines TMREFL and TEREFL on a Sun 41380workstation vs. the number of groovemodes used in the solution 176B.8 Combined execution time ofsubroutines TMREFL and TEREFL on a Sun 41380workstation vs. the number of free spacemodes used in the solution 176B.9 Execution profile of subroutines TMREFLand TEREFL 177C.1 Normalized response of a left circularpolarization selective reflector as a functionof the polarization state of the incidentwave 191C.2 Normalized response of a right circularpolarization selective reflector as a function of the polarization state of the incidentwave191xvC.3 A proposed implementationof a circular polarization selectivereflector using atransmission circular polarizer and alinear polarization selective reflector.C.4 A parallel plate transmissioncircular polarizerC.5 A proposedimplementation of a circular polarizationselective reflector using acircular polarization selective surface and atrihedral twist reflectorC.6 Scattering by aright circular polarization selectivetrihedral corner reflector.C.7 A linear polarizationselective gridded trihedral and itsco-polar and cross-polarazimuthal response patternsC.8 A bilaterally symmetrictrihedral corner reflector with triangular panelsC.9 Elements of a circularpolarization selective surfaceD.1 Block diagram of theradar cross section measurement rangeD.2 Photograph of the radarcross section measurement rangeD.3 Block diagram of the CWradar transmitterD.4 Photograph of the CW radartransmitterD.5 Mounting arrangement forthe transmitting and receiving antennasD.6 Proffle view of the RCSmeasurement rangeD.7 Arrangement for mechanicallyaligning the transmitting and receivingantennas.D.8 Radar cross section measurementrange link budgetD.9 Photograph of thedigital pattern recorder, positionercontrol unit, and portablemicrowave receiverD.10 Digital pattern recorder: synchrotest screenD.11 Digital pattern recorder:receiver calibration screenD.12 Equipment configuration forperforming relative and absolutecalibration of theCW radar apparatusD.13 Digital pattern recorder: parameterentry screenD.14 Digital pattern recorder: data acquisitionscreenD.15 Error model for radar cross section measurement• 193• 194• 195195196• 197198201202• 204204205• 206206• 207208211• 211• 212• 213213214xviD.16 The effective aperture of a trihedralcorner reflector with triangular panelsforincidence along the boresight and at an azimuth angleof 30 degrees 216D.17 Polarization response of thereceiving horn at a range of 11 m218D.18 Boresight response of a conventionaltrihedral corner reflector vs. range 218D.19 Azimuthal response pattern of aconventional trihedral corner reflector at rangesof 10 and 12 m219D.20 Azimuthal response patternof the model tower at a range of 11 m219xviiAcknowledgementsI wish to express my appreciation to my advisor, Prof. E.V.Jull; the members of my thesissupervisory committee, including Professors H.W. Dommel, W.McCutcheon, and L. Young;and Professors R.W. Donaldson, D.L. Pulfrey, L.M. Wedepohi, andM.M.Z. Kharadly for theirinterest and support during the past few years. Thanks arealso due to Helmut H. Lanziner,President of Offshore Systems Ltd., and Prof. Simon Haykin, Directorof the CommunicationsResearch Laboratory at McMaster University, for encouraging me to pursuea thesis projectconcerning the design of passive radar targets.I gratefully acknowledge the scholarships and awardswhich I have received during mycourse of study including a Science Council of British Columbia Graduate Scholarship,a BritishColumbia Advanced Systems Institute (ASI) Graduate Scholarship,an International Union ofRadio Science (URSI) Young Scientist Award, aProvince of British Columbia Graduate Scholarship, and a British Columbia Telephone Company Graduate Scholarship. Additionalfunding and support were provided by the Natural Sciences and EngineeringResearch Council ofCanada (NSERC) under Operating Grant 5-88571,Transport Canada (Transportation Development Centre, Montreal, P.Q.) under contract T8200-9-9560/O1-XSB, the NationalResearchCouncil of Canada (NRC) under Industrial Research Assistance Program (TRAP)Project No.9-8323-L-19, Offshore Systems Ltd. (North Vancouver, B.C.),Motorola Canada Ltd. (MDIDivision, Richmond, B.C.), and the University of British Columbia (Department ofElectricalEngineering).The assistance provided by Jean Liu, Kirk Jong, Robert Laing, David Clarke,and DonMacNeil in setting up the radar cross section measurement range, implementingthe digitalpattern recorder, constructing the prototype radar reflectors, and conducting theexperimentalmeasurements is greatly appreciated.xviiiChapter 1INTRODUCTION1.1 Background and MotivationRadars function by radiating electromagnetic energy and detecting thepresence and characterof the echoes returned by reflecting objects or targets. If these echoes are correlatedwiththe original transmitted signal, many of the characteristics of the targets can beestimatedincluding their range, bearing, apparent size or reflectivity, and certain aspects of theirphysicalgeometry. If a series of radar measurements is processed over time, the original estimateof thesecharacteristics can be refined and the future kinematic behaviour of the target can bepredicted.Radar has traditionally been associated with the detection and navigation of shipsand aircraft.More recently, it has become an important tool for remote sensing of the environment.Thechief problem of radar is to detect targets of interest and estimate their positionand physicalcharacteristics in the presence of interference from clutter returns and noise [1], [2].The scattering cross section a of a target is defined as the areaintercepting that amountof power which, when scattered isotropically, would produce anecho equal to that actuallyreturned by the target. Thus,E2HS2a lim 4irr2 lim 4irr2 (1.1)IE’12r—cc’HI2where E and W are the incident electric and magnetic fields, E8 andHsare the scatteredelectric and magnetic fields, and r is the range at which the scatteredfield is measured.Radar cross section or RCS refers to that portion of the scattering crosssection which isassociated with a specified polarization component of the scattered wave and is afunction ofthe size, shape, composition, and orientation of the target, the frequency of the incidentwave,and the polarization state of the radar transmitting and receiving antennas [3], [4].1Chapter 1. Introduction2The tendency for both natural and man-made objects to depolarize radar echoesin characteristic ways has been recognized since the earliest days of radar.It is convenient to describethe relationship between the polarization states of the incident and scatteredfields by a polarization scattering operator expressed in matrix form. Several representations arein commonuse. The polarization scattering matrix relates incident and scattered fields which havebeenexpressed as complex polarization vectors while the Mueller matrix relatesfields which havebeen expressed as Stokes vectors. Other forms, such as the Stokes scatteringoperator and thecovariance matrix, are used in certain methods forsynthesizing arbitrary polarization responsesfrom sets of experimental data. Measurement of the complete polarizationresponse of a targetrequires a radar which is capable of antenna polarizationcontrol or agility during transmissionand polarization diversity on reception. Although such radarsare considerably more complexthan their conventional counterparts, polarimetric radar signal processinghas become an important tool for target detection and classification in several fields including radar meteorology,geophysical remote sensing, and certain specialized forms of radar navigation[5]—[8j.It is often necessary to enhance the radar cross section of a cooperative targeteither toincrease the maximum range at which the target can be reliablydetected or to provide atarget with a known response which may be used to assist in radarcalibration and performanceverification. Although it is sometimes possible to achieve the desiredresult simply by makingminor modifications to the natural shape of the target or by disturbingthe current distributionon the surface of the body with discrete impedance loading, it is usuallymore convenient to makeuse of auxiliary devices such as corner reflectors, dielectric lenses, and retrodirective antennaarrays which have been designed specifically to present a large radar cross section overwideangular ranges. The characteristics and relative merits of the various types ofRCS enhancementdevices have been widely discussed in the literature[9]—{12]. The response characteristics whichare required of such devices are determined by several factorsincluding the distance betweenthe radar and the target, the reflectivity of the surrounding clutter, theresolution of the radarin range and azimuth, the combination of transmit and receive polarizationsemployed bythe radar, and the nature of the target detection or radar calibration algorithm. BeamwidthChapter 1. Introduction 3and pulse width limited radar resolution cells are shown in Figure 1.1. The probabilities ofdetection PD and false alarmFFAare determined by the value of the detection threshold andthe probability distribution functions of the target and clutter returns as suggested by theexample presented in Figure 1.2.xFigure 1.1: Detection of a point target in ground clutter where the radarresolution cell is eitherbeamwidth limited (entire ellipse) or pulse width limited (shadedportion of ellipse), a is theradar cross section of the target, r is the range tothe target, h is the height of the radar,8Hand8are the half-power beamwidths of the radar antenna in azimuthand elevation, t is thepulse duration, C is the speed of light, and c is thedepression angle. (after [2,p.84])P(x)Figure 1.2: Probability distribution functionsof clutter C and a signalembedded in clutterS+C and the corresponding probabilities of detectionPD and false alarmPFAfor a givendetection threshold. (after [2,p.42])hIcsin a°fAThresholdChapter 1. Introduction 4In recent years, the development of airborne and spaceborne imaging radar systems forgeophysical remote sensing and radar-assisted positioning systems for marine navigation hasled to a requirement for rugged yet inexpensive calibration targets and location markers whichpresent both a very large scattering cross section and a specified polarization response over awide angular range. The NASA/JPL synthetic aperture radar (SAR) calibration site shownin Figure 1.3 is similar to the calibration ranges which have been established by several otherresearch organizations to assist in the geometric, radiometric, and polarimetric calibration ofSAR imagery [13], [14]. The experimental radar-assisted positioning system depicted in Figures1.4 and 1.5 is being developed by Transport Canada as a method for allowing large vesselsnavigating in inland waterways, harbours, and harbour approaches to accurately determinetheir position with respect to cooperative shore-based targets in real time [15]—[19]. While therequirements for targets used in the calibration of airborne radars can be estimated with afair degree of confidence [12]—[14], scattering by terrain at grazing incidence has not been wellcharacterized [20]—[22] and it is not yet possible to give reliable estimates of the size of targetsrequired to achieve specified probabilities of detection and false alarm in applications such asradar-assisted positioning.Although active targets are physically compact and their scattering cross section, angularcoverage, and polarization response can be modified with relative ease, their usefulness is limitedby several factors including their requirement for an external power source, interaction betweentheir transmitting and receiving antennas which may lead to regenerative feedback and distortion of their response patterns, and the limited reliability and stability of active components.Despite their greater physical size and finer mechanical tolerances, passive targets provide abetter and more reliable alternative when it is necessary to install devices in remote locationsfor extended periods of time. In Figure 1.6, the relative sizes of corner reflectors which presentthe same maximum radar cross section (4500 m2) at a frequency of 10 GHz are compared. Inthe face of conflicting requirements for a target which presents a large response and wide angular coverage while retaining mechanical ruggedness and ease of manufacture, trihedral cornerreflectors frequently represent the best compromise.Chapter 1. IntroductionA A 6,T,8 TrlfledralsI.. .‘ 0,45 degree Dihedrals•0 0,45,90 degree PARC’s (LandTone Generators (L and C)0 Passive Receivers (L only)5Figure 1.3: NASA/JPL syntheticaperture radar(from [14,p.227])calibration site at Goldstone, California.Figure 1.4: Radar-assisted positioning with respect(from [15, front cover])to cooperative shore-based targets... pathGoldstone SARCaubratlon Site E4Chapter 1. Introduction6Building/I.-< Shore Station<- Reflector/// ////Radar AntennaFigure 1.5: Use of range and azimuthgates to isolate shore station reflectors fromsurroundingclutter.(a)(b) (e)3m.(c) 2m. -Figure 1.6: Relative size of corner reflectors which present thesame maximum radar crosssection (4500 m2) atf= 10 GHz. A spherical target whichpresents an equivalent responsewould have a diameter of over 75 metres. (a) Trihedral corner reflectorwith triangular panels,(b) Trihedral corner reflector with square panels. (c) Dihedral cornerreflector. (d) Bruderhedral(a cylindrical sector attached to a flat plate). (e) Top hat reflector.(f) Biconical reflector.Range Gate//Azimuth Gate/RadarBeamwidth/Shore StationReflector////(d)(t)Chapter 1. Introduction71.2 OutlineThis study considers several problems related to the analysis, design, andimplementation ofpassive targets including transformation of polarization descriptors betweencoordinate frames,modification of the angular coverage and polarization response of conventional trihedral cornerreflectors, design of top hat reflectors with specified response characteristics, and design ofreflection polarizers derived from conducting gratings with rectangular grooves.In Chapter 2, the problem of transforming representations of polarizationstate and polarization scattering operators between coordinate frames is considered. It is shown thatsuchtransformations correspond to rotation of the polarization basis by a prescribed angle which isa function of both the transformation matrix which relates the two coordinate frames and thedirection of propagation. Two methods for determining the angle of rotation are derivedforthe case in which the local vertical is defined by the direction in each frame. Algorithms fortransforming common polarization descriptors are presented.In Chapter 3, the problem of predicting the response and angular coverage of a trihedralcorner reflector with panels of completely arbitrary shape is considered. Formost purposes,only triple-bounce reflections from the interior of the reflector need be accounted for since theycompletely dominate the response for most directions of incidence. A simple yet robust RCSprediction algorithm which overcomes many of the limitations of previous work is described. Theresponse patterns of three-fold symmetric and bilaterally symmetric trihedral corner reflectorswith panels of various shapes are compared and design curves for realizing bilaterally symmetrictrihedral corners with specified response characteristics are given.In Chapter 4, a method for altering the polarization response of aconventional trihedralcorner reflector by adding conducting fins or corrugations of appropriate dimensions and orientation to one of its interior surfaces is proposed. Design curves for twist-polarizing andcircularly polarizing trihedral corner reflectors are given. Methods for realizing linear polarization selective trihedral corner reflectors using similar techniques are proposed. An algorithmfor predicting the contribution of triple-bounce reflections to the response of a depolarizingtrihedral corner reflector as a function of the direction of incidence and the orientation of theChapter 1. Introduction 8reflector is derived. Experimental results are presented which show that prototype depolarizingreflectors respond essentially as predicted.In Chapter 5, the results of this study are summarized and recommendations for furtherwork are offered.In Appendix A, the problem of designing top hat reflectors with specified response characteristics is considered. Expressions for the elevation response pattern, maximum scatteringcross section, angle of maximum response, 1 and 3 dB beamwidths, and angle of median response for 1 and 3 dB beamwidths are derived and used to generate design curves. The resultsare used to solve a sample design problem.In Appendix B, the problem of scattering by a conducting grating with rectangular groovesis ccrnsidered. Analytical solutions are derived for the cases of TM- and TE-polarized incidentwaves by mode-matching between the free space and groove regions. Procedures for determiningthe validity of numerical results are described and the problem of determining the minimumnumber of modes required to accurately represent the fields in each region is studied. Animplementation of the analytical solutions as a pair of subroutines coded in Fortran 77 ispresented.In Appendix C, the problem of modifying a conventional trihedral corner reflector to presenta circular polarization selective response is considered. It is shown that such a response cannotbe realized simply by using the techniques of Chapter 4 because the corresponding polarizationscattering matrix cannot be diagonalized. Alternative methods for obtaining such a responsebased on the addition of a suitable transmission polarizer to reflectors which present either alinear polarization selective or a twist polarizing response are proposed.In Appendix D, the experimental facility which was developed to measure the responseof prototype trihedral corner reflectors is briefly described. Details of its physical layout, thedesign and implementation of the CW radar apparatus and digital pattern recorder, and theresults of tests used to verify its suitability for use in the measurement program are given.Recommendations for future modifications and improvements are offered.References[1] F.E. Nathanson, Radar Design Principles: Signal Processing and the Environment, 2nd ed.New York: McGraw-Hill, 1991, pp. 1—48.[2] N. Levanon, Radar Principles. New York, Wiley, 1988,pp.1—97.[3] E.F. Knott, J.F. Shaeffer, and M.T. Tuley, Radar Cross Section: Its Prediction, Measurement, and Reduction. Norwood, MA: Artech House, 1985,pp.47—83.[4] A.K. Bhattacharyya and D.L. Sengupta, Radar Cross Section Analysis and Control. Nor-wood, MA: Artech House, 1991,pp.10—32.[5] D. Giuli, “Polarization diversity in radars,” Proc. IEEE, vol. 74, pp. 245—269, Feb. 1986.[6] A. Macikunas and S. Haykin, “Polarization as a Radar Discriminant,” in Selected Topics inSignal Processing, S. Haykin, Ed. Englewood Cliffs, NJ: Prentice-Hall, 1989,pp.251—286.[7] D. Atlas, Ed., Radar in Meteorology. Boston, MA: American Meteorology Society, 1990.[8] F. T. Ulaby and C. Elachi, Eds., Radar Polarimetry for GeoscienceApplications. Norwood,MA: Artech House, 1990.[9] D. R. Brown, R. J. Newman, and J. W. Crispin, Jr., “RCS Enhancement Devices,” inMethods of Radar Cross Section Analysis, J. W. Crispin, Jr. and K. M. Siegel, Eds.New York: Academic Press, 1969, pp. 237—280.[10] C. T. Ruck, “Radar Cross Section Enhancement,” in Radar Cross Section Handbook,G.T. Ruck, Ed. New York: Plenum Press, 1970, pp. 585—601.[11] F. T. Ulaby, R. K. Moore, and A. K. Fung, Microwave Remote Sensing: Active and Passive,vol. 2. Norwood, MA: Artech House, 1982, pp. 766—779.9References10[12] D. R. Brunfeldt and F. T. Ulaby, “Active reflector for radar calibration,”IEEE Trans.Geosci. Remote Sensing, vol. GE-22,pp.165—169, Mar. 1984.[13] S. H. Yueh, J. A. Kong, R. M. Barnes, andR. T. Shin, “Calibration of polarimetric radarsusing in-scene reflectors,” J. Electrornagn. Waves Appi., vol. 4,pp. 27—48, Jan. 1990.[14] A. Freeman, Y. Shen, and C.L. Warner, “PolarimetricSAR. calibration experiment usingactive radar reflectors,” IEEE Trans. Geosci.Remote Sensing, vol. GE-28,pp.224—240,Mar. 1990.[15] C. Stiles, R.O. Hewitt, and C.O. MdHale, “PRANS Trials: Evaluationof a Precise RadarNavigation System.” Transport Canada Pubi. No.TP 2800E. Montreal: TransportationDevelopment Centre, Aug. 1981.[16] D. Kalnicki and R. Harrs, “Performance Evaluation andDemonstration of the RANAV(Radar-Assisted Precise Navigation) System in the St. LawrenceRiver.” Transport CanadaPubi. No. TP 11326E. Montreal: Transportation DevelopmentCentre, April 1992.[17] S. Haykin, “Polarimetric radar for accurate navigation,”Can. J. Elect. Comp. Eng., vol. 17,pp.130—135, July 1992.[18] H. Lanziner, D. Michelson, S. Lachance, and D. Williams,“Experiences with a commercialECDIS,” mt. Hydrogr. Rev., vol. 67, no. 2,pp.69—86, July 1990.[19] D.C. Michelson et al., “Use of circular polarization in a marine radar positioningsystem,”Proc. IGARSS’89, (Vancouver, B.C.), July 1989.[20] M.W. Long, Radar Reflectivity of Land and Sea. Norwood, MA: ArtechHouse, 1983.[21] D.K. Barton, “Land clutter models for radar design andanalysis,” Proc. IEEE, vol. 73,pp. 198—204, Feb. 1985.[22] F.T. Ulaby and M.C. Dobson, Handbook of Radar Scattering Statisticsfor Terrain. Norwood, MA: Artech House, 1989.Chapter 2TRANSFORMATION OF POLARIZATIONDESCRIPTORSBETWEEN COORDINATEFRAMES2.1 IntroductionThe polarization state of anelectromagnetic wave is a vectorquantity which refers to thebehaviour with time of the electric fieldas observed at a fixed point in space.If the wave ismonochromatic, the tipof the electric field vector will trace anellipse in the plane orthogonalto the direction of propagation.Such a wave is said to be completelypolarized. If the wavecoiltains a random componentin amplitude or phase, it will occupy afinite bandwidth andthe polarization ellipse will tendto change shape and orientationwith time. Such a wave issaid to be partially polarized. Inthe extreme case, the tip of the electricfield vector will traceout a figure which istotally random in shape and the wave is saidto be randomly polarized(or unpolarized). The fundamentalaspects of wave polarization havebeen reviewed by severalauthors, e.g., [1]—[6].zAntenna coordinate framex”-\ /Local or body-fixed coordinate frameGlobal coordinatefra:Figure 2.1: Coordinate frames usedto define polarization statein radar scattering problems.11Chapter 2. Transformation of Polarization DescriptorsBetween Coordinate Frames 12Since polarization is a vector quantity, its descriptionmust be referred to a particularcoordinate frame. When problems involving propagationor scattering in the vicinity of theearth’s surface are considered, it is oftenconvenient to describe the polarization state of anantenna or the polarization response of a targetwith respect to a global coordinate frame inwhich the earth’s surface is coincident with the z-y planeand the local vertical is parallel tothe z axis. If the radiation or scattering characteristics of anobject can be determined moreefficiently in a different frame or if the object isfree to rotate about one or more axesas inthe case of airborne or spaceborneplatforms, it may be preferable to define thepolarizationcharacteristics of the object withrespect to a local or body-fixed coordinate frameinstead,as suggested by Figure 2.1. In turn, the responseof the device is measured in yet anothercoordinate frame which is defined by the antenna. However,with the exception of the specialcases considered by Mott [6] and Krichbaum [7], theproblem of transforming polarizationdescriptors between coordinate frames has received littleattention in the literature.In section 2.2, the concept of polarizationstate is briefly reviewed. In section 2.3, it isshown that transformation of a polarizationdescriptor between coordinate frames correspondsto rotation of its basis by a prescribed angle which isa function of both the transformationmatrix which relates the two coordinate frames andthe direction of propagation. Two methodsfor determining the angle of rotation for the case inwhich the local vertical is defined by thedirection in each frame are derivedusing spherical trigonometry and vector algebra, respectively.In section 2.3, methods for determining the elementsof the coordinate transformation matrixare reviewed. Although the matrixcan be determined from either the relativedirections ofthe three principal axes in each coordinateframe or the Euler angles which define a seriesof rotations which will transform onecoordinate frame into the other, in practice it maybedifficult to obtain these parameters. A thirdmethod is derived which permits theelements ofthe transformation matrix to be determined fromany pair of arbitrary directions whichhavebeen expressed in terms of both coordinateframes. In section 2.4, algorithms for rotatingthebasis of several common polarization descriptors arepresented.Chapter 2. Transformation of Polarization DescriptorsBetween Coordinate Frames 132.2 Wave PolarizationThe electric field vector E of a plane wave travelling in the direction k may becharacterized interms of a horizontally polarized component EHI and a vertically polarizedcomponent Ef’.Thus,E = (EH& + EvI)e_jk.r(2.1)where k = 27r/A is the propagation constant in free space and the timedependenceejwthasbeen suppressed. The amplitudes Ejq and Ev are complex quantitiesgiven byEH = aj, (2.2)= av, (2.3)where a and av are the magnitudes of EH and Ev, respectively, and H and 5vare their phaseangles. The coordinate frame (I, , k) can be specified in terms of thetriad (i,,)defined bya spherical coordinate system such thatk . .k = r = cosqsin8x+sinqsinOy+cos8z,(2.4)h = = —sinqx+cosqy,(2.5)xkj= kxh — = —cosqcos6—sinq5cos0+sin8,(2.6)as depicted in Figure 2.2. The definitions of I and Igiven by (2.5) and (2.6) are somewhatarbitrary and any pair of orthogonal directions which form a right hand triadwith k may besubstituted. In cases where the direction of propagation coincideswith the z axis, e.g., theantenna coordinate frame of Figure 2.1, the definitions of Ii and13 given by (2.5) and (2.6) areambiguous. If(2.7)it is common to define h and 13 such thatI. = = ±, (2.8)15 = kxh = +. (2.9)Chapter 2. Transformation of Polarization Descriptors BetweenCoordinate Frames 14When propagation occurs mainly in the vicinity of the z axis, it is often more convenient tosimply redefine polarization in terms of a new coordinate frame in which ‘, ‘, and ‘ correspondto ,, and , respectively, in the original frame and z’ defines the new local vertical.zIf the wave is monochromatic, the tip of the electric fieldvector will trace an ellipse in theplane orthogonal to the direction of propagation. The shape, sense of rotation,and orientationof the effipse are sufficient to specify the polarization state of the wave. Considera polarizationellipse with semi-major axis OA and semi-minor axis OB, as depicted in Figure 2.3. Themagnitude of the axial ratio R is given byIRH=g, (1IRIoo),According to the IEEE convention, the axial ratio is positive for right-hand polarizedwavesand negative for left-hand polarized waves while the reverse is true for the ellipticity angle.The tilt angle r is defined as the angle between the horizontal and the semi-major axis ofthepolarization ellipse and is valid over the range _90° r+900.The entire range of possiblepolarization states can be mapped onto the surface of a Poincaré sphere as shown in Figure 2.4.kryxFigure 2.2: Coordinate system for a plane wave propagating in the direction k.hwhile the ellipticity angle E is defined as= — cot1R, (—45° e +45°).(2.10)(2.11)Chapter 2. Transformation of PolarizationDescriptors Between CoordinateFrames 15A polarization state withellipticity angle E and tilt angle T correspondsto a point havinglongitude 2T and latitude 2c. Linear and circularpolarization states map onto the equatorandpoles, respectively, while left andright elliptical polarization states map onto theupper andlower hemispheres.Figure 2.3: A polarization ellipse with semi-majoraxis OA, semi-minor axis OB, and tilt angle r.EQUATORREPRESENTSLOWER HEMISPHERERIGHT-HAND SENSELONGmJOE =2X TILT ANGLEAhPolarizationellipseLATmJDE =2X ELLIPTICITYANGLECIRCULAR POLARIZATIONSUPPER HEMISPHERELEFT-HAND SENSELINEARFigure 2.4: Mapping of polarization states onto aPoincar sphere. (from [8],p.82)Chapter 2. Transformation ofPolarization Descriptors Between Coordinate Frames 162.3 Transformation of CoordinatesConsider two coordinate frames which are related by acombination of translation and rotationand let the directions and 2 define the localvertical in each frame. In cartesian coordinates,the position vectors r’ and R in the two frames are related by acoordinate transformationmatrix [T] such thatLm1 X — Xo=2m2 n2 Y —, (2.12)z’ £3 m3 n3 Z—Zowhere the origin 0’ of the x’y’z’ coordinate frame is locatedat (Xo, Y0,Zo) relative to the XYZcoordinate frame and l, m1,n1; 12, m2,n2; and 13,m3,n3 are the direction cosines of the x’, y’, z’axes relative to the X, Y, Z axes, respectively. Sincethe coordinate transformation matrix is aunitary matrix, its inverse is identicalto its transpose and the reverse transformation is simplygiven byXI XOy’+y0 . (2.13)z’ Z0The polarization state of a propagatingwave can be transformed from one coordinate frameto another by direct application ofeither (2.12) or (2.13) to the components of theelectricfield vector as described by Mott [6,pp.212—219). However, this technique is cumbersomeand cannot be easily generalized tothe many different’ methods which are used to representthe polarization state of a wave or the polarizationresponse of a scatterer. The limitations areparticularly apparent in cases where polarizationstate is described with respect to an ellipticallyor circularly polarized basis or where thewave is partially polarized. A more general approachis suggested by considering transformation ofthe corresponding polarization ellipse betweencoordinate frames. The ellipticity angle c isinvariant under either translation or rotation sinceit depends only on the magnitude of theaxial ratio and the polarization sense ofthe wave.Although the tilt angle r is invariant undertranslation, it will not be preserved under rotationxYZm1 m2 m3ni2fl3Chapter 2. Transformation of PolarizationDescriptors Between Coordinate Frames17cosb sini& 0[Tj=—sin cosb 00 0 1a = T — T, (—90° a+900), (2.15)where a corresponds to the angle between theunit vectors and ‘ which define the horizontalplane in each coordinate frame with respect tothe direction of propagation, as depicted inFigure 2.5. In general, it can beshown that transformation of any polarization descriptorbetween coordinate frames may be regardedas a change of basis transformation correspondingto rotation of the polarization basis by an angle aabout the direction of propagation.Figure 2.5: A polarization ellipse showing therelationship between the tilt angles rand r’ inthe xyz and x’y’z’ coordinate frames.unless the horizontal planes in bothcoordinate frames are parallel to each other. Thisconditionwill be satisfied only if the coordinatetransformation matrix is of the form(2.14)which corresponds to rotation about the z aidsby an angle b. Otherwise, the differencebetweenthe tilt angles of the polarization ellipse in thetwo coordinate frames is given by6Chapter 2. Transformation ofPolarization Descriptors Between CoordinateFrames 18The angle c between the unit vectorsand ‘ is a function of both the direction of propagation and the coordinate transformation matrixwhich relates the two coordinate frames. Thisvariation is apparent in Figures 2.6 and 2.7 where the parallels of twocoordinate frames whichare related by pure rotation are plotted on perspective andMercator projections, respectively.Since the Mercator projection is conformal, theangles between the parallels are accurately depicted at all points on the grid and the variationin the angle with the direction of propagationcan be easily visualized. In the context of radarcross section measurement, Krichbaum [7] hasderived an expression for this angle1for the special case in which the direction of propagation iscoincident with the z axis and the coordinate transformationcorresponds to rotation about thex and y axes. Here, two methods fordetermining the angle cr for any direction of propagationand coordinate transformation are derivedusing spherical trigonometry and vector algebra,respectively. Since translation between the coordinate framescan be neglected, it is convenientto definex=X—X0, y=Y—Yo, z=Z—Z0,(2.16)and reduce (2.9) and (2.10) toLm1 n1 x£2 m2n2 y, (2.17)z’ £3 m3 n3 zandx4£23y m1 m2m3 y’, (2.18)Zi 2713where the origins of the zyz and z’y’z’ coordinate frames arecoincident and l, m1,ni;12, m2,n2;and l3 m3,n3 are the direction cosines of thex’, y’, z’ axes relative to the x, y, z axes, respectively. Expressions for transforming directionexpressed in terms of the elevation angle 6 andazimuth angle 4 between coordinate frames can bederived from the coordinate transformations1which he refers to as the polarization angle T.Chapter 2. Transformation of Polarization Descriptors Between Coordinate Frames19Figure 2.6: A perspectiveby pure rotation.projection of the parallels of two coordinate frames which are related5.00,:2, 90.0175.0-180 -135 -90 -45 o 45 90 135 180•(deg)Figure 2.7: A Mercator projection of the parallels of two coordinate frames which are related bypure rotation. Since the projection is conformal, the angles between the parallels are accuratelydepicted at all points on the grid.I II I,I I/_________I______2Z5°-o(7 I\____I I I7.50I••_ I IChapter 2. Transformation of Polarization DescriptorsBetween Coordinate Frames 20(2.17) and (2.18) and the relationsr = 1, 8 = cos1z,4)= tan(y/x),219x = sin 6 cos4),y = sin9 sin4),z = cos 9.The forward transformation is given by9’ = cos1{sinS cos4) +m3 sinS sin4)+n3 cos 8] , (2.20)4)’— tan1[2sin9cos4)+msin9sin4)+ncos6]221— ii sin8 cos4)+m1 sin 9 sin4)+n1 cos 9] ‘while the reverse transformation is given by8 = cos1 [ni sin 0’ cos4)’+n2 sin 9’ sin4)’+n3 cos 9’] , (2.22)4)— tan’m1 sinS’ cos4)’ +m2 sin 0’ sin4)’+m3 cos ‘(2 23)L sin 9’ cos4)’+ £2sin 8’ sin4)’+£3 COS 6’Method IA general expression for the angle a may bederived using spherical trigonometry. Consider twocoordinate frames which are related by pure rotationas shown in Figure 2.6. Let the point 0 bethe common origin of the coordinate frames andlet the points Z, Z’, and P be the intersectionof the z axis, the z’ axis, and the direction of propagation witha unit sphere centered aboutthe origin. The points Z, Z’, and P define the spherical triangle ZZ’P.The great circle anglesdefined by the arcs Z’P, ZP, and ZZ’ are designated by z, z’, and p.The angles defined bythe vertices of the triangle opposite to arcs z, z’,and p, are designated by Z, Z’, and P. Sincethe vectors which are tangential to and 1 at P correspondto the directions O and ‘ andthe vectors 1, and,are orthogonal in both the x’y’z’ and xyz coordinate frames, it canbeshown that the angle P between the unitvectors 8 and ‘ is congruent to the angle a betweenthe unit vectors andChapter 2. Transformation ofPolarization DescriptorsBetween Coordinate Frames21In terms of sphericalcoordinates, the great circle anglesz, z’, and p are given byz = ,(2.24)=Op,(2.25)p ==,(2.26)while the vertex angles Z andZ’ are given byZ= —(2.27)= (2.28)If the point P lies on thegreat circle defined by the arc ZZ’,the vertex angles Z and Z’ vanishand a is given bya= 00 , 00<(Op — O) <180°,(2.29)= ±180°, —180°<(Op— 8) <00However, if the point P is coincidentwith either the point Z or its antipode,the local horizontalin the xyz coordinate framewill be undefined and the angle acannot be determined. Similarconsiderations apply if the pointP is coincident with either thepoint Z’ or its antipode. Inall other cases, an expression forthe angle a is obtained by applyingthe laws of sines to thespherical triangle ZZ’P togive eithersinpsinZsina ,(2.30)sin z— sinO sin(q5p— c5z’)— sinO’porsinp sinZ’ama = ,(2.31)sin z’— sinO’sin(qY— q5)— sinOpChapter 2. Transformation of PolarizationDescriptors Between CoordinateFrames 22then applying the law of cosines to givecosa = —cosZcosZ’+sinZsinZ’cosp,(2.32)=— cos(bp— ‘)cos(qV — qYp) +sin@p— )sin(—p)Cos 8.Since both the sine and cosine of the angle are known, it isa simple matter to determine theangle a using a four-quadrant arctangent function, e.g.,/sinp sin Z/sinzN— I\_cosZcosZl+sinZsinZlcosp)— tan’(sin8 sin(4p— c5’z)/sinO,— cos(p— ‘)cos(4/— ,)+ sin(qp— ‘)sin( —ç1,)cosMethod IIAn alternative expression for the angle a may bederived using vector algebra. Let cp definethe horizontal plane of the xyz coordinate framewith respect to the direction of propagationand be given byp=—sin4+cos,(2.34)and let74similarly define the horizontal plane of the x’y’z’coordinate frame and be given by(2.35)Let ip be the outward normal to the unit sphereat the point P in the xyz coordinate frameand be given by= — sin 8 cos+sin6 sin+cos 0 . (2.36)The expression for the unit vector in thex’y’z’ coordinate frame given by (2.35) must betransformed to the xyz coordinate frame. This canbe accomplished by determining qY in termsof 9 and4>using (2.21) and converting the basis of the vector using (2.18).The scalar tripleproduct of the unit vectors p, ,, and ipgivessina=p.(px4),(2.37)Chapter 2. Transformation of Polarization Descriptors Between Coordinate Frames 23while the dot product of the unit vectors c7p with , givescosa=p.. (2.38)Since both the sine and cosine of the angle are known, a can be determined using a four-quadrant arctangent function, i.e.,a = tan’(‘(2.39)If ip is coincident with z axis then the definition of p given by (2.34) is ambiguous. In suchcases, the direction of the horizontal plane with respect to the direction of propagation mustbe defined arbitrarily. Ifpis coincident with the z’ axis, similar considerations apply to thedefinition of, given by (2.35).2.4 Evaluation of the Coordinate Transformation MatrixIn order to apply the results presented in the previous section, it is necessary to determine theelements of the coordinate transformation matrix which relates the xyz and x’y’z’ coordinateframes according to (2.17) and (2.18). This can be accomplished if either the relative directionsof the basis vectors defined by the three principal axes in each coordinate frame or a series ofEuler angle rotations which will transform one coordinate frame into the other are known [9],[101.In the first case, the coordinate transformation matrix is given by4m1 n1[T]£2 m22 = • . ., (2.40)£3 m3 ri3 ‘• ‘• I’..where the two sets of basis vectors(i’, ‘, ‘)and(,,)have been expressed with respect to acommon coordinate frame. In the second case, the transformation is described in terms a seriesof angles through which the first frame can be rotated in order to bring it into coincidence withthe second. A maximum of three Euler angle rotations is sufficient to bring any two framesChapter 2. Transformation of Polarization DescriptorsBetween Coordinate Frames24into coincidence. One possibility is to begin the sequencewith rotation about the y axis by anangle4C,cost 0 sine x= 0 1 0, (2.41)z11 —sine 0 cos zfollow by rotation about the x” axis by an angle ,1 0 0= 0 cos i sin, (2.42)z” 0 — sin cos z”and conclude with rotation about the zaxis by angle ,cosC sinC0 x”=—Co, (2.43)z’ 0 0 1 z”as suggested by Krichbaum [7]. The product of the threerotation matrices given in (2.41),(2.42), and (2.43) yields the transformation matrixwhich relates the xyz and x’y’z’ coordinateframes,cos cos — sin sin ii sin sin(cos cosCsin + sin(sin7 cos[T]=—sinCcos— cos(sinisin cos(cos —sinCsin+ cosCsiniicos(2.44)—cosqsin —sing coscosSince rotation is not commutative, thisformulation is not unique and there are several othercombinations of Euler angle rotations which will yieldan equivalent coordinate transformation,e.g., [6], [9],[101.In practice, it may be difficult to obtain the parametersrequired by the basis vector andEuler angle methods for determining the elements of thecoordinate transformation matrix. Analternative method is derived here which allows thematrix to be determined from two arbitrarydirections (8k, qf1) and (02,2)which have been expressed in terms of both coordinateframes.Chapter 2. Transformation of PolarizationDescriptors Between Coordinate Frames 25The unit vectors which correspond to these two directionscan be determined by applying (2.19)to the ordered pairs (O,q4),(62,2),and(6,2)to yield= (xi,yi,zi) (a4,y,z) , (2.45)= (x2,y,z) E = (x,y,z) . (2.46)A third direction, r3, can be determined fromthe normalized cross products,— 1x2 ,______T3— , r3= ,2.47Ir Xr21 IrX rIin order to obtain two sets of vectors(,i, r) and(,,f) which form a linearly independentset in their respective coordinate frames, i.e.,X1 X2 X3 xc xYi Y2 J30, 0. (2.48)Z1 Z2Z3ZcZ ZFrom (2.17), it can be shown thatx x x £i m1 n1 z1 X2 X3c Y Y!3£2 m2fl2Yi Y2 /3(2.49)zcz z £3 m3 fl3 z1 z2 z3Solving for the coordinate transformation matrix gives—1L m1 n1 xc xx x1 x2 x3£2 m2 n2c Y2 Y!3 Yi Y2J3(2.50)£3 m3 n3zcZz!3zi z2 z3Since the coordinate matrices satisfy (2.48), theywill always have an inverse and (2.50) willalways have a valid solution.Chapter 2. Transformation of Polarization Descriptors Between Coordinate Frames262.5 Rotation of the Basis of Common Polarization DescriptorsIn this section, methods for rotating the basis of polarization coordinatesand the complexpolarization vector e.g., [1], [3], [11], are extended to moregeneral cases. Methods for rotatingthe basis of the Stokes vector and the coherency basis are derived. Someof the results presentedhere have recently been confirmed using a different approach by Mott[6].Polarization Coordinates and Complex Polarization RatiosThe entire range of possible polarization states of a completely polarized wave can be mappedonto the surface of a Poincaré sphere so that each polarization state is representedby itscoordinates in either latitude and longitude or elevation and azimuth. For example,a sphericalcoordinate system can be devised in which a polarization state W with ellipticity angle E andtilt angle r is represented by a point having longitude 2r and latitude 2e, as shownin Figure2.8. If the basis vectors Ii. and £‘ of (2.5) and (2.6) are rotated about the propagation vectorby an angle cr, the coordinates of the polarization state W’ in the new frame are givenbyC’ = C, (2.51)= r—c. (2.52)Alternatively, the coordinates of the polarization stateW may be described in terms of thepolarization angle y and phase angle 6 which are derived from the expression of thecorresponding plane wave as the weighted sum of orthogonally polarized basis states. For example, theelectric field vector E may be characterized in terms of horizontally and verticallypolarizedcomponents withE (EH,+ Evi’)e_ulcr,(2.53)where the time dependence €i has been suppressed andEH and Ev are complex amplitudesgiven byEH = aj-e3SH, (2.54)= av e5” , (2.55)Chapter 2. Transformation of Polarization Descriptors BetweenCoordinate Frames 27135°LINEARHORIZONTALLINEARLHCFigure 2.8: Poincaré sphere representation of the polarization state Wof a plane wave. (after[8],p.81)as described in section 2.2. The quantities aH and av are themagnitudes of EH and Ev,respectively, and6Hand5vare their phase angles. The polarization angle7Land phase angleare given by7L= tarf’(av/aH),=(—90°7L+90°)(—90°‘5L+900).(2.56)(2.57)Although any pair of orthogonally polarized basis states may be employed, the mostcommonare horizontal and vertical, 45° and 135° linear, and left and right circular. In thelatter cases,the polarization angle and phase angle are defined in a similar mannerto7Land6Lwith7D= tan1(a135/a45)613 6135645,-Ic= tan1(aR/aL),‘5C(—90° +90°),(—90° 5J3 +90°),(—90°7c+90°),(—90°6c+90°).(2.58)(2.59)(2.60)(2.61)LINEAR450LINEARRHCandChapter 2. Transformation of Polarization Descriptors Between Coordinate Frames 28In each case, a polarization state W with polarization angle y and phase angle is represented by a point having elevation 27 and azimuth 5. The orientation of the axes from whichand cc are measured is shown in Figure 2.8. The phase reference for orthogonal circularcomponents is defined in Figure 2.9.VERTICAL8CH OR ZON TA LFigure 2.9: Definition of the phase reference for orthogonal circular components. (from [8],p.82)Rotation of the basis of polarization coordinates is easily accomplished if they are expressedin terms of the ellipticity angle E and tilt angle r of the corresponding polarization ellipse.Expressions for transforming the polarization coordinates(7L, 6L), (7D, 6D),or(7c, 6c)into thepolarization coordinates (e, r) and back can be derived using spherical trigonometry. Considerthe right spherical triangle defined by the points corresponding to horizontal polarization, thepolarization state W, and the polarization state (0, T). For polarization coordinates expressedwith respect to a horizontally aiid vertically polarized basis, the forward transformation is givenby=sin1 (sin27LSIfl ccL) , (2.62)r = tan’(tan27L cos6L) , (2.63)while the reverse transformation is given by7L =cos1 (cos 2e cos 2T) , (2.64)= tan’(tan2c csc2T) (2.65)Chapter 2. Transformation of Polarization DescriptorsBetween Coordinate Frames 29Next, consider the right sphericaltriangle defined by the points corresponding to 45°linearpolarization, the polarization state W(€, r),and the polarization state (0, r). Forpolarizationcoordinates expressed with respect to a 45°and 135° linearly polarized basis, the forwardtransformation is given by= sin (sifl27D SiflD) ,(2.66)T= _tan_1(cot27Dsec6D), (2.67)while the reverse transformation is given by=cos1(cos 2€ sin 2r) ,(2.68)tan1 (tan2E sec2r) .(2.69)Finally, consider the right spherical triangledefined by the points corresponding to left circularpolarization (45°,0), the polarization state W(€,r),and the polarization state (€,0). For polarization coordinates expressed with respect to a leftand right circularly polarized basis, theforward transformation is given by= tan1(1)— -Ic, (2.70)r=t5cj/2, (2.71)while the reverse transformation is given by7c= tan’(l)—c, (2.72)= 2r. (2.73)The surface of the Poincar spherecan be mapped onto the entire complexplane by astereographic projection. Each polarizationstate is represented by a complex numberwhich isreferred to as the complex polarizationratio p. Rotation of the basis of the complexpolarizationratio is easily accomplished if the ellipticityangle e and tilt angle i- of the correspondingChapter 2. Transformationof Polarization Descriptors Between CoordinateFrames 30polarization ellipse are known.The complex polarization ratio is relatedto the polarizationcoordinates(7L, 6L), (7D,5I3), andQyc, c)byPL = PL(2.74)PD = PD, (2.75)pc = ,°c, (2.76)where the polarization ratio p isdefined byPL= tan7L av/aH,(2.77)PD= tan’yDa135/a45, (2.78)Pctan7c aR/aL,(2.79)and the phase angle is definedby (2.57), (2.59), and (2.61). Togetherwith (2.62)—(2.73), thisis sufficient to define the transformationof the complex polarization ratiosPL, PD, pcto thepolarization coordinates (€, r) and back.Complex Polarization VectorsThe pair of complex amplitudeswhich arise from the representation of a planewave as theweighted sum of orthogonally polarizedbasis states may be arranged to yield a complexpolarization vector, e.g.,EHE = ,(2.80)where EH and Ev are the complex amplitudeswhich correspond to a horizontally and verticallypolarized basis. Rotation of the polarizationbasis about the propagation vector k by ananglea is accomplished by applicationof a rotation operator {R] to the complexpolarization vectorE to yield= [R] E.(2.81)Chapter 2. Transformation ofPolarization Descriptors Between Coordinate Frames31If the basis of the complex polarizationvector is linear, the corresponding rotation operator[RUis simply given bycosa —sina[RL]. (2.82)sina cosaThus,E cosa —sina EH, (2.83)sina cosaandE5 cos a — sin a E45• (2.84)E35 sin a cos a E135If the basis of the complex polarization vector is eitherelliptical or circular, it must be transformed to linear before the rotation operator given by (2.82)is applied. After it has beenrotated, the polarization basis can be restoredto its original effipticity by the reverse transformation. This transformation can be performed byan ellipticity operator,coSE jsinc[H] = ,(2.85)jsinc coscwhich will modify the ellipticity angle of the polarizationbasis states which define the upperand lower element in the polarization vector by € and-, respectively, without affecting theirrespective tilt angles [11]. Let be the ellipticityangle of the polarization basis state whichdefines the first element in the polarizationvector. Since the basis states are orthogonal,the ellipticity angle of the polarization basis state whichdefines the second element in thepolarization vector is given by -e. Thus, thegeneral rotation operator [R] for the complexpolarization vector is given by[R] = [H(c)] [RL(a)] {IJ(—€)](2.86)cos e j sin c cos a — sin a cos e—jsinEj sin c cosE sin a cos a—jsin c cos eChapter 2. Transformation of Polarization Descriptors Between Coordinate Frames 32If the basis is circularly polarized, e = ir/4 and the corresponding rotation operator [Rc] isgiven bye 0[Rd = . (2.87)0 e2Thus,Eeic0 EL=. (2.88)E 0 eiERThe Stokes VectorThe .Stokes vector representation of a completely polarized wave is given bySo I(IEvI2+ IEHI2)1= S1= Q = *(IEvI2— IEHI2)= 1cos2ecos2T(2.89)S2 UIEvIIEHIcos6Lcos2csin2r53 VjIEvIjEHIsin5Lsin2cwhere EH, Ev, and6Lare defined in (2.54), (2.55), and (2.57),‘iis the impedance of free space,1 is the total power carried by the wave, and c and T are the ellipticity and tilt angles of thecorresponding polarization ellipse. For a completely polarized wave,s=s?+s+s.(2.90)The normalized Stokes vector is given by0181 cos2Ecos2rs = (2.91)cos2€sin2r83 sin2ewhere1 = s + s + s . (2.92)Chapter 2. Transformation ofPolarization Descriptors Between Coordinate Frames33Since all the elements of the Stokesvector are expressed in units of power, polarimetric datawhich are expressed in Stokes format can bespatially and temporally averaged with relativeease. Also, the elements of the Stokes vector are alwaysexpressed in real numbers so recourse tocomplex arithmetic is not required. Unlikethe complex polarization vector, the Stokesvectorcan also represent the polarization state of quasi-monochromaticor partially polarized wave.In such cases, the Stokes vector can be resolved into a completelypolarized component S, andan unpolarized component S, such thatS=S+S.(2.93)Then,S0 l—d1 10 cos2€cos2rdcos2fcos2rS= =10 +d =1, (2.94)S2 0 cos2Esin2Tdcos2Esin2TS3 0 sin2€dsin2cwhere d is the ratio of the power carried by thepolarized component of wave S, to the totalpower carried by the wave and is given by+ S + Sd= , 0dl.(2.95)soThe equivalent normalized Stokes vector is given by1—d 1 10 cos2ecos2T dcos2€cos2rs =, (2.96)0 cos 2E sin 2T dcos 2e sin 2r83 0 sin2€ dsin2Ewhere the depolarization ratio d is given byd=/4+$+s, 0dl.(2.97)Chapter 2. Transformation of PolarizationDescriptors Between Coordinate Frames 34Consider rotation of the polarization basisabout the propagation vector by an angle a.Let S represent the Stokes vector in theoriginal coordinate frame and let 5’ represent theStokes vector in the new coordinate frame. From (2.94), S’is given byS61S dcos2cos2r’= Io. (2.98)dcos2sin2r’dsin2€From (2.15), r’ = r — a, and1S dcos2cos2(T — a)1(2.99)S dcos2csin2(r — a)dsin2c1d cos 2e(cos 2T cos 2a + sin 2r sin 2a)10d cos 2E(sin 2r cos 2a — cos 2T sin 2u)dsin2cBy inspection, (2.99) can be factored to yield1 0 0 0 10 cos2a sin2a 0dcos2€cos2r= 1(2.100)0 — sin 2a cos 2a 0 dcos 2e sin 2T0 0 0 1dsin2cThus, S’ is related to S by1 0 0 0 So0 cos 2a sin 2a 0(2.101)0 —sin2a cos2a 0 S20 0 0 1 S3Chapter 2. Transformation of Polarization Descriptors Between Coordinate Frames35It is convenient to express (2.101) in the form[5’]= [Rs] [5], (2.102)where[Rs]is a rotation operator given by1 0 00o cos2c sin2o 0[Rs] = . (2.103)0 —sin2c cos2c 00 0 0 1A modified form of the Stokes vector is sometimes used to simplify theformulation of radiative transfer problems, e.g., [4], [12]. From (2.89), expressionsfor the vertically and horizontallypolarized intensity,Il-I= _IEvI2(I+Q)/2, (2.104)Iv= _IEHI2= (I—Q)/2, (2.105)are obtained. These expressions are substituted for the first twoelements of the Stokes vectorto yield the modified Stokes vectorSm,SmO -1H (jEvI2) (1+dcos2ecos2r)Sml — Iv — —(IEHI2)(1—dcos2Ecos2r), (2106)5m2UIEvIIEHlcos6L dcos2Ecos2TSm3VIEvWEHIsinL dsin2Ewhich is related to the conventional Stokes vector S byC, 1 1 i Cm0U U“0Smi 0 0 Si(2.107)5m20 0 1 0 S25m30 0 0 1 S3Consider rotation of the polarization basis aboutthe propagation vector A by an angle . LetSmrepresent the modified Stokes vector in theoriginal coordinate frame and let S representChapter 2. Transformation of PolarizationDescriptors Between Coordinate Frames 36the modified Stokes vector in the newcoordinate frame. From (2.101) and (2.107), it canbeshown that the rotation operator [Rm] isgiven by00 1 0 0 000_i1 _10 0 0 cos2a sin2a 01 _10 0— 2 22 2114’mi —0 0 1 0 0 —sin2acos2a 0 0 0 1 00 001 0 0 01 0 001cos2 a sin2 a sin 2a 0sin2 a cos2 a_lsin 2a 02(2.108)— sin 2a sin 2a cos 2a 00 0 0 1Thus, S is related to S byS0 cos2 a sin2 asin 2a 0SmO= sin2a cos2a —sin2a 0Smi(2.109)— sin 2a sin 2a cos 2a 0 Sm2C! 1 C‘m3U U U m3A third variant of the Stokes vector is usedto define the data format employed by the JPLpolarimetric imager [4]. The JPL Stokesvector S3 is given by*(JEHI2+ lEvi2)1S31 = j(lEHl2— lEvi2)= 10—dcos2€cos2r(2.110)5j2 —1IEvilEHicost5L—dcos2Esin2r5j3 EvllEHIsinbLdsmn2Eand is related to the conventional Stokes vector S bySjO1 0 0 0SOSi’0—1 0 0Si=(2.111)5j20 0—10 S2S33 0 0 0 1 53Chapter 2. Transformation ofPolarization Descriptors Between CoordinateFrames 37Consider rotation of thepolarization basis about the propagationvector k by an angle a.Let S3 represent the JPL Stokesvector in the original coordinate frame and letS representthe JPL Stokes vector in thenew coordinate frame. From (2.103)and (2.111), it can be shownthat the rotation operator [R3]is given by—11 000 1 0 0 01 0000 —1 0 0 0cos2a sin2a 0 0 —1 00[R3]=0 0 —1 0 0 —sin2acos2a 0 0 0 —1 00001 0 0 01 00011 0 0 00 cos2a sin2a 0V(2.112)0 —sin2a cos2a 00 0 0 1Thus, S is related to S3 byi 0 0 0 s30S’ 0 cos 2a sin 2a 0 Si1(2.113)S’2 0 — sin 2a cos 2a 0 Sj2S’3 0 0 0 1The Coherency MatrixThe coherency matrix is anothermethod for representing the polarizationstate of a partiallypolarized wave which is sometimesused [5], [6], [13]. In terms of the Stokes vector, theelementsof the coherency matrix aregiven by[J}11 J12 = (So+S1)(S2+jS3). (2.114)J21 J22(S2—iS3) (SoSi)Substitution of the expressions for S0, S, S2,and S3 given in (2.94) into (2.114)yields anexpression for the coherency matrix in termsof the wave intensity Io, the ellipticity anglee andChapter 2. Transformation of Polarization Descriptors Between Coordinate Frames 38tilt angle T of the po’arized component of the wave and the depolarization ratio d,[J] 1(1+dcos2ecos2r) d(cos2Esin2r+jsin 2e)(2 115)d(cos2esin2r —jsin2c) (1 — dcos2ecos2r)The coherency matrix can be resolved into polarized and unpolarized components,{J]= {J] +[Jr] , (2.116)to give[J]—Io0 1—d+Iod(1+cos2os2T) (cos2Esin2r+jsin2E)1—d 0 (cos2esin2T—jsin2) (1—cos2€cos2r)(2.117)Consider rotation of the polarization basis about the propagation vector k by an angle o.Let [J] represent the coherency matrix in the original coordinate frame and let [J’] representthe coherency matrix in the new coordinate frame. Since the unpolarized component[J]isinvariant under rotation of the polarization basis,J1l J,’412= Jul1 Jul2, (2.118)J2l,‘L22J2lJ22and only transformation of the polarized component [Jr] need be considered. In the newcoordinate frame, [J] is given byJ71 J72 (1+cos2€cos2r’) (cos2Esin2r’ +jsin2E)JfJ= 10d(cos2csin2r’—jsin2E) (1—dcos2ccos2r’). (2.119)From (2.15), r’ = r — cr andJ1 J12 (1+cos2ccos2Qr—cr))(cos2Esin2(r—c) +jsin2c)=10dJ21 J22 (cos2esin2(r—a) —jsin2) (1 — cos2ccos2(T—a))(2.120)With a little effort, (2.120) can be factored to yieldJ1 J12 — cosa sina J,11 J12 cosa —sina(2121)J,21 J22 —sina cosa J21422sina cosaChapter 2. Transformation of PolarizationDescriptors Between Coordinate Frames39Thus, from (2.118) and (2.121), it can be shown that [J’jis related to [J] byil ‘i2= cosa sina J11 J12 CO5O. (2.122)J1 J2 —sinacosa J21 J22 sina cosaPolarization Scattering OperatorsThe scattering cross section a of a target is defined as the area interceptingthat amount of powerwhich, when scattered isotropically, would produce an echo equal to thatactually returned bythe target. Thus,1Es12a = urn 4irr2 , (2.123)r—ooIE’Iwhere E is the incident electric field,ESis the scattered electric field, andr is the range atwhich the scattered field is measured. Radar cross section refers to thatportion of the scatteringcross section which is associated with a specified polarizationcomponent of the scattered waveand is a function of the size, shape, composition,and orientation of the target, the frequency ofthe incident wave, and the polarization state of the radartransmitting and receiving antennas.The relationship between the polarization statesof the incident and scattered fields canbe described by a polarization scattering operatorexpressed in matrix form. Following thedefinition of scattering cross section presented in (2.123),the polarization scattering matrix [S]relates incident and scattered fields which have been expressed as complexpolarization vectors,i.e.,Ef S11 S12 El= . (2.124)E S21522EWhen defined with respect to a horizontally and vertically polarizedbasis, the polarizationscattering matrix can also be used to relate incidentand scattered fields which have beenexpressed as coherency matrices, i.e.,Jr1 Jf2 = SHH SHyhJ25IHT(2.125)47rr2SVH SvvJ J2Sjj S17Chapter 2. Transformation of Polarization DescriptorsBetween Coordinate Frames 40The Mueller matrix [L] relates incident and scattered fieldswhich have been expressed as Stokesvectors, i.e.,L11 L12 L13 L14 SS= 1L21 L22 L23 L24S• (2.126)s4irrL31 L32 L33 L34 SS L41 L42 L43 L44 SThere are three principal variants of the Mueller matrixwhich correspond to the conventional,modified, and JPL Stokes vectors, respectively. Inpractice, polarization scattering matrices andMueller matrices are often normalized by factoringout the scattering cross section of the targetand the range dependence of the response. Otherpolarization scattering operators which arederived from the polarization scattering matrixand the Mueller matrix, such as the covariancematrix and the Stokes scattering operator, are usedin computationally efficient methods forsynthesizing arbitrary polarization responses from experimentaldata [4].it is convenient to describe scattering problems with respectto a coordinate frame which iscentered on the scatterer. The local coordinate systemsused to define the polarization state ofthe incident and scattered fields are specified in a mannersimilar to that presented in section 2.2for the case of a single propagating wave. However, itis also necessary to specify the relationship between the local coordinate systems.According to the forward scatter alignment (FSA)convention, the propagation vectors of the incidentand scattered fields are aligned with the direction of propagation while accordingto backscatter alignment (BSA) convention, they alwayspoint towards the scatterer. While the expressions forthe incident field are identical underboth conventions, the expressions for the scatteredfield and, by extension, the correspondingpolarization scattering operator, are not.The coordinate systems and scattering geometriescorresponding to the FSA and BSAconventions are depicted in Figures 2.10 and 2.11, respectively. The subscripts i and s refer to fields expressedwith respect to the FSA convention whilethe subscripts t and r refer to fields expressedwith respect to the BSA convention. The unitChapter 2. Transformation ofPolarization Descriptors Between CoordinateFramesAxAAy41Figure 2.10: Coordinatesystems and scattering geometry for theforward scattering alignment(FSA) convention. (after[4], p. 18])Figure 2.11: Coordinate systemsand scattering geometry for the backscatteralignment (BSA)convention. (after [4], p. 18])AvsAht/ krIVr/ I/ _JChapter 2. Transformation of Polarization Descriptors Between CoordinateFrames 42vectors in the local coordinate system for the incident (or transmitted) field aregiven by= = E — = —cosqsin02 —sinsin0 —cos92., (2.127)xk= —,. = —sinx+cosqy, (2.128)Izxklfitf kxh —Ô = —cosq5cos8—sinqcos8+sin8, (2.129)while the corresponding unit vectors for the scattered (orreceived) field are given by= —•k = E — = —cos5sin68—sin4sinO3—cos68, (2.130)• xkhr =—h3 = —q =sinq53x—cos4i’y, (2.131)Izxkl= k x h E —8 = —cosq5cos68—sin4cosO+sin68. (2.132)Since the coordinate systems lii,fit)and(kr, ‘zr,fi,.) are coincident when the transmittingand receiving antennas are collocated, the BSA convention is a particularlyconvenient choicefor use in radar scattering problems. Unless otherwise stated, the BSAconvention will be theconvention used in the remainder of this study.Once the local coordinate systems for the incident andscattered fields have been defined,it is a simple matter to apply the results derived earlier in this section to theproblem ofrotating the basis of either the polarization scattering matrix or the Mueller matrixabout theradial vector i by an angle c. Since the propagationvector = —1, this transformation isequivalent to rotating the polarization basis about k by an angle —tx. Let [S]represent thenormalized polarization scattering matrix in the original coordinate frame andlet [S’] representthe normalized polarization scattering matrix in the new coordinate frame. Thecorrespondingscattering equations are given byES= [Sj E , (2.133)= [S’]Et’. (2.134)Chapter 2. Transformationof Polarization DescriptorsBetween Coordinate Frames43Since the rotationoperator [R] is a unitary matrix,= [R(—a)](2.135)andEs’= [RJ’ES,(2.136)E” = [RI—’ E’ ,(2.137)where expressionsfor the rotation operator [RI whichare appropriate for use withcomplexpolarization vectors having linearlypolarized, circularly polarized,or arbitrarily polarized basisstates are given by (2.82), (2.87),and (2.86), respectively. Multiplyingboth sides of (2.136)and (2.137) by [RI yields= [RIE’,(2138)Et = [Rj E’ .(2.139)Substituting (2.138) and 2.(139)into (2.133) gives[RIEs’= [SI [RIEt’.(2.140)Multiplying both sides by [R1’gives[Rj’ [R]Es’= [R]’ [SI [RIEt’, (2.141)which simplifies to= [RI’ [SI [RIEt’.(2.142)Thus, by equating (2.134) and (2.142),it can be shown that [S’I and[SIare related by[S’I= [Rj’ [5] [RI .(2.143)In a similar fashion, itcan be be shown that a Muellermatrix [L’j in the new coordinateframeis related to the Muellermatrix [L] in the original coordinateframe by[L’}= [RsI1 [LI[Rs](2.144)where expressions for the rotation operator [RsIwhich are appropriate for usewith conventional,modified, and JPL Stokes vectors aregiven by (2.94), (2.106), and (2.110),respectively.Chapter 2. Transformation ofPolarization Descriptors Between CoordinateFrames 442.6 ConclusionsIt has been shown that transformation ofa polarization descriptor betweencoordinate framescorresponds to rotation of its polarizationbasis by a prescribed anglewhich is a function ofboth the transformation matrixthat relates the two coordinate framesand the direction ofpropagation. Two methods fordetermining the angle of rotation for thecase in which thelocal vertical is defined by thedirection in each frame have been derivedusing sphericaltrigonometry and vector algebra,respectively. Both methods are robustand will yield thecorrect result but the method basedon vector algebra is more compactand would be easierto implement in software. Althoughthe elements of the coordinatetransformation matrix canbe determined from either the relativedirections of the three principal axesin each coordinateframe or the Euler angles which define aseries of rotations which will transformone coordinateframe into the other, in practice itmay be difficult to obtain these parameters.A third methodhas been derived which overcomes thislimitation by allowing the elements of thecoordinatetransformation matrix to be determinedfrom any pair of directions whichhave been expressedin terms of both coordinate frames.Algorithms for rotating the basisof several commonly usedpolarization descriptors, includingpolarization coordinates, the complexpolarization ratio, thecomplex polarization vector, the Stokesvector and several of its variants,the coherency matrix,the polarization scattering matrix,and the Mueller matrix have beenderived.References[1] R.M.A. Azzamand N.M. Bashara,Ellipsometry and PolarizedLight. Amsterdam:North-Holland, 1977.[2] S.R. Cloude,“Polarimetric techniquesin radar signalprocessing,” MicrowaveJ., vol. 26,no. 7,pp.119—127, July1983.[3] J.D. Kraus,Antennas, 2nd ed. NewYork: McGraw-Hill,1988,pp.70—81.[4] F.T. Ulabyand C. Elachi, Eds.,Radar Polarimetryfor Geoscience Applications.Norwood,MA: Artech House,1990,pp.1—52.[5] W.L. Stutzman,Polarization in ElectromagneticSystems.Norwood, MA: ArtechHouse,1992.[6] H.A. Mott,Antennas for Radarand Communications:A PolarimetricApproach. NewYork: Wiley, 1992.[7] C.K. Krichbaum,“Radar cross-sectionmeasurements,” inRadar Cross SectionHandbook.(G.T. Rucket al., Eds.) vol. 2,New York: Plenum,1970,pp.893—896.[8] IEEE StandardTest Procedums forAntennas. (ANSI/IEEEStd 149—1979.)New York:IEEE, 1979.[9] W.T. Thomson,Introduction to SpaceDynamics. NewYork: Wylie, 1961,pp.33—37.[10] R.R. Bate,D.D. Mueller, andJ.E. White, Fundamentalsof Astrodynamics.New York:Dover, 1971,pp.74—83.[11] S.H. Bickel,“Some invariantproperties of thepolarization scatteringmatrix,” Proc. IEEE,vol. 53,pp.1070—1072, Aug. 1965.45References46[12] S. Chandrasekhar, Radiative Transfer. New York:Dover, 1960,pp.24—35.[13] M. Born and E. Wolf, Principles of Optics. New York: Pergamon, 1965, pp.545—553.Chapter 33.1 IntroductionTRUNCATION AND COMPENSATION OFTRIHEDRAL CORNER REFLECTORSA trihedral corner reflector is a reentrant structure formed by the intersection of threemutuallyorthogonal reflecting panels. In general, a ray incident upon one of its interior surfaces willundergo reflection from each of the others in succession and will be returned tothe source.Although other scattering mechanisms contribute to the response, triple-bounce reflections fromthe interior of the reflector dominate over most directions of incidence. Sincetrihedral cornerreflectors present a large scattering cross section over a wide angular range, are mechanicallyrugged, and can be manufactured with relative ease, they are widely used in radar navigationand remote sensing as location markers and calibration targets. The relative sizes of trihedralcorner reflectors in common use are compared in Figure 3.1.(a)3m.2m±1m0mFigure 3.1: Relative sizes of trihedral corner reflectors used as location markers andcalibrationtargets in (a) radar navigation, (b) radar-assisted positioning, and (b/c) radar remote sensing.(b)(c)47Chapter 3. Truncation and Compensation of Trihedral Corner Reflectors 48The dependence of the scattering cross section and angular coverage of a trihedral cornerreflector on the size and shape of its reflecting panels has been recognized since the adventof radar. In the 1940’s and 50’s, closed-form expressions for the contribution of triple-bouncereflections to the response of trihedral corner reflectors with triangular, elliptical or rectangularpanels with corners of arbitrary length were derived [1]—[3]. Although it was apparent that awide variety of response characteristics could be obtained by appropriate shaping of the reflecting panels, a procedure which Robertson [4], [5] referred to as truncation and compensation,work in this area was not pursued due to the lack of either suitable methods for determining theresponse of a trihedral corner reflector with panels of completely arbitrary shape or a need forphysically large targets which would benefit from such modifications. In recent years, interestin altering the response of trihedral corner reflectors in this manner has been renewed by arequirement for physically large targets to serve as location markers and calibration targetsin radar navigation and remote sensing [6]—[9]. However, very little design data and relatedmaterial to guide the development of such reflectors are available in the literature.In section 3.2, the problem of predicting the response of a trihedral corner reflector withpanels of completely arbitrary shape is considered and an efficient and robust numerical methodfor solving Robertson’s model for the equivalent flat plate area of a trihedral corner reflectoris proposed. In section 3.3, the response characteristics of a selected set of trihedral cornerreflectors which present three-fold symmetry are compared. In section 3.4, the response characteristics of trihedral corner reflectors which present bilateral symmetry are considered andthe possibility of increasing the beamwidth of the response of such reflectors in one principalplane relative to the beamwidth in the orthogonal plane by modifying the size and shape ofthe reflecting panels in a suitable manner is examined. Design curves for bilaterally symmetricreflectors which are composed solely of triangular, elliptical, or rectangular reflecting panels aregiven. The response characteristics of a selected set of bilaterally symmetric reflectors whichare composed of combinations of panels with various shapes including triangular, circular, andsquare are compared. A related problem, the design of top hat reflectors with specified responsecharacteristics, is considered in Appendix A.Chapter 3. Truncation and Compensation of Trihedral Corner Reflectors 493.2 AnalysisA rigorous solution for the scattering cross section of a trihedral corner reflector must accountfor the contributions of single, double, and triple-bounce reflections from the interior of thereflector, deviations of the reflecting panels from perfect flatness and mutual orthogonality, anddiffraction by the panel edges. Numerical techniques such as the finite-difference time-domain(FD-TD) and the shooting and bouncing ray (SBR) methods have been successfully applied tothe problem and can account for most contributions to the response. However, calculating theresponse of a large target is extremely demanding and access to some type of supercomputeror massively parallel processor is generally required [10], [11]. In the case of an ideal reflectorwith reflecting panels which are perfectly flat and mutually orthogonal, the problem can besimplified considerably. A reasonably complete solution can be obtained by using physicaloptics (P0) to account for the contribution of reflections from the interior of the reflector whileusing the method of equivalent currents (MEC) to account for first order diffraction from theedges [12]. Alternatively, a hybrid approach which permits application of the Uniform Theoryof Diffraction (UTD) to the problem can be employed [13]. However, neither of these techniquescan be easily applied to reflectors with panels of completely arbitrary shape.For the purposes of designing trihedral corner reflectors with specified response characteristics, it is usually sufficient to account for the contribution of triple-bounce reflections from theinterior since they completely dominate the response for most directions of incidence. If thereflector is ideal, a ray which is incident upon one of the interior surfaces will generally undergoreflection from each of the others in succession and will be returned to the source. However,the reflecting panels are of finite extent and some rays will fail to intercept one or more of thepanels and will be lost. The equivalent flat plate area A of the reflector can be determined bylaunching a set of parallel rays towards the target, tracing each ray as it is reflected by eachof the interior surfaces, and projecting that portion of the reflector which contributes to thebackscatter response onto a view plane which is normal to the direction of incidence, as suggested by Figure 3.2. The scattering cross section a of the reflector is related to its equivalentChapter 3. Truncation and Compensation of Trihedral Corner Reflectors 50(a) (b)Figure 3.2: Scattering by a trihedral corner reflector with triangular panels.(a) Alternativeray paths. Path 1-1’ represents a ray which has undergone triple-bounce reflection back to thesource while path 2-2’ represents a ray which has undergone double-bounce reflection and beenscattered in a different direction. (b) The equivalent flat plat area of the reflector for incidencealong the symmetry axis. (c) The equivalent flat plat area for incidence off the symmetry axis.(from [5],p.13-11)xyFigure 3.3: Problem geometry and coordinate system forscattering by a trihedral corner reflector composed of triangular, elliptical, or rectangular panels withcorners of arbitrary length.equivalent flat plate area(c)zChapter 3. Truncation and Compensation of Trihedral Corner Reflectors 51A’,encepupilFigure 3.4: Spencer’s model for the equivalent flat plate area of a trihedral corner reflector.In this example, the intersection of the entrance pupil ABC and exit pupil A’B’C’ gives theequivalent flat plat area of a trihedral corner reflector with triangular panels for incidence offthe symmetry axis.flat plate area A by the physical optics relation,A2(3.1)where A is the wavelength of the incident wave.It is convenient to describe scattering by a trihedral corner reflector with respect to thecoordinate frame shown in Figure 3.3. Spencer [1] empirically derived a simple geometricmodel for predicting the equivalent flat plate area of an ideal trihedral corner reflector based onexperiments that he conducted with reflectors fabricated from optical mirrors. In the model,the polygon which defines the outside edges of the reflecting panels and its inverted image areprojected onto a view plane which is normal to the direction of incidence. The inverted imageis obtained by projecting the original polygon through the apex of the reflector. The projectionand its inverted image are referred to as the entrance pupil and exit pupil of the reflector,respectively. An example for the case of a reflector with triangular panels is shown in Figure3.4. According to the model, the equivalent flat plate area of the reflector is the area commonto the two pupils. Using this model, Spencer derived closed-form expressions for the response oftrihedral corner reflectors composed of triangular and square panels with equal corner lengths.exit pupilAChapter 3. Truncationand Compensation of TrihedralCorner Reflectors52Later, these were extended to thecase of trihedral corner reflectorscomposed of triangular,elliptical, or rectangular panelswith unequal corner lengths by Siegelet al. [2], [3]. Intheirformulation, the procedure for calculatingthe equivalent flat plate areaof a reflector beginswith determination of the intermediatequantities p, q, and r from therelations= sin6sinqS(3.2)q= sin6cos(3.3)r= sin6:inqS(3.4)where a, b, and c are thecorner lengths of the reflector alongthe x, y, and z axes,respectively,and the direction of incidence isgiven by the angles 6 andqas depicted in Figure 3.3.Thevalues given to p, q, and r are thenreassigned in order of increasingmagnitude such thatjpj JqJjr. For a reflector with triangularpanels, the equivalent flatplate area is given byp2 + q2 + r2abcp+q+r—2, p+qT,p+q+r(35)4abc[p++rj, p+qr,while for a reflector withelliptical panels, the area is givenby(M N)2 L2M N2—Labc [tan_1(4L)+tan’( )+vtan_1((M4L2)], L+MN,abc[tan_1(2)+tan’(LN)]L +M N,(3.6)where L=p2, M = q2,and N = r2. For a reflector withrectangular panels, the equivalentflatplate area is simply givenbyabcp(4 — r/q) , q r/2,(3 74abc(pq/r), q r/2.Chapter 3. Truncation and Compensation of Trihedral Corner Reflectors 53These closed-form expressions permit rapid and efficient computation of the equivalent flatplate area of trihedral corner reflectors with panels having certain specific shapes. In order toapply Spencer’s model to trihedral corner reflectors with panels having more general shapesKeen [14), [15) devised a numerical method for determining the size and shape of the polygondefined by the intersection of the entrance and exit pupils.Although Spencer’s model accurately predicts the equivalent flat plate area of trihedral corner reflectors with many different panel geometries, it may fail without indication when appliedto reflectors with panels of completely arbitrary shape. This was first noticed by Robertson [4]who proposed an alternative geometric model which will always yield the correct solution.Robertson’s model is based on the observation that the absolute values of the coordinates atwhich a ray incident from a given direction intersects the three planes defined by the trihedralaxes are identical to the coordinates which define the points of reflection of the ray when it isincident upon one of the interior surfaces of a trihedral corner reflector. A physical implementation of Robertson’s model is presented in Figure 3.5. Consider the trihedral corner reflectorwith panels of arbitrary shape which is shown in Figure 3.5(a). First, the panels of the reflector are replaced by complementary apertures which are derived from each reflecting panel byreflection about the trihedral axes as shown in Figure 3.5(b). The optical model which resultsis shown in Figure 3.5(c). To the observer, the polygons defined by the three complementaryapertures are projected onto a view plane which is normal to the direction of incidence. Thearea common to all three polygons is the equivalent flat plate area of the reflector.Examples of trihedral corner reflectors for which Spencer’s model gives accurate predictionsof the equivalent flat plate area are shown in Figure 3.6 while examples for which the predictionsare erroneous are shown in Figure 3.7. From these cases, it appears that Spencer’s model willgive the correct result if the projections of all three complementary apertures in Robertson’smodel are convex polygons. However, if one or more of the aperture polygons is concave, it ispossible that Spencer’s model will fail without indication and give an incorrect result which willbe larger than the actual value. Although a more formal study of the conditions under whichSpencer’s model fails was not pursued, it is clear that Keen’s numerical implementation is notsuitable for use with trihedral corner reflectors with panels of completely arbitrary shape.Chapter 3. Truncation and Compensation of Trihedral Corner Reflectors54zxz(a)(c)Figure 3.5: A physical implementation of Robertson’s model for the contribution oftriple-bounce reflections to the response of a trihedral corner reflector. (a) A trihedral corner reflector with panels of arbitrary shape. (b) Aperture planes derived from each panel byreflection about the trihedral axes. (c) An optical model for the equivalent flat plate area of atrihedral corner reflector. (from [16,p.240])(b)Chapter 3. Truncation and Compensation of Trihedral Corner Reflectors55zIprojection ofy-z aperturepolygonprqection ofx-y aperturepolygonexit pupil/prection ofz-x aperturepolygonx_____yzIxprojection ofprojection ofentrance pupilz-x aperturex-y aperturepolygonpolygon(a) (b) (c)Figure 3.6: Examples of trihedral corner reflectorsfor which Spencer’s model gives accuratepredictions of the equivalent flat plate area. (a) Reflectorgeometry. (b) Spencer’s model:projection of the entrance and exit pupils of the reflectoronto the view plane and determinationof their intersection. (c) Robertson’s model: projectionof the aperture planes onto the viewplane and determination of their intersection.Chapter 3. Truncation and Compensationof Trihedral Corner Reflectors 56(b)Figure 3.7: Examples of trihedral corner reflectors for which Spencer’s model giveserroneouspredictions of the equivalent flat plate area. (a) Reflector geometry. (b) Spencer’smodel:projection of the entrance and exit pupils of the reflector onto the view plane anddeterminationof their intersection. (c) Robertson’s model: projection of the aperture planes onto theviewplane and determination of their intersection.zentrance pupPxexit pupilprojection ofx-y aperturepolygonprojection ofz-x aperturepolygonzIexit pupil(xNy(a)projection ofz-x aperturepolygonprojection ofx-y aperturepolygon(c)Chapter 3. Truncation and Compensationof Trihedral Corner Reflectors 57An algorithm for solving Robertson’smodel for the contribution of triple-bouncereflectionsto the equivalent flat plate area of an ideal trihedralcorner reflector can be devised using anapproach similar to the one followed by Keen[14], [15] in his solution of Spencer’s model.The four geometric primitives usedin this algorithm are defined as follows [17]—[19j: Apointis specified by its coordinates P(x, y, z). A line segmentis specified by giving its end pointsPi(xi,yi,zi) and P2(x,Y2,z2). A polyline is a chain of connected line segmentswhich isspecified by giving a list of the vertices F1,. . .,F that define the line segments. The firstvertex is called the initial or starting point whilethe last vertex is called the final or terminalpoint. A polygon is a closed polyline inwhich the initial and terminal points coincide. The linesegments P1 F2,F2 F3,. . . , FN P1 are called the edges ofthe polygon. The vertex list for theexterior boundary of the polygon istraversed in a counterclockwise direction and theenclosedregion has a positive vector area. Ifthe polygon contains interior boundaries (or holes), thecorresponding vertex lists are traversed in a clockwisedirection and the enclosed regions havea negative vector area. Once the polygons which represent thepanels of the reflector and thedirection of incidence have been specified, the predictionalgorithm is executed in four steps:1. The polygons which represent thex-y, y-z, and z-x reflecting panels are converted intocorresponding aperture polygons by reflectionabout the principal axes of the trihedral.2. The x-y, y-z, and z-x aperturepolygons are projected onto a view plane which containsthe origin and is normal to the direction of incidence.3. The polygon which representsthe region that is common to the projection of allthreeaperture polygons is determined. This is accomplishedby calculating the intersection ofthe projection of the x-y aperture polygonand the projection of the y-z aperture polygonthen calculating the intersection of theresult and the projection of the z-x polygon.4. The area of the polygon which representsthe region that is common to the projection ofall three aperture polygons is calculated.Chapter 3. Truncation and Compensation of Trihedral Corner Reflectors 58The first step in the prediction algorithm, conversion of the polygons which represent thereflecting panels into aperture polygons by reflection about the principal axes of the trihedral,can be performed by inspection. The second step, projection of the aperture polygons ontoa view plane which contains the origin and is normal to the direction of incidence, may beaccomplished by a transformation of coordinates through pure rotation. If the reflector frameis defined by x,y,and z axes of the trihedral corner reflector, let the x1y” plane define theview plane and let the z’ axis be coincident with the direction of propagation of the reflectedwave, as suggested by Figure 3.8. If the direction of propagation is given by the elevation andazimuth angles 8 and4,the view plane is defined by—++--=O, (3.8)and the vector ‘ is given by(3.9)where the direction cosines a, 3, and y are given bya = sin8sinq, (3.10)/3 = sinOcos, (3.11)= cosO. (3.12)zx,/4// Z// ——/ —/ ——/ ——/ ——/—Figure 3.8: Transformation of coordinates by pure rotation.Chapter 3. Truncation and Compensation ofTrihedral Corner Reflectors 59From (3.9), it can be shown that the view plane and reflectorcoordinate systems are relatedby a transformation matrix of the form4 4 4X £11 £12 £13 X= t21 t22 t23, (3.13)cosc cosf3 cos zwhere t3 corresponds to elements of the coordinate transformation matrix withunknown values.Since the angle of rotation of the x’-y’ plane with respect to the z’ axescan be defined arbitrarily,it is convenient to set 113 to zero so that a vector parallel to the z axiswill have only acomponent in the view plane. Then, from the unitary property of thetransformation matrix,it can be shown that 123 = SlIt7.Thus,29 t11 t12 0 x= t21 122 51117. (3.14)z’ cosa cos/3 cos zIn order to determine the values of the remaining elements,the coordinate transformationmatrix may be compared to a prototype transformationmatrix which corresponds to rotationabout the z-axis through angle ‘çb followed by rotation about thex’-axis through angle ç whichgives29 cosib sin& 0 x= —cosCsint’ —cosCcos& sin( y. (3.15)z’ sin c sin — sin (cos cos(zBy inspection of (3.14) and (3.15), it can be shown thatsin( sin,(3.16)cos( E COS7 , (3.17)sin’çb cosa/sin( , (3.18)cos& E —cos/3/sin(. (3.19)Chapter 3. Truncation and Compensationof Trihedral Corner Reflectors 60From (3.16)—(3.19), it can be shownthat the view plane coordinate frame is related to thereflector coordinate frame by—cos/3/sin7 cosc/sin’y 0 x= —coscr/tan7 —cos/3/tan7 sin7 y. (3.20)cosc cos,6 cos zwhere , 8, and7are given by (3.1O)—(3.12). The projection of each of theaperture polygonsonto the view plane can be determined by transforming the coordinates of eachvertex from thereflector frame into the view plane frame using (3.20) and setting theirz’ coordinates to zero.The third step in the prediction algorithm, determiningthe region of the view plane whichis common to the projection of all three aperture polygons, is moredifficult. A variety ofalgorithms for determining the intersection of overlappingpolygons have been developed for usein computer graphics applications and are widely used. Most of these,including the SutherlandHodgman and Liang-Barskey polygon-clipping algorithms, are unsuitablefor use in predictionalgorithm derived from Robertson’s model because they requireat least one of the polygons tobe convex. However, the Weiler-Atherton polygon-clipping algorithmovercomes this limitationand is capable of clipping a concave polygon with interiorholes to the boundaries of anotherconcave polygon with interior holes [18]—[21j.In the Weiler-Atherton polygon-clipping algorithm,the subject and clip polygons are described by circular lists of vertices S1,52,. . . , SMand C1,C2,. . . , CN, respectively. Before theactual clipping is performed, the points at which the subject and clippolygons intersect aredetermined. The coordinates of the intersection points areinserted into both the subject andclip polygon vertex lists in the appropriate sequence. In order to establisha bidirectional linkbetween the vertex lists, each intersection point in thesubject polygon vertex list is given apointer to the location of the same intersection point in the clip polygonlist and vice versa. Theactual clipping is performed as follows: The subject polygon is traversedin a counterclockwisedirection until an intersection is reached. If thisseries of points lies in the interior of the clippolygon, they are added to the result list. If thenext vertex of the subject polygon lies insideChapter 3. Truncation and Compensation ofTrihedral Corner Reflectors61the clip polygon, the subject polygonvertex list is followed. Otherwise, thealgorithm jumps tothe clip polygon vertex list and follows it to thenext intersection. This process continues untilall the intersections have beentraversed and the algorithm has returned to thefirst point inthe result polygon. Methods for implementingthe algorithm and enhancing its efficiency androbustness have beendiscussed in the literature [18)—[21].Application of the Weiler-Atherton polygon-clippingalgorithm to the problem of determining the equivalent flat plate area of a trihedralcorner reflector with triangular panels of equalcorner length for incidence along the symmetryaxis is demonstrated in Figure 3.9. Theprojections of the x-y, y-z, and z-x aperturepolygons onto the view plane are shown in Figure3.9(a).The shaded region represents the areawhich is common to the projections of all three aperturepolygons. Symbolically, this area maybe described byA = ((Z fl X) fl Y), (3.21)where Z, X, and Y are the projections of thex-y, y-z, and z-x aperture polygons, respectively.In Figure 3.9(b), polygon Z is clipped againstpolygon Y to yield a first result. In Figure3.9(c),the first result is clipped against polygon X to yieldthe final result. The polygon which definesthe region common to the projectionof all three aperture polygons is shown in Figure 3.9(d).The last step in the prediction algorithmis calculation of the area of the polygonwhichdefines the region common to theprojection of all three aperture polygons. If thevertices ofthe polygon are given by {S}, the z-directedvector area of the polygon is simply given byA = S x S1. (3.22)Once the equivalent flat plate area hasbeen calculated from (3.22), the scatteringcross sectionof the reflector can be determined from(3.1).Chapter 3. Truncation and Compensation of Trihedral Corner Reflectors62projection ofz-x aperturepolygon_____ ,j \\___________1::/\ / \ ,,/\- VI’ —projection ofy-z aperture(a)polygonZ polygon list X polygon listS1/1/(C) C1/1/(S)52N12C212C3/I/(S)S3/1/(C)5421414 C4Figure 3.9: Execution of the Weiler-Atherton polygon clippingalgorithm. (a) The x-y, y-z, andz-x apertures are projected onto the view plane toyield polygons Z, X, and Y, respectively.(b) Polygon Z is clipped against polygon X to yield the first resultlist. (c) The polygon definedby the first result list is clipped against polygon Y toyield the final result list. (d) The finalresult list defines the polygon A which containsthe area common to polygons Z, X, and Y. Ineach case, S refers to points in the subject polygon,C refers to points in the clip polygon,and I, refers to the points at which the polygons intersect.C4projection ofx-y aperturepolygon14polygon V(b)S5S4polygon ZS3/I3/C3polygon Xfinal result(d)Second Clipping Sequence12(c)First Clipping SequenceFirst result list Ypolygon listSi ClIl C252 1212 1353 C313 C4S4 1414 IChapter 3. Truncation and Compensation of Trihedral Corner Reflectors633.3 Reflectors with Three-Fold SymmetryIn this section, the extent to which the response characteristics of trihedralcorner reflectorswhich present three-fold symmetry can be altered by appropriateshaping of their reflectingpanels is considered. In each case, the contribution of triple-bounce reflectionsto the scatteringcross section of the reflector was calculated over the entire quadrantdefined by the axes of thetrihedral. The resulting array of values in 8 and were then convertedto contours expressedin decibels with respect to the maximum response of thereflector. Although the reflectorcoordinate frame of Figure 3.10(a) is a convenientchoice for analyzing the response of a trihedralcorner reflector, it is not the most natural choice for presentingresponse patterns. Instead, theresults were transformed to the global coordinate frameof Figure 3.10(b) in which the z’ axisis aligned with the local vertical and the direction of maximumresponse (or the boresight ofthe reflector) lies in the horizontal plane and is aligned with thex’ axis.boresightboresightz,(b)Figure 3.10: (a) Reflector coordinate frame and (b) globalcoordinate frame.z(a)z,x,Chapter 3. Truncationand Compensation of Trihedral CornerReflectors64The boresight of a trihedralcorner reflector which presentsthree-fold symmetry isusuallycoincident with the symmetry axis definedby 8 54.74° and4)= 45°. Using the method ofevaluating coordinate transformationmatrices that was derived in section2.4, it can be shownthat the global coordinate frameis related to the reflector coordinateframe by1 1 1X— 7x— 1 1 1 (323— 7 77•/.0 ZThus, direction expressed interms of the angles 9’ and4)’in the global coordinateframe isrelated to direction expressed interms of the angles 0 and4)in the reflector coordinateframeby0’ = cos sin8 cos4)+Vhicosej , (3.24)—3=sin8cos4)+—=sin8sin4)— -3=cos04)’= tanV3 V2V6(3.25)*smn0cos4)_ *slnosln4)_*cos8while the reverse transformation isgiven by0 = cos1[_*sino’cos4)’— çsin0’sin4)’+Jjcoso’] , (3.26)—+sin9’cos4)’+ ±sin0’sin4)’4)= tan1v2 v2, (3.27)sin 8’ cos4)’+sin 8’ sin4)’+cos 9’Although the response contourscould be plotted on a conventionalrectangular grid, distortionof the pattern can be minimizedby using a grid derived from anequal area projection of thetype used in geodesy and cartography[22], [23]. Several common projectionsare suitable, butthe sinusoidal projection is theeasiest to implement and wasselected for use here. Whilethe horizontal axis ofthe response patterns is expressed in termsof the azimuth angle4)’,thevertical axis is expressed in termsof the co-elevation or altitude angle‘ given by8’ = 90 — 9(3.28)in order to simplify interpretationof the results.Chapter 3. Truncation and Compensation of Trihedral CornerReflectors 65The angular coverage of symmetrical trihedral corner reflectorswith triangular, circular,and square panels were predicted using the algorithm described inthe previous section. Thecircular panels were represented by a twenty-sided polygon ofequivalent area. The results arepresented in Figures 3.11, 3.12, and 3.13 and summarized in Table 3.1. It is generallyfoundthat attempts to increase the angular coverage of a symmetrical trihedral cornerreflector witha fixed corner length by modifying the shape of its reflecting panelsare generally accompaniedby a corresponding decrease in the ratio of the scattering cross sectionto the physical size ofthe reflector, as previously noted by Robertson [5].Robertson [4] proposed a method for altering the angular coverage of a symmetricalreflectorwith triangular panels by removing notches of prescribed width aild depth from the outsideedgeof each panel, as shown in Figure 3.14. The depth of the notch canbe expressed by a parameterd/d which is the ratio of the depth of the notch to the length of the medianto the outside edgeof the panel. The angular coverage of a Robertson reflector withshallow notches (d/d = 0.25)is depicted in Figure 3.15. After notching the panels, three equalmaxima appear at about12 degrees off the symmetry axis. Although the angularcoverage of the reflectors increasesdramatically, the maximum response is correspondingly smaller.The angular coverage of areflector with notches of intermediate depth (d/d = 0.50) is depictedin Figure 3.16. The threemaxima are more pronounced and have shifted outward to about 21 degrees off thesymmetryaxis. The angular coverage of a reflector with deep notches (d/d = 0.75) isdepicted inFigure 3.17. The null in the response which has formed along the symmetryaxis is -25 dB withrespect to the maximum response. The results are summarized in Table3.1.The cumulative probability distribution of the response of symmetricaltrihedral corner reflectors with triangular, circular, and square panels, and Robertson panelswith shallow notches(where d/d = 0.25) over the quadrant defined by the axes of the trihedral is shown in Figure3.18. Of the four targets, Robertson’s reflector withshallow notches presents the most uniformresponse. Siegel et al. [2] claimed that the cumulative distribution of atrihedral corner reflectorwith circular panels exceeds that of reflectors with either triangular or square panels.In fact,the results presented here show that its cumulative distribution is onlyintermediate.Chapter 3. Truncation and Compensationof Trihedral Corner Reflectors6660453015C)ci)v 0-15-30-45-60dB-60 -30 0 30 60s’(deg)Figure 3.11: Angular coverage of a trihedral corner reflector withtriangular panels.60453015C)ci).0-15-30-45-60-60 -30 0 30 600’(deg)Figure 3.12: Angular coverage of a trihedral corner reflector with circular panels.Chapter 3. Truncation and Compensation of Trihedral Corner Reflectors 67jBFigure 3.14: Boresight view and panel geometry of trihedral corner reflectors with Robertsonpanels where the ratio of notch depth to panel depth, d/d, is (a) 0.25, (b) 0.50, and (c) 0.75.The shaded portion indicates the equivalent flat plate area of the reflector for incidence alongthe boresight.ci604530150-15-30-45-60-60 -30 0 30 600’(deg)Figure 3.13: Angular coverage of a trihedral corner reflector with square panels.BORESIG HTVIEWPANELGEOMETRYd(a)(b) (c)Chapter 3. Truncation and Compensation of Trihedral Corner Reflectors 68I L I SFigure 3.15: Angular coverage of a trihedral corner reflector with Robertson panels whered/d = 0.25.60453015G)V 0-15-3045-60qS’(deg)a)-D604530150-15-3045-60-60 -30 0 30 60çS’ (deg)-60 -30 0 30 60Figure 3.16: Angular coverage of a trihedral corner reflector with Robertson panels whered/d = 0.50.Chapter 3. Truncation and Compensation of Trihedral Corner Reflectors 69Figure 3.18: Cumulative probability distribution of the response of trihedral corner reflectorswith triangular, circular, square, and Robertson panels (where d/d = 0.25) over the quadrantdefined by the axes of the trihedral.604530150)a)00-15-30-45-60‘ (deg)coverage of a trihedral corner reflector with Robertson panels where-60 -30 0 30 60Figure 3.17: Angulard/d2 = 0.75.1.00.810*6O.4EC)0.20.0-10 -8 -6 -4 -2 0O/L (dB)maxChapter 3. Truncation and Compensation ofTrihedral Corner Reflectors70Table 3.1:Response Characteristics of SelectedTrihedral Corner Reflectors with Three-foldSymmetryMaximum ResponseFigure Panel Geometry Umax(1)Umax(2)3.11 Triangular (4ir/3)a4/A20.0 dB3.12 Circular 15.6a/5.7dB3.13 Square12ira4/)2 9.5 dB3.15 Robertson) 0.39a4/A2 -4.1 dB3.16 Robertson) 0.18a4/)i2 -7.0 dB3.17 Robertson(s)0.12a4/)2 -9.4 dB1 and 3 dB Elevation and AzimuthalBeamwidths of the Main Response LobeFigure Panel Geometry0iB 41dB 03dB 3dB3.11 Triangular 24°24° 39° 39°3.12 Circular 18° 17°31° 30°3.13 Square 8° 8°22° 20°3.15 Robertson(s) 41° 40°52° 50°3.16 Robertson() n/a n/an/a n/a3.17 Robertson)n/a n/a n/a n/a6 and 10 dB Elevation andAzimuthal Beamwidths of the Main ResponseLobeFigure Panel Geometry O6 0o 41odB3.11 Triangular 52° 51°63° 61°3.12 Circular44° 43° 57° 55°3.13 Square 36° 35°50° 50°3.15 Robertson(s)61° 59° 70° 66°3.16 Robertson) 65°64° 73° 68°3.17 Robertson(5)n/a n/a n/a n/aNotes: (1) where a is the corner length ofthe reflector.(2) relative to Umax of a trihedral cornerreflector with triangularpanels of the same corner length.(3) where the ratio of notchdepth to panel depth, d/d, is 0.25.(4) where d/d = 0.50.(5) where d/d = 0.75.Chapter 3. Truncation and Compensation of Trihedral Corner Reflectors 713.4 Reflectors with Bilateral SymmetryThe response characteristics of bilaterally symmetric trihedral corner reflectors with triangular,elliptical, and rectangular panels can be predicted using the closed-form expressions for theirequivalent flat plate area which are given by (3.5), (3.6), and (3.7), respectively. Although theform of the expressions makes it difficult to derive closed-form solutions for either the size anddirection of the maximum response or the azimuthal and elevation beamwidths of the mainresponse lobe, it is a relatively simple matter to determine these quantities using a numericalapproach. The problem geometry and the reflector coordinate system are shown in Figure 3.3.It is convenient to define bilateral symmetry by a mirror plane which contains the z axis andbisects the x-y plane at an azimuthal anglepf45 degrees. If the length of the corner along thez axis is given by c and the length of the corners along the x and y axes are given by a, thenthe reflector aspect ratio is defined as c/a.If a reflector presents a single main response lobe and is bilaterally symmetric, the directionof its maximum response must lie in its mirror plane. While the azimuthal anglecmaxof thedirection of maximum response is a fixed quantity, the corresponding elevation angle6mis afunction of the reflector aspect ratio c/a. A golden section search [24] was used to determine theelevation angle of maximum response of bilaterally symmetric trihedral corner reflectors withtriangular, elliptical, and rectangular panels. The results are presented in Figure 3.19. Althoughthe directions of maximum response are coincident when the reflector aspect ratio is unity,they diverge as the reflector becomes increasingly asymmetric. Trihedral corner reflectors withtriangular panels are of special interest since they present a planar aperture which facilitatesthe attachment of either a transmission polarizer or a protective dielectric cover. Since theintercepts of the plane which defines the reflector aperture are given by the corner lengths ofthe reflector, it can be shown that the elevation angle of the normal to the reflector aperture isgiven by= cos1(2a+C2). (3.29)Chapter 3. Truncation and Compensationof Trihedral Corner Reflectors 72The angle of maximum response8mand the angle of the normal aperture 8 of bilaterallysymmetric trihedral corner reflectors with triangular panels areplotted as a function of thereflector aspect ratio c/a in Figure 3.20. Although the direction of maximumresponse and thenormal to the reflector aperture are coincident if the reflectoraspect ratio is unity, they divergeas the reflector becomes increasing asymmetric. If thereflector aspect ratio is less than unity,the normal to the aperture is lower than the direction of maximum response.If the ratio isgreater than unity, the reverse is true.If a trihedral corner reflector has cornersof equal length £, it is convenient to express itsscattering cross section in the formumax— K. (3.30)The parameter K is a figure of merit which can be used to comparethe maximum responseof trihedral corner reflectors with similar dimensions but differentpanel shapes. Although itis desirable to employ a similar scheme to express the scattering crosssection of a bilaterallysymmetric reflector, the equivalent corner length of’ such a reflector must bedefined. Althoughthe equivalent corner length could be defined in severaldifferent ‘ways, it is convenient to simplytake the arithmetic mean of the three corner lengths, i.e.,a+2c(3.31)Using this definition, the parameter K for bilaterally symmetrictrihedral corner reflectors withtriangular, elliptical, and rectangular panels is plotted as a function of the reflectoraspect ratioin Figure 3.21.The azimuthal and elevation beamwidths of bilaterallysymmetric trihedral corner reflectorswith triangular, elliptical, and rectangular panels were determined by applying a bracketingandbisection algorithm [24] to (3.5), (3.6), and (3.7) in a global coordinate framesimilar to thatdefined in Figure 3.10(b) where the z’ axis is aligned with the local vertical andthe x’ axisis aligned with the direction of maximum response. Theresults are presented as a functionof the reflector aspect ratio c/a in Figures 3.22, 3.23, and 3.24,respectively. The increase inChapter 3. Truncation and Compensation of Trihedral Corner Reflectors 73and the converse for large reflector aspect ratios is evident. The angular coverage of bilaterallysymmetric trihedral corner reflectors with triangular panels and reflector aspect ratios of 4.0and 0.25 are presented in Figures 3.25 and 3.26, respectively.In certain cases, it has been found useful to realize bilaterally symmetric trihedral cornerreflectors which are composed of combinations of triangular, elliptical, and rectangular panels.For example, the trihedral corner reflector developed by the European Space Agency for useas a calibration target in the SAR-580 program is composed of triangular side panels and asquare center panel [8]. The angular coverage of the SAR-580 calibration target and five otherreflectors which have been realized in a similar manner are presented in Figures 3.27 through3.32. The results are summarized in Table 3.2.During the course of this study, an alternative approach to the design of bilaterally symmetric trihedral corner reflectors for use in radar-assisted positioning systems was suggested byHelmut Lanziiier of Offshore Systems Ltd. [25]. The primary design objective was to produce aset of reflectors which would support themselves on a horizontal surface with their direction ofmaximum response in the horizontal plane. A secondary objective was to provide most of thereflectors with a planar aperture in order to facilitate the attachment of a protective cover or atransmission polarizer. Three variations were devised as shown in Figure 3.33. The truncatedreflector of Figure 3.33(a) is simply a symmetrical trihedral corner reflector with triangular panels which has been inverted and had its side panels truncated iii such a way that structure is selfsupporting when placed on a horizontal surface but the effective flat plate area for incidencealong the boresight is unaffected. The compensated reflector of Figure 3.33(b) is developedfrom the truncated reflector by increasing the size of the side panels until maximum apertureefficiency is obtained. The extended reflector of Figure 3.33(c) is developed from the compensated reflector by extending the side panels further still. While the aperture of the reflectoris no longer planar, the result is a large increase in the azimuthal beamwidth. The angularcoverage of the truncated, compensated, and extended reflectors is presented in Figures 3.34,3.35, and 3.36, respectively. The results are summarized in Table 3.2.Chapter 3. Truncation and Compensationof Trihedral Corner Reflectors749075604530150I0.1 10Reflector aspect ratio - c/aFigure 3.19: Angle of maximum response6mof bilaterally symmetrictrihedral corner reflectors with triangular, elliptical, and rectangularpanels vs. the reflector aspectratio c/a.90756045301500.1 10Reflector aspect ratio - c/aFigure 3.20: Angle of maximum response °mand angle of the normal to the aperture 8 ofbilaterally symmetric trihedral corner reflectorswith triangular panels vs. the reflectoraspectratio c/a.1Chapter 3. Truncation and Compensation of Trihedral Corner Reflectors 75Figure 3.21: Maximum response of a bilaterally symmetric trihedralcorner reflectors with triangular, elliptical, and rectangular panels vs. the reflector aspect ratio c/awhere = K £ /A2and £ = (a+2c)/3.Reflector aspect ratio - c/aFigure 3.22: Azimuthal and elevation beamwidths of a bilaterally symmetric trihedral cornerreflector with triangular panels vs. the reflector aspect ratio c/a.2010I10010I0.100.1 110Reflector aspect ratio - c/a-10907560.45E301500.1 I 10Chapter 3. Truncation and Compensation of TrihedralCorner Reflectors 76Reflector aspect ratio - c/aFigure 3.23: Azimuthal and elevation beamwidthsof a bilaterally symmetric trihedral cornerreflector with elliptical panels vs. thereflector aspect ratio c/a.Reflector aspect ratio - c/aFigure 3.24: Azimuthal and elevation beamwidths of a bilaterally symmetric trihedral cornerreflector with rectangular panels vs. the reflector aspect ratio c/a.90750.1 1 10907545301500.1 1 10Chapter 3. Truncation and Compensationof Trihedral Corner Reflectors 7760453015D)a)-15-30-45-60qS’ (deg)Figure 3.25: Angular coverage of a bilaterally symmetric trihedraicorner reflector with triangular panels and reflector aspect ratio c/a = 0.25.A..N-60 -30 0 30 60‘ (deg)Figure 3.26: Angular coverage of a bilaterally symmetric trihedralcorner reflector with triangular panels and reflector aspect ratio c/a = 4.0.-60 -30 0 30 6060453015C). 0I-15-30-45-60Chapter 3. Truncation and Compensation of Trihedral Corner Reflectors 7860453015C)a)-15-30-45-60-60 -30 0 30 60‘ (deg)Figure 3.27: Angular coverage of a bilaterally symmetric trihedral corner reflector with triangular side panels and a circular center panel.60453015C)a)DOI-15-30-45-60-60 -30 0 30 60q5’(deg)Figure 3.28: Angular coverage of a bilaterally symmetric trihedral corner reflector with triangular side panels and a square center panel.Chapter 3. Truncation and Compensation of Trihedral Corner Reflectors 7960453015C)a)DO-15z-60-3003060‘ (deg)Figure 3.29: Angular coverage of a bilaterally symmetric trihedral corner reflector with circularside panels and a triangular center panel.60453015C)DO-15-30-45-60qS’ (deg)Figure 3.30: Angular coverage of a bilaterally symmetric trihedral corner reflector with circularside panels and a square center panel.I /VtE.N\-30-45-60-60 -30 0 30 60Chapter 3. Truncation and Compensation of Trihedral Corner Reflectors 801Figure 3.31: Angular coverage of a bilaterally symmetric trihedral corner reflector with squareside panels and a triangular center panel.604530Figure 3.32: Angular coverage of a bilaterally symmetric trihedral corner reflector with squareside panels and a circular center panel.60453015C)-15-30-45-60V.-60 -30 0 30 60‘ (deg)0)a,•0150-15-30-45-60I ••.i / 1 ..• I..••-60 -30 0 30 600’(deg)Chapter 3. Truncation and Compensation of Trihedral Corner Reflectors81flyzflBORESIGHTVIEW22a,‘ a--\‘ /‘ /\ /“I(a) (b) (C)SIDE PANELGEOMETRYFigure 3.33: Evolution of Lanziner’s bilaterally symmetric trihedral corner reflector through(a) truncation, (b) compensation, and (c) extension of the triangular side panels.The shadedportion indicates the equivalent flat plate area of the reflector for incidence along theboresight.604530Figure 3.34: Angular coverage of a bilaterally symmetrictrihedral corner reflector with truncated triangular side panels and a triangular center panel.i115-15-30-45-60-60 -30 0 30 600’(deg)Chapter 3. Truncation and Compensation of Trihedral Corner Reflectors 82r1604530._ 150)a)0‘—0I-15-30-45-60....1 •..•‘ 1•.•’ I\i//-s’-..\(///1d-60 -30 0 30 60ç5’ (deg)Figure 3.35: Angular coverage of a bilaterally symmetric trihedral corner reflector with truncated and compensated triangular side panels and a triangular center panel.flifl604530150)a)--15-30-45-60•...1 ..-60 -30 0 30 60‘ (deg)Figure 3.36: Angular coverage of a bilaterally symmetric trihedral corner reflectorwith truncated, compensated, and extended side panels and a triangular center panel.Chapter 3. Truncation and Compensation of Trihedral Corner Reflectors 83Table 3.2:Response Characteristics of Selected Trihedral Corner Reflectors with Bilateral SymmetryMaximum ResponseFigure Side Panels Center Panel(1)max(2)9m3.27 Triangular Circular 6.4a4/A2 1.8 dB 58.5°3.28 Triangular Square 8.7a/A 3.2 dB 61.4°3.29 Circular Triangular 10.0a4/A2 3.8 dB 51.6°3.30 Circular Square 21.3a4/A2 7.1 dB 56.9°3.31 Square Triangular 18.8a4/A2 6.5 dB 49.8°3.32 Square Circular 28.6a4/A2 8.3 dB 52.8°3.34 Truncated Triangular (4K/3)a4/A2 0.0 dB 54.7°3.35 Compensated Triangular (167r/3)a4/A2 6.0 dB 54.7°3.36 Extended Triangular (l6ir/3)a4/A2 6.0 dB 54.7°1 and 3 dB Elevation and Azimuthal Beamwidths of the Main Response LobeFigure Side Panels Center PanelOi 4’ii 03dB 3d]33.27 Triangular Circular 21° 22° 35° 37°3.28 Triangular Square 17° 17° 31° 31°3.29 Circular Triangular 21° 19° 35° 32°3.30 Circular Square 13° 14° 27° 26°3.31 Square Triangular 14° 11° 28° 25°3.32 Square Circular 11° 10° 25° 24°3.34 Truncated Triangular 16° 23° 30° 38°3.35 Compensated Triangular 8° 8° 20° 24°3.36 Extended Triangular 8° 29° 20° 41°6 and 10 dB Elevation and Azimuthal Beamwidths of the Main Response LobeFigure Side Panels Center PanelO6 4’6dB OiodB 4’lOdB3.27 Triangular Circular 48° 49° 60° 60°3.28 Triangular Square 45° 45° 58° 57°3.29 Circular Triangular 48° 45° 60° 57°3.30 Circular Square 41° 39° 55° 53°3.31 Square Triangular 41° 40° 54° 53°3.32 Square Circular 38° 39° 52° 52°3.34 Truncated Triangular 43° 51° 56° 61°3.35 Compensated Triangular 33° 40° 45° 55°3.36 Extended Triangular 33° 52° 45° 62°Notes: (1) where a is the corner length of the center panel of the reflector.(2) relative toUmof a trihedral corner reflector with triangularpanels of the same corner length.Chapter 3. Truncation and Compensation of Trihedral Corner Reflectors 843.5 Effect of Errors in Construction on Reflector PerformanceA trihedral corner reflector will present an optimum response if its three reflecting panels areperfectly flat and its corner angles are exactly 90 degrees. Approximate methods for determiningthe extent to which the response degrades as the reflecting panels deviate from perfect flatnessand mutual orthogonality have been presented by Spencer [1], Keen [26], and Trebits [27]. Thereduction in the scattering cross section of a symmetrical trihedral corner reflector with triangular panels due to errors in all three corner angles is presented in Figure 3.37. In general, thetolerances on the corner angles generally decrease as the size of the reflector increases. It hasalso been found that the tolerances on the corner angles depend on the shape of the reflectingpanels and the size of the equivalent flat plate area. Although the approximate methods described by Spencer, Keen, and Trebits may be used to determine the effect of panel deviationson the response of trihedral corner reflectors with modified panel geometries, numerical techniques such as the finite-difference time-domain (FD-TD) and the shooting and bouncing ray(SBR) methods could also be used to perform the necessary calculations if sufficient computational resources are available [10], [11]. Alternatively, experimental techniques may be used toassess the effect of construction errors on the response of trihedral corner reflectors [28].0wFigure 3.37: Effect of errors in all three corner angles on the response of a trihedral cornerreflector with triangular panels for incidence along the the symmetry axis. (from [5,p.13-12])40w3300z 20-Jw>w-J89 90 91ANGLEBIN DEGREESChapter 3. Truncation and Compensation of Trihedral CornerReflectors 853.6 ConclusionsMethods for predicting the scattering cross section and angular coverage of a conventional trihedral corner reflector with panels of completely arbitrary shapehave been considered. Althoughnumerical techniques such as the finite-difference time-domain (FD-TD)and the shooting andbouncing ray (SBR) methods can predict both the contributionsof single, double, and triple-bounce reflections from the interior and the effect ofdeviations of the reflecting panels fromperfect flatness or mutual orthogonality, they are extremely demandingcomputationally. If thereflector is sufficiently large, it is usually sufficient to account for the contribution oftriple-bounce reflections from the interior since they completely dominate the response for most directions of incidence. Spencer’s model [1] for the contribution of triple-bounce reflections to theequivalent flat plate area of an ideal trihedral corner reflector with either triangular or squarepanels, and the prediction algorithm subsequently derived from it by Keen [14], [15], mayfailwhen applied to reflectors with panels of completely arbitrary shape but analternative modelproposed by Robertson [4] will always provide the correct solution. An efficient and robustnumerical method for solving Robertson’s model has been presented.The response patterns of ideal trihedral corner reflectors which present three-foldsymmetryhave been plotted on grids derived from equal area projections and compared. If the cornerlength of the reflector is fixed, it is generally found that attempts to increase thebeamwidthof the response by modifying the shape of the panels are accompanied by a reductionin theamplitude of the maximum response. If three-fold symmetry is broken so thatthe reflectorsimply presents bilateral symmetry about a mirror plane which contains one of the trihedralaxes and bisects the opposite panel, the beamwidth of the response in one principalplanecan be increased relative to the beamwidth in the orthogonal plane by modifyingthe shape ofthe reflecting panels in a suitable manner. The additional degree of freedom mayalso proveuseful when a reflector must be designed subject to a constraint such as arequirement thatthe modified reflector present a planar aperture in order to facilitate the attachment ofeithera transmission polarizer or a protective cover. A set of design curves for bilaterally symmetricChapter 3. Truncation and Compensation ofTrihedral Corner Reflectors 86reflectors which are composed solely oftriangular, elliptical, or rectangular reflectingpanels hasbeen presented. The curves present the elevation angleof the direction of maximumresponse,the amplitude of the maximum response,and the elevation and azimuthal beamwidthsas afunction of the reflector aspect ratio. Forthe case of a bilaterally symmetric trihedralcornerreflector with triangular panels, the elevationangles of the direction of maximum responseandthe normal to the reflector aperturehave been compared as a function of the reflectoraspectratio. The response patterns of selectedbilaterally symmetric reflectors which are composed ofcombinations of panels with various shapesincluding triangular, circular, and square have alsobeen presented and compared. The EuropeanSpace Agency’s SAR-580 calibrationtarget isincluded in this set. A related problem, thedesign of top hat reflectors with specified responsecharacteristics, is considered in AppendixA.Degradation of the response of a trihedralcorner reflector caused by deviation of its reflecting panels from perfect flatness ororthogonality is an important consideration in thedesignand fabrication of calibrationtargets and location markers. Although the approximatemethodspresented by Spencer [1], Keen [26], andTrebits [27] may be used to determine theeffect ofpanel deviations on the response of trihedralcorner reflectors with modified panel geometries,numerical techniques such as the finite-differencetime-domain and the shooting and bouncing ray methods mentioned above can alsobe used to perform such calculations ifsufficientcomputing resources are available.References[1] R. C. Spencer, Optical Theory of the Corner Reflector. Cambridge, MA: MIT Rad.Lab.Tech. Rep. 433, 2 Mar. 1944.[2] K.M. Siegel et al., Studies in Radar Cross Sections XVIII - Airborne PassiveMeasuresand Countermeasures. Ann Arbor, MI: Univ. Michgan, Jan. 1956. (cited in C.G. Bachman,Radar Targets. Lexington, MA: Heath, 1982,p.71—101.)[3] A.L. Maffett, Topics for a Statistical Description of Radar Cross Section. New York:Wiley,1989,pp.190—201.[4] 5. D. Robertson, “Targets for microwave radar navigation,” Bell Syst. Tech. J.,vol. 26,pp.852—869, 1947.[5] W. C. Jakes, Jr., and S. D. Robertson, “Passive Reflectors,” in AntennaEngineeringHandbook, H. Jasik, Ed. New York: McGraw-Hill, 1961, chap. 13.[6] P. 0. Gillard and K. B. Whiting, Ground Plane Corner Reflectors for NavigationandRemote Indication. U.S. Patent No. 4,104,634, dated Aug. 1, 1978.[7] S. Haykin, “Polarimetric radar for accurate navigation,” Can. J. Elect. Comp.Eng., vol. 17,pp.130—135, July 1992.[8] D.G. Corr, A.D. Woode, and S. Bruzzi, “The SAR-580calibration site at RAE, Bedford,”mt. J. Remote Sensing, vol. 3,pp.223-227, 1982.[9] A. Freeman, Y. Shen, and C. L. Werner, “Polarimetric SARcalibration experiment usingactive radar reflectors,” IEEE Trans. Geosci. Remote Sensing, vol. GE-28,pp.224—240,Mar. 1990.87References88[101A. Taflove and K.R. Umashankar, “Review of FD-TDnumerical modeffing of electromagnetic wave scattering and radar cross section,”Proc. IEEE, vol. 77, pp. 682—699, May1989.[11] J. Baldauf, S-W Lee, L. Lin, S-K Jeng, S.M. Scarborough,and C.L. Yu, “High frequencyscattering from trihedral corner reflectors and otherbenchmark targets: SBR versus experiment,” IEEE Trans. Antennas Propagat.,vol. 39,pp.1345—1351, Sept. 1991.[12] A.C. Polycarpou, C.A. Balanis, and P.A. Tirkis,“Radar cross section evaluation of thesquare trihedral corner reflector using P0 and MEC,” in IEEE/AP-S Symp. Dig. (AnnArbor, MI), June 1993,pp.1428—1431.[131Z.0. Al-hekail and W.D. Burnside, “Scattering fromcorner reflectors: a hybrid approach,”in IEEE/A P-S Symp. Dig. (Ann Arbor, MI), June 1993,pp.1432—1435.[14] K. M. Keen, “New technique for the evaluationof the scattering cross-sections of radarcorner reflectors,” lEE Proc. H, vol. 130,pp.322—326, Aug. 1983.[15] K. M. Keen, “Fast algorithm for the exact determination of themapped effective areas oftrihedral radar reflectors,” Electron. Lett., vol. 19,pp.1014—1015, Nov. 24, 1983.[16] D.R. Brown, R.J. Newman, and J.W. Crispin, Jr., “RCSEnhancement Devices,” in Methods of Radar Cross Section Analysis, J.W. Crispin, Jr. and K.M. Siegel,Eds. New York:Academic Press, 1969,pp.237—280.[17) R. Sedgewick, Algorithms, 2nd ed. Reading,MA: Addison-Wesley, 1988,pp.347—356.[18] D.F. Rogers, Procedural Elements for Computer Graphics.New York: McGraw-Hill, 1985,pp.111—188.[19] J.D. Foley, A. van Dam, S.K. Feiner, and J.F.Hughes, Computer Graphics: Principles andPractice, 2nd ed. Reading, MA: Addison-Wesley, 1990,pp.124—127,References89[20] K. Weiler and P. Atherton,“Hidden surface removal using polygonarea sorting,” ComputerGraphics, vol. 11,pp.214—222, Summer 1977.[21] K. Weiler, “Polygoncomparison using a graphrepresentation,” Computer Graphics,vol. 14,pp.10—18, Spring 1980.[22] P. Richardus and R.K. Adler, MapProjections for Geodesists,Cartographers, and Geographers. Amsterdam: North-Holland, 1972.[23] J.P. Snyder and P.M. Voxiand,An Album of Map Projections. U.S. Geol.Surv. Prof. Pap.,no. 1453, 1989.[24] W. H. Press, B. P. Flannery, S.A. Teukolsky, and W. T. Vetterling,Numerical Recipes,2nd ed. Cambridge Univ. Press: Cambridge,1992.[25] D. G. Michelson and H. H. Lanziner,Radar Reflector to Enhance Detection. U.S.PatentNo. 4,990,918, dated Feb. 5, 1991.[26] K. M. Keen, “Predictionof scattering cross-section reductionsdue to plate orthogonalityerrors in trihedral radar reflectors,”Electron. Lett., vol. 19,pp.115—117, Feb. 3, 1983.[27] R.N. Trebits, “Radar Cross Section,”in Radar Reflectivity Measurement:Techniques andApplications, N.C. Currie, Ed. Norwood, MA:Artech House, 1989,pp.46—49.[28] D. Kähny and J. van Zyl,“How does corner reflector construction affectpolarimetric SARcalibration?” in Proc. IGARSS’90 (CollegePark, MD), May 1990,pp.1093—1096.Chapter 4DEPOLARIZING TRIHEDRAL CORNER REFLECTORS4.1 IntroductionA conventional trihedral corner reflector returns linearly polarized incident waves without modification but reverses the sense of elliptically or circularly polarized waves. This is often referredto as a regular polarization response since it is also characteristic of spheres and flat plates.Although this is an ideal response for targets intended for use with conventional radars whichemploy either horizontally or vertically polarized antennas for both transmission and reception,radars which employ same-sense circular polarization to assist in rain clutter suppression or,more recently, various forms of polarization diversity to assist in target classification and identification often require calibration targets and location markers with twist-polarizing, circularlypolarizing, or linear polarization selective responses [1], [2]. Several methods for altering thepolarization response of conventional trihedral corner reflectors have been developed over theyears [3]—[12]. The three basic approaches are shown in Figure 4.1.Figure 4.1: Methods for altering the polarization response of a conventional trihedral cornerreflector. (a) Removal of one reflecting panel and loading of the interior of the reflector witha low-loss dielectric material (2.3 Er 4). (b) Installation of a transmission polarizer acrossthe reflector aperture. (c) Replacement of one of the reflecting panels by a reflection polarizer.(a) (b)(c)90Chapter 4. Depolarizing Trihedral Corner Reflectors 91Since it is often necessary to deploy large numbers of calibration targets and location markersin the field for extended periods of time where they are exposed the effects of sun, wind, andrain, the mechanical ruggedness of trihedral corner reflectors which incorporate depolarizingelements and the ease with which such reflectors can be manufactured are important considerations. Schemes which require that a transmission polarizer be installed across the reflectoraperture become increasingly difficult to implement as the reflector grows larger and cannot beapplied at all if the reflector doesn’t present a planar aperture [4]—[6]. Schemes which involvereplacing or augmenting one of the reflecting panels with a reflection polarizer oriented suchthat the axis of the grating is parallel to one of the principal axes of the trihedral have attractedconsiderable interest, particularly in recent years, because they avoid many of the mechanicalproblems associated with the use of transmission polarizers and can easily be applied to reflectors with panels of arbitrary size and shape [7J—[12]. However, trihedral corner reflectors whichincorporate reflection polarizers that make substantial use of wire grids and dielectric materialsin their construction are relatively fragile [13]—[14j. Although mechanical and environmentaldamage can be prevented through the use of protective covers and weatherproof seals, suchmeasures substantially increase manufacturing costs and are not always effective.Many of the limitations of previous schemes for altering the polarization response of trihedral corner reflectors can be overcome [15j—[16J by utilizing a reflection polarizer derived fromconducting fins or corrugations [17]—{22]. In section 4.2, the scattering properties of conducting gratings with rectangular grooves are reviewed and design curves for twist polarizing andcircularly polarizing trihedral corner reflectors are presented. Methods for realizing conductinggratings with linear polarization selective responses are proposed. In section 4.3, an algorithmfor predicting the contribution of triple-bounce reflections to the polarization scattering matrixof a trihedral corner reflector which has been modified by the addition of conducting fins orcorrugations to one of its interior surfaces as a function of the direction of incidence and orientation of the reflector is described. In section 4.4, the response patterns of prototype reflectorswith regular, twist-polarizing, and circularly polarizing responses are compared to theoreticalpredictions. The results show that the prototype reflectors respond essentially as expected.Chapter 4. Depolarizing Trihedral Corner Reflectors924.2 Scattering by a Conducting Gratingwith Rectangular GroovesConsider a plane wave incident on a conducting gratingwith rectangular grooves as shown inFigure 4.2. Let the direction of the grooves define the grating axis. Ifthe grating is of infiniteextent, the scattered field will consist of a finite number of propagating planewaves or diffractedorders and an infinite number of nonpropagating or evanescentwaves. If the plane of incidenceis perpendicular to the grating axis, the direction inwhich the mth diffracted order propagatesis given by the grating equation,S1flm =sin+m,m ...,—1,0,1,2,... , (4.1)wherekmis the angle of reflection of the mth diffracted order, 4j is the angle of incidence,A isthe wavelength of the incident wave, and d is the periodof the grating. IfIsinq51, the mthdiffracted order will propagate away from the surface at angle çbwith respect to the x axis.If sinmI> 1, the mthdiffracted order takes the form of an evanescent field which decaysexponentially with increasing distance from the surface of the grating.The number of diffracted orders which are visible is dependent on theperiod of the gratingand the angle of incidence. The chart shown in Figure 4.3 isderived by solving the gratingequation at each angle of incidence for the point at which eachdiffracted order becomes visibleas the period of the grating is gradually increased.Diffracted orders with negative indices willappear atq=90° while those with positive indices will appearat4m= 900.If the gratingperiod and the angle of incidence are chosen such thatsin — < —1, (4.2)then only the specularly reflected order will propagate.However, if the grating period and theangle of incidence are chosen such that the following conditions are satisfied simultaneously,sin — > —1 , (4.3)sin—2 <—1, (4.4)sin+ > +1,(4.5)Chapter 4. Depolarizing Trihedral Corner Reflectors 93V0I-a)00)C0Figure 4.3: Visible diffracted orders as a function of the grating periodincidence j where the plane of incidence is normal to the grating axis.d/\ and the angle ofFigure 4.2: Problemgeometry for scattering by a conducting grating with rectangular grooveswhere d, a, and h are theperiod of the grating and the width and depth of the grooves,respectively, is the angleof incidence, and1mis the angle of reflection of the mth diffractedorder.2.01.51.00.50.00 15 30 45 60Angle of Incidence-q5, (degrees)75 90Chapter 4. Depolarizing Trihedral Corner Reflectors94then both the m = 0 and the m = —1 diffracted orders will propagate.The Bragg anglecondition corresponds to propagation of the m = —1 diffracted order backtowards the sourceand is given byd=2sinq(4.6)It has been found that a grating which satisfies this condition can be made to diffract all thepower in an incident wave into the m = —1 diffracted order by appropriate shaping of its profile.Since the grating is uniform in the z-direction, the spatial derivatives associated with theincident and scattered fields vanish in z and Maxwell’s equations divide into two independentsets,ôE ôE= jwe0E, = —jwç,E9, —= —jwi0H , (4.7)Oy ox Ox 9yOE.OIIOH-h-—= -jwH,-h-- = jwH, -h-— — = jwe0E. (4.8)The first, which consists oflIe,E,and E components, is called a transverse magnetic (TM)or E-polarized field while the second, which consists of E, H, and H components, is called atransverse electric (TE) or H-polarized field. Let [Sm] be the normalized polarization scatteringmatrix which is associated with diffraction into the mth order. Since the TM- and TE-polarizedfields are decoupled, the off diagonal elements of the matrix vanish and[Sm]reduces to0 c nEErc’ 1EE,m EH,m1m[JmJ = = ,C c n rHHHE,m HH,m mwherepErepresents the complex reflection coefficient for diffraction of a TM-polarized incident wave into the TE-polarized component of the mth diffracted order and[Sm]is expressedaccording to the forward scattering alignment (FSA) convention. If the amplitude of the complex reflection coefficients1JE = f0HH1= 1, then the total power in the specular reflectedorder will be independent of the polarization state of the incident wave. Such a response is saidto be polarization-operative. However, if the amplitudeJEE1 JIH1,then the total power inthe specular reflected order will be dependent on the polarization state of the incident wave.Such a response is said to be polarization-selective.Chapter 4. Depolarizing Trihedral Corner Reflectors 95The complex reflection coefficients PEmE and are functions of both the dimensions ofthe grating and the direction and wavelength of the incident wave. Since the free space andgroove regions are defined by separable coordinate systems, a rigorous solution to the problemcan be obtained by representing the fields in each region as the weighted sum of orthonormalbasis functions and determining the relative amplitude and phase of the propagating diffractedorders by mode-matching at the boundary between the regions. A complete derivation ofthe solution and procedures for determining the validity of numerical results obtained by thismethod are presented in Appendix B. Although the solution will invariably converge to anessentially constant result as the number of modes used in the field expansions are systematicallyand gradually increased, it is shown that the manner in which the solution converges and thevalue of the final result will depend on both the number of modes used to represent the fields ineach region and their ratio. This phenomenon is commonly referred to as relative convergence.Numerical results are presented which suggest that the optimum ratio of groove modes to freespace modes is similar in value to the aspect ratio of the grating, aid. An implementation ofthe solution as a pair of subroutines coded in Fortran 77 is also presented.For the purpose of defining the polarization states of the incident and reflected waves, letthe z axis define the local vertical. If period of the grating is less than one-half wavelength,only the specular order will propagate. In such cases, the phase difference 6 between the TEand TM reflection coefficients can be exploited to yield a depolarizing response of the form,[S]=SHH SHy= pEE0 = 1 0(4.10)SVH Svy 0pHH0 e3A phase difference of 0 will yield a regular polarization response while phase differences of 180and 90 degrees will yield twist-polarizing and circularly polarizing responses, respectively. Theco-polar and cross-polar response of regular, twist, and circular polarizers with polarizationresponses defined by (4.10) are plotted as a function of the polarization state of the incidentwave in Figures 4.4, 4.5, and 4.6. If the grating is modified in such a way thatpE0, thelinear polarization selective response of Figure 4.7 will be obtained.Chapter 4. Depolarizing Trihedral Corner Reflectors96(a)b(b)Figure 4.4: Normalized response of a regular reflector as a function ofthe polarization state ofthe incident wave. (a) Co-polar response. (b) Cross-polar response.(a) (b)Figure 4.5: Normalized response of a twist-polarizing reflector as a functionof the polarizationstate of the incident wave. (a) Co-polar response. (b) Cross-polarresponse.Chapter 4. Depolarizing Trihedral Corner Reflectors97(a)Figure 4.6: Normalized response of a circularly polarizing reflector as a functionof the poiarizatiori state of the incident wave. (a) Co-polar response. (b) Cross-polarresponse.(a) (b)Figure 4.7: Normalized response of a vertical polarization selective reflector as a functionof thepolarization state of the incident wave. (a) Co-polar response. (b) Cross-polarresponse.a,(b)Chapter 4. Depolarizing Trihedral Corner Reflectors 98For a grating with a vanishingly small period and arbitrarily thin fins, the groove depthsrequired to realize twist and circularly polarizing responses for normal incidence are simply A/4and .X/8, respectively. In practice, however, account must be taken of both the finite dimensionsof the grating and the angle of incidence. The groove depths required to realize twist andcircularly polarizing responses were determined as a function of the period and aspect ratio ofthe grating and the direction and wavelength of the incident wave by applying a bracketing andbisection algorithm to a subroutine which calculates the phase difference between the TE andTM specular reflection coefficients. The calculations were performed using the mode-matchingformulation presented in Appendix B. Twist polarizer design curves for essentially normalincidence and 45 degree incidence are presented in Figures 4.8 and 4.9, respectively. (Theresponse of the grating could not be calculated for incidence at 4j= 00due to the nature ofthe analytical formulation so the design curves for normal incidence were calculated at qSj= 10instead.) As the grating period approaches zero, the groove depth required to yield a twist-polarizing response converges to A/4 for all aspect ratios and angles of incidence. As the periodof the grating increases, the required groove depth increases for large aspect ratios and decreasesfor small aspect ratios. Circular polarizer design curves for essentially normal incidence and45 degree incidence are presented in Figures 4.10 and 4.11. In this case, the ideal groove depthis satisfactory only for very thin fins and at normal incidence. The difference between thegroove depths required to realize twist polarizers and circular polarizers with identical aspectratios is nearly constant as a function of the grating period. A large difference between therequired groove depths implies that the grating will present the desired polarization responseover a wide operating bandwidth while a small difference implies that the grating will presentthe desired response over a relatively narrow operating bandwidth. For normal incidence, thedifferences between the required groove depths for gratings with aspect ratios a/d of 0.3333,0.7500, and 0.9999, are approximately 0.05 A, 0.10 A, and 0.13 A, respectively. Similar resultsare obtained for 45 degree incidence. This suggests that twist and circular polarizers derivedfrom corrugated surfaces should be designed with the largest possible aspect ratios in order toobtain optimum performance.Chapter 4. Depolarizing Trihedral Corner Reflectors99Three methods for realizing a grating which presents alinear polarization selective responseare depicted in Figure 4.12. A grating with a smallperiod and large aspect ratio which has beenloaded with lossy media, as shown in Figure 4.12(a), will presenta TE-polarization selectiveresponse since the groove region functions as a resonantabsorber for TM-polarized incidentwaves [23], [24]. This approach has several disadvantages including the limitations onthe minimum value of the TM reflection coefficient that can be obtained,the frequency selective natureof the response, and the vulnerability of the lossygroove media to mechanical or environmental damage. Alternatively, aconventional grating which satisfies the Bragg condition can bemade to diffract the TM-polarized component of an incidentwave into the m = —1 diffractedorder while it specularly reflects the TE-polarized componentby appropriate shaping of itsprofile [25]. Although such a grating will act as a retro-reflector forTM-polarized waves andcannot be used to realize a linear polarizationselective trihedral corner reflector, it has beenshown that gratings can be perfectly blazed for other angles of incidence as well, as suggestedby Figure 4.12(b) [26]. The phenomena of perfect blazingfor non-Bragg angle incidence isnot well understood and relatively few examples areknown. Also, the direction of the higherdiffracted order is a function of both the wavelength anddirection of the incident wave andmust be accounted for when designing an apparatus which incorporatessuch a polarizer. Athird approach is shown in Figure 4.12(c). A grating with a smallperiod, large aspect ratio,and a sloped profile which subtends an angle 3, will reflectTM-polarized incident from thebottom of the groove region at an angleqowith respect to the grating normal but will reflectTE-polarized incident waves from the top of the corrugations at anangleo— 23. Since theresponse is independent of the wavelength of the incident wave and thestructure is inherentlyrugged, this approach overcomes many of the limitations of the othertwo schemes. Althoughthe angle 3 should be as large as possible so that the directions of the TMand TE polarizedreflected waves are separated by a wide angle, the maximum height of thefins is given byhm= r sin /3 where r is the maximum dimension of the polarizer in the planeperpendicularto the grating axis. Thus, a practical limit on the slope anglewhich can be accomodated isimposed when a physically large polarizer of this type is required.Chapter 4. Depolarizing Trihedral Corner Reflectors1000.400.35x.0.300.25— y0.200.150.0 0.1 0.2 0.3 0.4Period - d/AFigure 4.8: Twist polarizer design curves for normal incidence.0.400.35X0.300.25—y0.200.150.0 0.1 0.20.3 0.4 0.5Period - d/AFigure 4.9: Twist polarizer design curves for 45 degree incidence.0.5Chapter 4. Depolarizing Trihedral Corner Reflectors1010.250.200.15a,0.050.000.200.150.100.050.000.0 0.2 0.3 0.4Period - d/AFigure 4.11: Circular polarizer design curves for 45 degree incidence.x— ya/d— 0.9999- 0.90000.75000.50000.3333Circular polarizer=1°0.0 0.1 0.2 0.3 0.4 0.5Period - d/AFigure 4.10: Circular polarizer design curves for normal incidence.0.25x— ya/d— 0.9999- 0.90000.75000.50000.3333Circular polarizer0=4500.1 0.5Chapter 4. Depolarizing Trihedral Corner Reflectors 102TE,TM TM TEFigure 4.12: Linear polarization selective reflectors derived from corrugatedsurfaces. (a) Specular reflection of TE-polarized incident waves and dissipation ofTM-polarized incident wavesin lossy groove media. (b) Specular reflection of TE-polarized incident waves and diffractionofTM-polarized incident waves into the m = —1 order.(c) Specular reflection of TM-polarized incident waves with respect to the base of a polarizer with a slopedprofile and specular reflectionof TE-polarized incident waves with respect to the upper surface.TE(a)__J_ _ _(b)TE,TM TE TM(c)Chapter 4. Depolarizing Trihedral Corner Reflectors103Thus far, it has been assumed that the plane ofincidence is perpendicular to the gratingaxis with the phase of the incident wave given byexp(jk(2)= exp[jk(xcosq5 — ysinq1)] .(4.11)Since the grating is uniform along its axis, theresults are easily extended to the case of obliqueincidence [27]. In this case, the phase of the incidentwave is given byexp(jk() = exp[jk(x sin 6 cos—y sin6sin ç+z cos 92)J , (4.12)which is obtained from (4.11) by replacing k by k sin 8: and multiplying by exp(jkz cos 8:).Sincethe scattered field also varies with z as exp(jkzcos82),the spatial derivative ô/ôz inMaxwell’sequations may simply be replaced with jk cos8.As in the case of perpendicular incidence, theequations divide into TM-polarized and TE-polarizedsets. In the TM-polarized case, E 0and the field components are expressed in terms of H whereÔ2H Ô2HZ2 22+2+k sin 8 H = 0 , (4.13)ax ayZ0 (ÔHZ t9HEx+Ey= )ksin28:——z——WY)(4.14)cos8 (aHZ aHHx+Hy= Jksin281—-—X+--—--Y), (4.15)while in the TE-polarized case, li = 0 and the field components areexpressed in terms of EwhereÔ2E Ô2E+(4.16)cos 8 (dEZ OE‘\ErX+Eyy= Jksin28j---—x+--—Y) ,(4.17)1 (0E2 OEZIJx + H‘3kZ0sin2 8X —Y}. (4.18)Since the boundary conditions are identical, the solution to theproblem of scattering by aconducting grating at oblique incidence to the grating axis is identicalto the solution forperpendicular incidence if k is replaced by k sin 8 and all fields are multiplied byexp(jkz cos 8).Chapter 4. Depolarizing Trihedral Corner Reflectors 1044.3 Scattering by a Depolarizing Trihedral Corner ReflectorThe polarization response of a trihedral corner reflector can be modified by replacing one ofits panels by a reflection polarizer derived from a conducting grating with rectangular grooves,as depicted in Figure 4.13. Projections of a typical ray path in the x-y, z-z, and y-z planesare also shown. The coordinate frame is identical to that used in the previous section todescribe scattering by a free-standing grating. In order for this type of depolarizing trihedralcorner reflector to function correctly, the grating must be oriented in such a way that all of therays which are incident from a given direction undergo identical polarization transformationsregardless of the sequence in which they are reflected by each of the interior surfaces. It can beshown that this condition will be satisfied if the grating axis is parallel to one of the principalaxes of the trihedral. If the incident field is resolved into orthogonal components which areTM- and TE-polarized with respect to the grating axis, the components will not be coupleddue to reflection from either the grating or any of the three dihedral corners which comprise thetrihedral corner reflector. As a result, reflection from each interior surface can be represented bya diagonal polarization scattering matrix and the cumulative transformation due to the threereflections necessary to return a incident ray to the source will be independent of the sequencein which the reflections occur. Furthermore, it can be shown that the angles of incidence of therays with respect to both the grating normal and the grating axis will be equivalent regardless ofthe sequence in which the rays are reflected by the interior surfaces of the reflector. Thus, all ofthe rays which are incident on the reflector will undergo identical polarization transformationson reflection from the grating and the required condition will be met.For the purpose of defining the polarization states of incident and reflected waves, hencethe polarization response of the reflector, let the z axis define the local vertical in the reflectorcoordinate frame. For an arbitrary direction of incidence, it is convenient to specify horizontaland vertical in terms of the triad(,Ô,)associated with a conventional spherical coordinatez///(a)(b)Chapter 4. Depolarizing Trihedral Corner Reflectors105z/xyxyx.—x—.y(c)(d)Figure 4.13: Replacement of one panel of a trihedral cornerreflector by a reflection polarizerderived from a corrugated surface. (a) View along the symmetry axis. (b)Projection of atypical ray path in the x-y plane. (c) Projection of a typical raypath in the z-x plane. (d)Projection of a typical ray path in the y-z plane. If thedirection of incidence is (6,4),theangle 0 is given by cos’(cos8cos4) and the angle 6 isgiven by cos’(cos6sin).z zy0.350.30z0.25-caChapter 4. Depolarizing Trihedral Corner Reflectors 1060.50aid0.45—- 0.90000.25x0.20Twist polarizer0 = 54.74°, = 45°0.150.00Figure 4.14: Twist polarizer design curves for incidence along the symmetry axis of a trihedralcorner reflector.0.10 0.20 0.30Period - d/)0.40 0.50:0.10 - 0.90000.7500x . .0.50000.05 Circularpolarizer. 00= 5474b,=4500.00. .1... .1 • ...I. • ..I. •0.00 0.10 0.20 0.30 0.40 0.50Period - d/AFigure 4.15: Circular polarizer design curves for incidence along the symmetry axis of a trihedralcorner reflector.Chapter 4. Depolarizing Trihedral Corner Reflectors 107system where= E= (4.19)= ZXkE= —sin+cos, (4.20)Izxkl= kxh — = —cosqcos8i—sinq5cos9+sin6, (4.21)as described in section 2.2. Here, the component of an incident wave which is TM-polarized withrespect to z corresponds to horizontal polarization while the component which is TE-polarizedwith respect to z corresponds to vertical polarization.The polarization scattering matrix of a depolarizing trihedral corner reflector accordingto the backscatter alignment (BSA) convention is identical to that of a free-standing gratingaccording to the forward scatter alignment (FSA) convention for the same direction of incidence.If only the specular reflected order propagates, the polarization response of the reflector is ofthe form[S]= SHH SHy= 0 = 1 0• (4.22)SVH Svv 0pHH0 e6Since TM- and TE-polarized waves are decoupled upon reflection from the grating even if theplane of incidence makes an oblique angle with the grating axis, the off diagonal elements ofthe polarization scattering matrix given by (4.19) will always vanish in the reflector coordinateframe of Figure 4.13. The polarization scattering matrix of a target relates the scattered fieldESat the receiver to the incident field E at the target according to— —jkrS11 S12 E4 2— rS21 S22 E‘3)where both the fields and the polarization scattering matrix have been expressed with respectto an arbitrarily polarized basis. Often, this expression is normalized by factoring out thescattering cross section of the target and the range dependence of the response. The polarizationmatch factor or polarization efficiency between the scattered field and a receiving antenna ofChapter 4. Depolarizing Trihedral Corner Reflectorsboresightboresightzz,Figure 4.16: (a) Reflector coordinate frame and (b) global coordinate frame.(a)(b)108Figure 4.17: Angle of rotation of the projection of the grating axis onto the view plane forincidence (a) along the symmetry axis and (b) 30 degrees off the symmetry axis.z(a)z,(b)Chapter 4. Depolarizing Trihedral Corner Reflectors109polarization state h is given byIEs.h12P OP1. (4.24)A reflection polarizer designed for incorporation into a trihedral corner reflector shouldpresent the desired polarization response along the same direction of incidence that the reflectorpresents its maximum scattering cross section. For reflectors composed of triangular, circular,or square panels with equal corner lengths, this direction falls along the symmetry axis of thereflector. In the reflector coordinate frame of Figure 4.13, this direction is given in sphericalcoordinates by 8= 54•740,= 45°. Design curves for twist polarizers and circular polarizersfor incidence along this direction are presented in Figures 4.14 and 4.15, respectively. In theprevious section, it was shown that the response of a grating for oblique incidence to the gratingaxis is identical to its response for perpendicular incidence to the grating axis if k is replaced byk sin 6 and all fields are multiplied by exp(jkz cos 8). In this case, sin 8 is given by and cos 8is given by i//. This transformation may also be applied in the design of linear polarizationselective reflectors which are based on the concepts described in the previous section.The response of the target can be altered to suit a particular application by rotating thetarget about its symmetry axis [28). After the target has been rotated by an angle c from thevertical, as shown in Figure 4.19, its polarization scattering matrix is given by[S’j = [R]’ [S] [R] . (4.25)If the response is expressed with respect to a linearly polarized basis, the rotation operator [R]is given bycosu —sina[RL] =, (4.26)sina cosawhile if the response is expressed with respect to a circularly polarized basis, the operator isgiven byeja0[Rc] = . (4.27)0 e3Chapter 4. Depolarizing Trihedral Corner Reflectors 110For example, if a twist-polarizing target is rotated about its symmetry axis to c= 00,±90°,or 180° degrees, it will present a maximum co-polar response for horizontally, vertically,andcircularly polarized incident waves. If the target is rotatedto c = ±45° or ±135°, it will stillpresent a maximum co-polar response to circularly polarized incident waves but willpresent amaximum cross-polar response to horizontally and verticallypolarized incident waves.In practice, it is often necessary to describe the polarizationresponse of a depolarizingtrihedral corner reflector with respect to a global coordinateframe in which the horizontalplane contains the boresight of the reflector and the verticalplane is parallel to the axis of thegrating, as suggested by Figure 4.16. By the methods describedin section 2.4, it can be shownthat the global coordinate framexyzis related to the reflector coordinate frame xyz byI 1 1 1x,/X= * * .—*, (4.28)0 zwhere the x’ axis in the global coordinate frame corresponds tothe symmetry axis of the reflectorand the z’ axis is parallel to vertical. Transformation of the polarization scatteringmatrix ofa depolarizing trihedral corner from the reflector coordinateframe to the global coordinateframe corresponds to rotation of its basis by a prescribed angle which is a functionof both thetransformation matrix which relates the two coordinate frames and the directionof propagation,as described in Chapter 2. Letp define the horizontal plane of the xyz coordinate frame withrespect to the direction of propagation and be given byp=—sin&+cosq, (4.29)and let , similarly define the horizontal plane of the x’y’z’ coordinate frameand be given by(4.30)Let p be the outward normal to the unit sphere at the point P in the zyz coordinateframeand be given by= —sin 6cosq +sinsin+ cos9. (4.31)Chapter 4. Depolarizing Trihedral Corner Reflectors 111The expression for the unit vector , in the x’y’z’ coordinate frame given by (4.30) must betransformed to the xyz coordinate frame. This can be accomplished by determining qY in termsof 6 and according to the relation—3-sin9cos+ -3-sin6sinq— —3-cos0= tan , (4.32)and converting the basis of the unit vector çYp according to the relation given by (4.28). Then,the angle of rotation can be determined from= tan’ . (4.33)\ c5p.14IIf ip is coincident with z axis then the definition ofpgiven by (4.29) is ambiguous. In suchcases, the direction of the horizontal plane with respect to the direction of propagation mustbe defined arbitrarily, as noted in section 2.2. If ip is coincident with the z’ axis, similarconsiderations apply to the definition of given by (4.30). Although several factors causethe polarization response of a depolarizing trihedral corner reflector to degrade asthe directionof incidence shifts away from the boresight, rotation of the projection of the grating axis inthe aperture plane with respect to the local vertical is the most important since it resultsin an effective rotation of the corresponding polarization scattering matrix, as suggested byFigure 4.17.The polarization scattering matrix corresponding to a circular polarization selective response cannot be diagonalized when expressed with respect to a linearly polarized basis. As aresult, it is not possible to modify a conventional trihedral corner reflector to present a circularpolarization selective response simply by adding conducting fins or corrugations of appropriatedimensions and orientation to one of its interior surfaces. Alternative methods for obtainingsuch a response based on the addition of a transmission circular polarizer to a linear polarization selective reflector or a circular polarization selective surface to a twist-polarizing reflectorare proposed in Appendix C.Chapter 4. Depolarizing Trihedral CornerReflectors1124.4 Numerical and ExperimentalResultsIn this section, the results of a testprogram that was conducted in order toevaluate theresponse characteristics of depolarizing trihedralcorner reflectors which incorporate a reflectionpolarizer derived from a corrugated surface are given.Prototype trihedral corner reflectorswhich present regular, twist-polarizing,and circularly polarizing responses were designedforuse at the standard marine radarfrequency of 9.445 GHz and assembled. Their polarizationand azimuthal response patternswere measured and the results were compared to theoreticalpredictions of the contribution of triple-bounce reflectionsto their response based on the analysispresented in the previous section. Details of the designand construction of the prototypereflectors are presented in section 4.4.1. Thepolarization and azimuthal response patterns ofthe prototype reflectors are presented in sections 4.4.2 and 4.4.3,respectively.The response patterns of the prototype reflectorswere measured using the microwave antenna range located on the roof of theElectrical Engineering building at the University ofBritish Columbia. A side view of the outdoorportion of the range is shown in Figure 4.18.The model tower supports the target at a heightof 2.2 metres and travels on a carriage along15 metres of track down the center of the roof of theeast wing of the Electrical Engineeringbuilding. The response pattern of the target will be distortedby near field effects if the rangeto the target is too short or by multipathpropagation effects if the range is too long. Testswee conducted and it was found that the optimum rangeat which to measure the response of aprototype reflector is 11 metres. The model tower canbe configured to either rotate the targetabout a vertical axis for conventional azimuthalpattern measurements or roll the target abouta horizontal axis for polarization response measurementsas shown in Figure 4.19. The CWradar apparatus and the digital pattern recorderwhich were used to measure and record theco-polar and cross-polar response of the target under test weredeveloped specifically for use inthis project. A description of the design andimplementation of the instrumentation, the resultsof tests performed to verify the suitability of the range foruse in the measurement program,and recommendations for future modifications and improvementsare presented in Appendix D.Chapter 4. Depolarizing Trihedral CornerReflectorsOm 5mFigure 4.18: Profile view of theRCS measurement range.0’+(a)(b)— a2mOm113Figure 4.19: Measurement of the response of a prototypetrihedral corner reflector. (a) Polarization response. (b) Azimuthal response.Horn AntennasPenthouse- 3 dBlOm 15m0’Chapter 4. Depolarizing TrihedralCorner Reflectors 1144.4.1 Design and Construction of the Prototype TrihedralCorner ReflectorsA prototype trihedrai corner reflector with a regular polarization responsewas assembled fromtriangular panels with equal corner lengths of 60 cm. At the design frequency of 9.445 GHz,theprototype reflector presents a maximum scattering cross sectionof approximately 540 squaremetres (27 dBsm) with a half-power beamwidth of approximately 38 degrees inboth azimuthand elevation. The reflecting panels were cut from 12 gauge (approximately 2 mm inthickness)utility grade (3003) sheet aluminum and secured to 50 cm lengths of 90°angle aluminum ofapproximate dimensions 1 in x 1 in x in with rivets spaced approximately 5 cm apart. Theright hand panel of the reflector (as viewed from the front) was securedto the angle aluminumwith#6flathead machine screws so that it could be removed easily. A mounting flangewhichmatches the corresponding adapter on the antenna range model tower was secured to therearof the reflector in line with the axis of symmetry. Prototype reflectors with twist-polarizingandcircularly polarizing responses were realized by removing the right hand paneland replacingit with a suitable reflection polarizer. Details of the design andconstruction of the prototypedepolarizing reflectors are shown in Figure 4.20. A photograph of a prototype reflectorwith atwist polarizer installed in place of the right hand panel is shown in Figure 4.21.The extent to which the response of a trihedral corner reflector degrades as the reflectingpanels deviate from perfect flatness and mutual orthogonality has been discussed by Trebits [29]and others. The tolerances on the corner angles of a trihedral corner reflector with triangularpanels are presented as a function of the corner length of the reflector in Table 4.1. Each timethe right hand panel of the prototype reflector was replaced, a set square with a corner lengthof 30 cm was used to verify that all three panels were mutually orthogonaland, if required,suitable adjustments were performed. Although commercially designed and fabricatedtrihedralcorner reflectors intended for use under field conditions are typically manufactured from castalumimum tooling plate between 10 and 15 mm thick in order to realize a structurewhich willmeet the required tolerances [11], [30], [31], the simpler construction technique employed herewas far easier to implement and was deemed adequate for use in the test program.Chapter 4. Depolarizing Trihedral Corner Reflectors115900Angle6JuminumShea ,8juminumMounting FlangeFigure 4.20: Construction details of the prototype depolarizing trihedral cornerreflector.(a) Front view. (b) Rear View.Figure 4.21: Photograph of the prototype twist-polarizingtrihedral corner reflector mountedon the antenna range model tower.Removable Panel(a) (b)Chapter 4. Depolarizing Trihedral Corner Reflectors116Table 4.1:Tolerances on the Corner Angles of a Trihedral Corner Reflector withTriangular Panels.(1)Error in a Single Corner Angle(2) Error in ThreeCorner Angles(2)Corner Length(-1 dB) (-3 dB) (-10 dB) (-1 dB) (-3 dB)(-10 dB)0.15 m (4.7).) ±5.1° ±8.4° ±14.4° ±2.4° ±4.2°±7.5°0.30 m (9.4).) ±2.5° ±4.2° ±7.3° ±1.2° ±2.1°±3.8°0.45 m (14.2).) ±1.7° ±2.8° ±4.9° ±0.8° ±1.4°±2.5°0.60 m (18.9).) ±1.3° ±2.1° ±3.7° ±0.6° ±1.1° ±1.9°1.00 m (31.4).) ±0.8° ±1.3° ±2.2° ±0.4°±0.6° ±1.1°3.00 m (94.3).) ±0.3° ±0.4° ±0.7° ±0.1° ±0.2° ±0.4°Notes: 1. )‘. = 3.18 cm.. 2. which leads to the indicated reduction in the maximum response of thereflector.-d-afl ri[40mmI150mmFigure 4.22: A single segment of a prototype reflection polarizer.Table 4.2: Dimensions of the Prototype Reflection PolarizersType of Polarizer Period d Groove Width a Groove Depth h Aspect Ratio aIdTwist 13.5 mm 12.0 mm 11.5 mm 0.89(0.426).) (0.379).) (0.355).)Circular 13.5 mm 12.0 mm6.0 mm 0.89(0.426).) (0.379).) (0.186).)2mmNote: A = 3.18 cmChapter 4. Depolarizing Trihedral Corner Reflectors 117Reflecting panels which present twist-polarizing and circularly polarizingresponses weredesigned using the methods described in section 4.3 and were assembled from extrusionsthatwere custom manufactured for the purpose by AlcanExtrusions (Richmond, B.C.) from marinegrade (6061) aluminum. A view of a single segment in cross sectionis depicted in Figure 4.22.The closely spaced row of fins along the top of the structure form thereflection polarizer whilethe three longer fins which form the base of the structure merely provide mechanicalsupportduring manufacture and points of attachment during assembly.The extrusion was designed topresent a twist-polarizing response. In order to realize a circularly polarizing response,it wasnecessary to mill the upper row of fins down to the height specified in Table4.2. During themilling procedure, the machinist found it difficult to secure the extrusion to the platformofthe milling machine and problems with vibration were encountered. As a result, theheight ofthe fins varied slightly along the length of the circularly polarizing reflectingpanel. This mayhave contributed to the discrepancies that were observed betweenthe measured response of thecircularly polarizing trihedral corner reflector and theoretical predictions.4.4.2 Polarization ResponseThe polarization response of each of the prototype reflectors was verified by rotating the targetabout its boresight and measuring the co-polar and cross-polar response of the target to avertically polarized incident wave as a function of the angle of rotation c, as suggestedbyFigure 4.19(a). The result corresponds to taking a cross section through the normalized copolar and cross-polar response patterns which are presented in Figures 4.4, 4.5, and 4.6 forellipticity angle e = 0°. as suggested by Figure 4.23. The tilt angle r of the incident wave inthe reflector frame is related to the rotation angle c byr a — 90°,where a = 0° corresponds to vertically polarized incidence. The response patternsof prototypereflectors designed to present regular, twist-polarizing, and circularly polarizing responses arepresented in Figures 4.24, 4.25, and 4.26, respectively. In the case of the regular and twist-Figure 4.23: Evaluation of the polarization response of a radar target. Measuringthe co-polarand cross-polar response of the target to a vertically polarized incidentwave as a function of theangle of rotation r, as suggested by Figure 4.19(a), corresponds to taking a crosssection throughthe normalized co-polar and cross-polar response patternswhich are presented in Figures 4.4,4.5, and 4.6 for ellipticity angle c= 00.Response0 dB-10 dB-20 dB— MeasuredTheoreticalAngle of rotation - a(a) (b)Figure 4.24: Polarization response of a prototype regular trihedral corner reflector as a functionof rotation about the boresight. (a) Co-polar response - vv. (b) Cross-polarresponse- HV.Chapter 4. Depolarizing Trihedral Corner Reflectors118Rigflt circular polanzationUwUHV0.-30dB.08kI’,.08[Chapter 4. Depolarizing Trihedral Corner Reflectors119—MeasuredTheoreticalAngle of rotation - a(a) (b)Figure 4.25: Polarization response of a prototype twist-polarizing trihedral corner reflector as afunction of rotation about the boresight. (a) Co-polar response - aVV. (b) Cross-polar response- HV•—MeasuredTheoreticalAngle of rotation - a(a) (b)Figure 4.26: Polarization response of a prototype circularly-polarizing trihedral corner reflectoras a function of rotation about the boresight. (a) Co-polar response - uvv. (b) Cross-polarresponse -5w0•UHV0•Responseo dB-10dB-20 dB-30dB.08k .08kUw0•ResponseHV0dB-10dB-20 dB-30dB.08kI’).08kChapter 4. Depolarizing Trihedral Corner Reflectors120polarizing reflectors, the measured responses and the theoreticalpredictions agree reasonablywell. In the case of the circularly polarizing reflector, the general form of the measuredresponsesis correct but the discrepancy between measurement andtheory is relatively large. As notedabove, this may have been due to the problems encountered in milling the fins of the circularpolarizer to the correct height.4.4.3 Azimuthal ResponseThe co-polar and cross-polar azimuthal response uvv and HV of each of theprototype tnhedral corner reflectors were measured as a function of the azimuthal angle/‘ for selectedangles of rotation a with respect to the local vertical axis as suggested by Figure4.19. Theresults are compared to theoretical predictions of the contribution oftriple-bounce reflectionsto their response based on the analysis presented in section 4.3. The co-polar andcross-polarazimuthal response patterns of the prototype regular reflector for a rotation angle of 0 degreesare presented in Figures 4.27 and 4.28. The measured responses and thetheoretical predictionsgenerally agree. The co-polar and cross-polar azimuthal response patterns of the prototypetwist-polarizing reflector after it has been rotated to 0 degrees for maximum co-polarresponseare presented in Figures 4.29 and 4.30 while the corresponding response patterns afterthereflector has been rotated to 45 degrees for maximum cross-polar response arepresented inFigures 4.31 and 4.32. Once again, the measured responses and thetheoretical predictionsgenerally agree. The rapid degradation in the polarization response for incidence off the bore-sight is apparent in both cases. The co-polar and cross-polar azimuthalresponse patterns ofa prototype circularly polarizing reflector for rotation angles of 0degrees and 45 degrees arepresented in Figures 4.33 through 4.36. Although the general form of the measured response iscorrect, the measured co-polar azimuthal beamwidth of the prototype reflector isslightly widerthan predicted in both cases. Also, the measured cross-polar azimuthalresponse pattern fora rotation angle of 0 degrees is substantially higher than predicted.However, the measuredcross-polar azimuthal response pattern for a rotation angle of 45 degrees agrees verywell withthe predicted values.121Chapter 4. Depolarizing Trihedral Corner Reflectors0.0-10.0-20.0-30.0-45.0 -30.0 15.030.0 45.0— MeasuredTheoretical** contribution oftriple-bouncereflecbonsAzimuth Angle - ‘(deg)Figure 4.27: Co-polar azimuthal response pattern of a prototype regular trihedral corner reflector for rotation angle a= 00and vertically polarized transmission.0.0-10.0a,>-20.0-15.0 0.0-15.0 0.0— MeasuredTheoretical*-30.0*contribution of-45.0 -30.0 15.0 30.0 45.0triple-bouncereflectionsAzimuth Angle - ‘(deg)Figure 4.28: Cross-polar azimuthal response patterns of a prototype regular trihedral cornerreflector for rotation angle a =00and vertically polarized transmission.Chapter 4. Depolarizing Trihedral Corner Reflectors1220.0-10.0Cl)C..).-20.0-30.0-45.0 -30.0 -15.0— MeasuredTheoretical** contribution of15.0 30.0 450triple-bouncereflectionsAzimuth Angle - (deg)Figure 4.29: Co-polar azimuthal response pattern of a prototype twist-polarizing trihedralcorner reflector for rotation angle c— 00and vertically polarized transmission.0.0-10.0Cl)C)4--20.0-30.0-45.0 -30.0 -15.0 0.0 15.0 30.0 45.0— MeasuredTheoretical**contribution oftriple-bouncereflectionsAzimuth Angle - ‘(deg)Figure 4.30: Cross-polar azimuthal response patterns of a prototype twist-polarizing trihedralcorner reflector for rotation angle c = 00and vertically polarized transmission.0.0Chapter 4. Depolarizing Trihedral Corner Reflectors1230.0-10.0Cl)C)-I-20.00.0? -10.0Cl)C)a,>-20.0-30.0-45.0 -30.0 -15.0 0.0 15.0 30.0 45.0— MeasuredTheoretical** contribution oftriple-bouncereflectionsAzimuth Angle - ‘(deg)Figure 4.32: Cross-polar azimuthal response patterns of a prototype twist-polarizingtrihedralcorner reflector for rotation angle a= 450and vertically polarized transmission.— MeasuredTheoretical*-30.0*contribution or-45.0 -30.0 -15.0 0.0 15.0 30.0 45.0triple-bouncereflectionsAzimuth Angle - ‘(deg)Figure 4.31: Co-polar azimuthal response pattern of a prototype twist-polarizingtrihedralcorner reflector for rotation angle a = 45° and vertically polarized transmission.Chapter 4. Depolarizing Trihedral Corner Reflectors124Cl)C)a,a,*contribution oftriple-bouncereflectionsFigure 4.34: Cross-polar azimuthal response patterns of a prototype circularly-polarizing trihedral corner reflector for rotation angle a= 00and vertically polarized transmission.0.0-10.0-20.0-30.0-45.0 -30.0 -15.0— MeasuredTheoretical**contribution of0.0 150 300 450triple-bouncereflectionsAzimuth Angle -V(deg)Figure 4.33: Co-polar azimuthal response pattern of a prototype circularly-polarizing trihedralcorner reflector for rotation angle a= 00and vertically polarized transmission.0.0-10.020.0—MeasuredTheoretical*-30.0-45.0 -30.0 -15.0 0.0 15.0 30.0 45.0Azimuth Angle - ‘(deg)Chapter 4. Depolarizing Trihedral Corner Reflectors 125Azimuth Angle - ‘(deg)Figure 4.35: Co-polar azimuthal response pattern of a prototype circularly-polarizing trihedralcorner reflector for rotation angle c =450and vertically polarized transmission.-30.0*contribution of-45.0 -30.0 -15.0 0.0 15.0 30.0 45.0triple-bouncereflectionsAzimuth Angle - ‘(deg)Figure 4.36: Cross-polar azimuthal response patterns of a prototype circularly-polarizing trihedral corner reflector for rotation angle o = 45° and vertically polarized transmission.:9-Cl)00.0-10.0-20.0-30.0-45.0 -30.0 -15.0 0.0 15.0— MeasuredTheoretical**contribution oftriple-bouncereflections30.0 45.00.0:9-i-10.0-20.0— MeasuredTheoretical*Chapter 4. Depolarizing Trihedral Corner Reflectors 1264.5 ConclusionsA method for altering the polarization response of a conventional trihedral corner reflector byreplacing one of its reflecting panels by a reflection polarizer derived from a conducting gratingwith rectangular grooves of appropriate dimensions and orientation has been proposed. Sincea depolarizing reflector of this type does not make use of wire grids or dielectric materials in itsconstruction, it is less vulnerable to mechanical or environmental damage than other schemeswhich have been proposed in recent years. The scattering properties of conducting gratings withrectangular grooves have been reviewed and analytical solutions to the problem of scatteringby a conducting grating with rectangular grooves have been derived for the cases of TM- andTE-polarized incident waves by mode-matching between the free space and groove regions. Theresults, including an implementation of the solution as a pair of subroutines coded in Fortran77, are presented in Appendix B.If the axis of the fins or corrugations are aligned with one of the axes of the trihedral thenall rays incident from a given direction will experience the same polarization transformationupon reflection from the grating regardless of the sequence in which each ray is reflected fromeach of the three interior surfaces of the reflector. If the period of the grating is less than one-half wavelength and the reflector is oriented so that the projection of the grating axis into theaperture plane is oriented vertically, the phase difference 5 between the TM and TE reflectioncoefficients of the specular reflected order can be exploited to yield a depolarizing responsealong the reflector boresight of the form,[S]=SHH SHy = 1 0SVH Svv 0 ewhere S is the ratio of the electric field components of the specular reflected order and the incident field. A phase difference of 0 yields a regular polarization response while phase differencesof 180 and 90 degrees yield twist polarizing and circularly polarizing responses, respectively.Since the grating is uniform in one dimension, the TM- and TE-polarized fields are decoupledChapter 4. Depolarizing Trihedral Corner Reflectors127and SHy =SVH = 0. Design curves for trihedral corner reflectorswhich present circularly polarizing and twist polarizing responses along their boresighthave been given. The polarizationresponse of the reflector can be modified by rotating the reflector about its boresightby anangle c, yielding a response of the form,cos sinc 1 0 coscr —sino—sina cosc 0 e6 sinc cosüMethods for realizing trihedral corner reflectors which present a linearpolarization selectiveresponse of the form,[Sj=SHH SHy= 1 0SVHSyy 0 0have been proposed. The problem of realizing trihedral corner reflectors which present acircularpolarization selective response is considered in Appendix C.An algorithm for predicting the contribution of triple-bounce reflectionsto the polarizationscattering matrix of a modified trihedral corner reflector as a function of thedirection of incidence and the orientation of the reflector has been described. Although severalfactors causethe polarization response to degrade as the direction of incidenceshifts away from the reflectorboresight, rotation of the projection of the grating axis in the aperture planewith respect to thelocal vertical is the most important since it results in an effective rotation of thecorrespondingpolarization scattering matrix. The polarization and azimuthal response patternsof prototypereflectors with regular, twist polarizing, and circularly polarizing responses were measuredusingthe experimental facility which is described in Appendix D.The results show that the reflectorsrespond essentially as predicted.References[1] J. Croney, “Civil Marine Radar,” in Radar Handbook, M.I. Skolnik, Ed., New York:McGraw-Hill, 1970, chap. 31.[2] S. H. Yueh, J. A. Kong, R. M. Barnes, and R. T. Shin, “Calibration of polarimetric radarsusing in-scene reflectors,” J. Electromagn. Waves Appl., vol. 4,pp.27—48, Jan. 1990.[3] E.M. Kennaugh, Dielectric Reflector. U.S. Patent No. 2,872,675, dated Feb. 3, 1959.[4] F.M. Weil, M.L. Ingaisbe, R.E. Stein, and J.G. McCann, Radar Reflector for CircularlyPolarized Radiation. U. S. Patent No. 2,786,198, dated Mar. 19, 1957.[5] S.J. Blank and L.H. Sacks, Radar Target for Circularly Polarized Radiation. U.S. PatentNo. 3,309,705, dated Mar. 14, 1967.[6] J.A. Scheer, “Radar Reflectivity Calibration Procedures,” in Radar Reflectivity Measurement: Techniques and Applications, N.C. Currie, Ed. Norwood, MA: Artech House, 1989,p.109.[7] E.M. Kennaugh, A Corner Reflector for Use with Circularly Polarized Radars. Columbus,OH: Ohio State Univ., Rep. 601-8, 30 Jan. 1956. (cited in C.G. Bachman, Radar Targets.Lexington, MA: Heath, 1982,pp.71—82.)[8] P.Münzer, Radar-Reflektor für zirkular, elliptisch oder in beliebiger Ebene linear polarisierte elektrornagnetische Wellen. (Radar reflector for circularly, elliptically, and mostlinearly polarized electromagtnetic waves.) Federal German Patent No. 1,196,255, datedJul. 8, 1965.[9] A. Macikunas, S. Haykin, and T. Greenlay, Trihedral Radar Reflector. U.S. Patent No.4,843,396, dated Jun. 27, 1989.128References129[10) A. Macikunas and S. Haykin, “Trihedral twist-grid polarimetricreflector,” lEE Proc. F,vol. 140, pp. 216—222, Aug. 1993.[11] E. Johansen and A. Fromm. Depolarizing Radar CornerReflector. U.S. Patent No.4,724,436, dated Feb. 9, 1988.[12] D.R. Sheen, E.L. Johansen, L.P. Elenbogen, and E.S. Kasischke, “The griddedtrihedral: a new polarimetric SAR calibration reflector,” IEEE Trans.Geosci. Remote Sensing,vol. GE-30,pp.1149—1153, Mar. 1990.[13] P.W. Hannan, Twistreflector. U.S. Patent No. 3,161,879, datedDec. 15, 1964.[14] L.G. Josefsson, “A broad-band twist reflector,” IEEE Trans.Antennas Propagat., vol. AP19,pp.552—554, July 1971.[15] D.G. Michelson and E.V. Jull, “A Depolarizing Calibration Targetfor Radar Polarimetry,”Proc. IGARSS’90 (College Park, MD), May 20—24, 1990,p.799.[16] D.G. Michelson and E.V. Jull, “Depolarizing Trihedral Corner Reflectors,”IEEE/AP-SSymp. Dig. (Ann Arbor, MI), June 28—July 2, 1993,pp.238—241.[17] J.D. Hanfling, G. Jerenic, and L.R.Lewis, “Twist reflector design using E-type and H-typemodes,” IEEE Trans. Antennas Propagat., vol. AP-29,pp.622—629, July 1981.[18] W.B. Offut and L.K. DeSize, “Methods of Polarization Synthesis,” in AntennaEngineeringHandbook, 2nd ed., R. C. Johnson and H. Jasik, Eds. New York: McGraw-Hill, 1984,chap. 23.[19] E.V. Jull, “Reflection circular polarisers,” Electron. Lett.,vol. 15,pp.423—424, July 5, 1979.[20] R. Kastner and R. Mittra, “A spectral-iteration technique for analyzinga corrugatedsurface twist polarizer for scanning reflector antennas,” IEEE Trans. AntennasPropagat.,vol. AP-30,pp.673—676, July 1982.References130[21] D.G. Michelson, P. Phu, and E.V. Jull, “Millimetre-wavegrating polarizers,” 22nd URSIGeneral Assembly, Tel Aviv, Israel, Aug. 1987.[22] A. Hessel, J. Shmoys, and D.Y. Tseng, “Bragg-angle blazing ofdiffraction gratings,”J. Opt. Soc. Am., vol. 65, pp. 380—384, Apr. 1975.[23] E.F. Knott, J.F. Shaeffer, M.T. Tuley, Radar Cross Section: Its Prediction, Measurement,and Reduction. Norwood, MA: Artech House, 1985,pp.244—247.[24] A.K. Bhattacharyya and D.L. Sengupta, Radar Cross Section Analysis and Control.Nor-wood, MA: Artech House, 1991,pp.208—212.[25] E.V. Jull and J.W. Heath, “Reflection grating polarizers,” IEEE Trans.Antennas Propagat., vol. AP-28,pp.586—588, July 1980.[26] E.V. Jull and J.W. Heath, “Radio applications of rectangulargroove corrugations,” inIEEE/AP-S Symp. Dig. (Seattle, WA), June 1979,pp.515-518.[27] E.V. Jull, Aperture Antennas and Diffraction Theory. Stevenage, UK:Peregrinus, 1981,pp.80—81.[28] S.H. Bickel, “Some invariant properties of the polarization scattering matrix,” Proc. IEEE,vol. 53,pp.1070—1072, 1965.[29] R.N. Trebits, “Radar Cross Section,” in Radar Reflectivity Measurement: TechniquesandApplications, N.C. Currie, Ed. Norwood, MA: Artech House, 1989,pp.46—49.[30] A. Freeman, Y. Shen, and C. L. Werner, “Polarimetric SAR calibrationexperiment usingactive radar reflectors,” IEEE Trans. Geosci. Remote Sensing, vol. GE-28,pp.224—240,Mar. 1990.[31] Precision Radar Corner Reflectors, Product Brochure No. 9187F. Norcross, GA: SpectrumTechnologies International, 1991.Chapter 5SUMMARY, CONCLUSIONS, AND RECOMMENDATIONS5.1 Summary and ConclusionsIt is often necessary to enhance the radar cross section of a cooperative target either to increasethe maximum range at which the target can be reliably detected or to provide a target witha known response which may be used to assist in radar calibration and performance verification. In recent years, a requirement has arisen for rugged yet inexpensive radar reflectorswhich present both a very large scattering cross section and a specified polarizationresponseover a wide angular range for use as calibration targets in airborne and spaceborne imagingradar systems for geophysical remote sensing and location markers in radar-assisted positioningsystems for marine navigation. Passive targets, such as trihedral corner reflectors,are oftenbetter suited for use in such applications than active targets since they do not require anexternal source of power and are inherently more reliable. This study considers several problemsrelated to the analysis, design, and implementation of passive targets including transformation of the polarization response of a target or the polarization state of an antenna betweencoordinate frames, modification of the angular coverage of a conventional trihedral corner reflector by appropriate shaping of its reflecting panels, modification of the polarization responseof a conventional trihedral corner reflector by the addition of conducting fins or corrugationsof appropriate dimensions and orientation to one of its interior surfaces, design of top hat reflectors with specified response characteristics, and design of reflection polarizers derivedfromconducting gratings with rectangular grooves.131Chapter 5. Summary, Conclusions, and Recommendations 132In Chapter 2, it is shown that transformation of a polarization descriptor between coordinate frames corresponds to rotation of its basis by a prescribed angle which is a functionof both the transformation matrix that relates the two coordinate frames and the direction ofpropagation. Two methods for determining the angle of rotation for the case in which the localvertical is defined by the direction in each frame are derived using spherical trigonometry andvector algebra, respectively. Both methods are robust and will yield the correct result but themethod based on vector algebra is more compact and would be easier to implement in software.Although the elements of the coordinate transformation matrix can be determined from eitherthe relative directions of the three principal axes in each coordinate frame or the Euler angleswhich define a series of rotations which will transform one coordinate frame into the other, itmaybe difficult to obtain these parameters in practice. A third method is derived which overcomes this limitation by allowing the elements of the coordinate transformation matrix to bedetermined from any pair of dirctions which have been expressed in terms of both coordinateframes. Algorithms for rotating the basis of several polarization coordinates in common use,including the complex polarization ratio, the complex polarization vector, the Stokes vectorand several of its variants, the coherency matrix, the polarization scattering matrix, and theMueller matrix and several of its variants, are presented.In Chapter 3, the problem of predicting the scattering cross section and angular coverage ofa conventional trihedral corner reflector with panels of completely arbitrary shape is considered.Although general purpose numerical techniques such as the finite-difference time-domain (FDTD) and the shooting and bouncing ray (SBR) methods can predict both the contributionsof single, double, and triple-bounce reflections from the interior and the effect of deviationsof the panels from perfect flatness and mutual orthogonality, they are extremely demandingcomputationally. For most purposes, it is sufficient to account for the contribution of triplebounce reflections from the interior since they completely dominate the response for mostdirections of incidence. The empirical model originally proposed by Spencer for predicting theequivalent flat plate area of ideal trihedral corner reflectors with either triangular or squareChapter 5. Summary, Conclusions, and Recommendations133panels (and the numerical implementation of the modelthat was devised by Keen) may failwithout indication when applied to reflectors with complex panel shapes. Asimple yet robustnumerical prediction algorithm which overcomes this limitation is formulatedby applying theapproach used by Keen and a polygon-clipping algorithm formulated by Weiler andAtherton toan alternative model originally proposed by Robertson. Theresponse patterns of ideal trihedralcorner reflectors which present three-fold symmetry about the boresight areplotted on equalarea projection grids and compared. If the corner length of the reflector is fixed, it isgenerallyfound that attempts to increase the beamwidth of the response by modifying theshape ofthe panels are accompanied by a reduction in the maximum response andvice versa. A setof design curves for bilaterally symmetric reflectors which are composedsolely of triangular,elliptical, or rectangular reflecting panels are derived and theresponse patterns of selectedbilaterally symmetric reflectors which are composed of combinations ofpanels with variousshapes including triangular, circular, and square are presented andcompared. It is found thatif three-fold symmetry is broken so that the reflector simply presentsbilateral symmetry abouta mirror plane containing the boresight then the beamwidth of theresponse in one principalplane can be increased relative to the beamwidth in the orthogonal plane bymodifying theshape of the reflector panels. The additional degree of freedom is shown to be usefulwhen areflector must be designed subject to a constraint such as a requirement that it present aplanaraperture in order to facilitate the attachment of either a transmission polarizer or aprotectivecover.In Chapter 4, a method for altering the polarization response of aconventional trihedralcorner reflector by adding conducting fins or corrugations ofappropriate dimensions and orientation to one of its interior surfaces is proposed. Since depolarizing reflectors of this typedo notmake use of wire grids and dielectric materials in their construction, they are lessvulnerable tomechanical or environmental damage than other schemes which have been proposed in recentyears. They also avoid many of the mechanical problems associated with the attachment ofatransmission polarizer across the reflector aperture. If the axis of the fins or corrugations areChapter 5. Summary, Conclusions, and Recommendations134aligned with one of the axes of the trihedral then all raysincident from a given direction willexperience the same polarization transformation upon reflectionfrom the polarizer regardless ofthe sequence in which each ray is reflected from each of the three interiorsurfaces. If period ofthe grating is less than one-half wavelength and the reflector isoriented so that the projectionof the grating axis into the aperture planeis oriented vertically, the phase difference t betweenthe TM and TE reflection coefficients of the specular reflected order can beexploited to yielda depolarizing response along the boresight of the reflector of theform,[S]=tSHH SHy= 1, (5.1)SVH Svv 0e’6where S is the ratio of the electric field components of the specularreflected order and theincident field. A phase difference of 180 degrees yields a twist-polarizing responsewhile a phasedifference of 90 degrees yields a circularly polarizing response. Sincethe grating is uniform inone dimension, the TM and TE polarizations are decoupled and SHy = SVH = 0.Design curvesfor trihedral corner reflectors which present twist-polarizingand circularly polarizing responsesalong their boresight are given. The polarizationresponse of the reflector can be modified byrotating the target about its boresight by an angle a, yielding aresponse of the form,cosa sina 1 0 cosa —sina[S9=. (5.2)— sin a cos a 0 e6 sin acos aMethods for realizing linear polarization selective trihedral cornerreflectors using similar techniques are proposed. An algorithm for predicting the contribution oftriple-bounce reflectionsto the polarization scattering matrix of a modified trihedralcorner reflector as a function ofthe direction of incidence and the orientation of thereflector is derived. Although severalfactors cause the polarization response to degrade as thedirection of incidence shifts awayfrom the boresight, rotation of the projection of the grating axis in the apertureplane withrespect to the local vertical is the most important since it results in an effectiverotation ofthe corresponding polarization scattering matrix.Prototype trihedral corner reflectors withChapter 5. Summary, Conclusions, and Recommendations 135regular, twist-polarizing, and circularly polarizing responses were designed and assembled. Experimental results are presented which show that the prototype reflectors respond essentiallyas predicted.In Appendix A, the problem of designing a top hat reflector with specified response characteristics is considered. Since the top hat reflector is a body of revolution, its response pattern isuniform in azimuth. Expressions for the elevation response pattern, maximum scattering crosssection, angle of maximum response, and 1 and 3 dB elevation beamwidths of a top hat reflectorare derived and design curves are presented. As the angle of maximum response becomes appreciably greater or less than 45 degrees, the elevation response pattern will become increasinglyasymmetrical. In such cases, it may be preferable to consider the angle of median response fora given elevation beamwidth rather than the angle of maximum response. Expressions for theangle of median response are derived and design curves are presented. The results are used tosolve a sample design problem.In Appendix B, the problem of scattering by a conducting grating with rectangular groovesis considered. Analytical solutions are derived for the cases of TM- and TE-polarized incidentwaves by mode-matching between the free space and groove regions. Procedures for determining the validity of numerical results which are obtained using this formulation are discussed.Although the solution will invariably converge to an essentially constant result as the numberof modes used in the field expansion are systematically and gradually increased, it is shownthat the manner in which the solution converges and the value of the final result will dependon both the number of modes used to represent the fields in each region and their ratio. Thephenomenon is commonly referred to as relative convergence. Numerical results are presentedwhich suggest that the optimum ratio of groove modes to free space modes is similar in valueto the ratio of the groove width to the grating period. An implementation of the analyticalsolutions as a pair of subroutines coded in Fortran 77 is presented.In Appendix C, the problem of modifying a conventional trihedral corner reflector to presenta circular polarization selective response is considered. It is shown that such a response cannotChapter 5. Summary, Conclusions, and Recommendations136be realized using the techniquesdescribed in Chapter 4 because the correspondingpolarizationscattering matrix cannot be diagonalized whenexpressed with respect to a linearly polarizedbasis. Alternative methods for obtainingsuch a response based on the addition of atransmissioncircular polarizer to a linear polarizationselective reflector or a circular polarizationselectivesurface to a twist-polarizing reflectorare proposed. The performance of suchreflectors andtheir usefulness in practice will dependon a number of factors that are not consideredhereincluding (1) degradation of the polarizationresponse of either of the polarizers forincidenceoff the reflector boresight,(2) possible degradation of the polarizationresponse due to multiplereflections between the transmission and reflection polarizers,and (3) the mechanical ruggednessof the modified reflector and the ease withwhich it can be fabricated.In Appendix D, the experimentalfacility which was set up to measure theresponse ofprototype trihedral corner reflectors is brieflydescribed. Details of the physical layoutand thedesign and implementation of the CW radarapparatus and digital pattern recorder aregiven.Test results show that the facilityis suitable for use in the prototype reflectormeasurementprogram. Before the facility is usedin future measurement programs, however,considerationshould be given to replacing the existing microwavereceiver with a newer model which is lessprone to drift, correcting the backlash inthe model tower positioner head, and evaluatingtheuse of either a berm or a clutterfence to reduce multipath effects.5.2 Recommendations for Further WorkAlthough depolarizing trihedral corner reflectorswhich incorporate conducting fins or corrugations along one of their interior surfaceswill find immediate use as location markersandcalibration targets in polarimetricradar systems used for radar-assisted positioningand geophysical remote sensing, consideration shouldalso be given to using reflectors of thistype asradar cross section enhancementdevices in conventional marine radar navigation. Anincreasing number of civil marine radars areequipped with provision for transmitting andreceivingcircularly polarized waves in order to suppressrain clutter. However, mariners must exerciseChapter 5. Summary, Conclusions, and Recommendations 137caution when using them in this mode because the conventional trihedral corner reflectors whichare used to augment the response of most buoys, channel markers, and navigations hazards areessentially invisible to a circularly polarized radar. Replacement of conventional reflectors bydepolarizing reflectors which respond equally well to horizontally, vertically, and circularly polarized incident waves would significantly enhance the safety of marine navigation. Furtherengineering studies and field trials should be conducted in order to verify the suitability of depolarizing reflectors of the type proposed in Chapter 4 for use in such applications. Since suchreflectors would be required in large numbers, it would also be desirable to study alternativemethods of manufacture with the aim of minimizing their unit cost.Detailed analysis of modified trihedral corner reflectors for engineering purposes would besimplified if the radar cross section prediction algorithms and software tools described in Chapters 3 and 4 were combined into a single integrated software package designed to run on a highperformance UNIX workstation or other platform suitable for numerically intensive computing. A package designed to assist in the evaluation of alternative schemes for truncating andcompensating trihedral corner reflectors should include implementations of the prediction algorithm for ideal trihedral corner reflectors based on Robertson’s model (as described in Chapter3), either the finite-difference time-domain (FD-TD) or the shooting and bouncing ray (SBR)methods for predicting the response of non-ideal trihedral corner reflectors, and suitable toolsfor entering the reflector geometry and displaying the results of the calculations. The FD-TDand SBR methods are sufficiently demanding of computing resources that it will probably benecessary to include provision for running them either as background tasks for extended periodsor remotely on a sufficiently powerful host. A package designed to assist in the design of depolarizing trihedral corner reflectors which incorporate conducting fins or corrugations would alsorequire code for predicting the response of a conducting grating with rectangular grooves. Thepackages should also include a facility for comparing numerical predictions with experimentalresults obtained by the user.Chapter 5. Summary, Conclusions, and Recommendations138Although recent trials conducted by Transport Canada (TransportationDevelopment Centre, Montreal, P.Q.) have shown that radar-assisted positioning systemsbased on measurementof the range and bearing to shore-based cooperative targets can providevessels navigating ininland waterways, harbours, and harbour approaches with veryaccurate position data underoptimum conditions, the lack of data concerning the scattering statisticsof terrain at grazingincidence makes it difficult to give reliable estimates of the size of targets requiredto achievespecified probabilities of detection and false alarm. Future developmentof radar-assisted positioning systems for marine navigation should assign a high priority tothe collection of suchdata in a variety of representative clutter environments in order toassist in system planning andevaluation. Future trials should include provision for determining theconfidence with whichthe cooperative targets are detected against the clutter background inorder to provide somemeasure of the reliability of such systems.Three methods for realizing linear polarization selective reflectivesurfaces derived fromconducting gratings with rectangular grooves are suggestedin Chapter 4 including (1) gratingsdesigned to reflect TE-polarized incident waves while diffracting TM-polarized incidentwavesinto the m = —1 diffracted order, (2) gratings designed to reflect TE-polarizedwaves whiledissipating TM-polarized incident waves in lossy groove media, and (3) gratings designedtoreflect TE- and TM-polarized waves in different directions by utilizing asloped grating profile.Each type of grating suffers from various limitations. The scattering characteristicsof eachscheme should be examined in further detail so that their relative performance andmerits canbe assessed in the context of both realizing linear and circular polarizationselective trihedralcorner reflectors and their possible use in other quasi-optical systems.The first scheme canbe analyzed simply using the mode-matching formulation presented in Appendix Band someresults have been presented in the literature. Analysis of the secondscheme would requirethat the representation for the fields in the groove region be modified in order to accountfor the effects of the lossy groove media. In both cases, the impact of relativeconvergenceon calculation of the amplitude of the diffracted orders should be assessed. Sincethe problemChapter 5. Summary, Conclusions, and Recommendations139geometry of the third scheme is not periodic and the free spaceand groove regions do not belongto separable coordinate systems, it will necessary to employ a different approach.Although apurely numerical technique could be used in such an investigation, it may be simpler to beginsimply by setting up a design parameter matrix, constructing a set of prototypereflectionpolarizers, then measuring and comparing the response characteristics of the prototypes.The task of coding solutions to antenna or scattering problems which involve theprediction of polarization dependent effects would be simplified if astandard library of functionsand subroutines was available for performing polarization calculations. Sucha library wouldinclude routines for transforming polarization descriptors betweencoordinate frames based onthe results presented in this study and routines for converting polarization descriptorsfrom onetype to another (e.g., polarization scattering matrix to Mueller matrix),converting polarization descriptors between the backscatter and forward scatter alignmentconventions, verifyingthe validity of polarization descriptor data, calculatingpolarization efficiency, and extractingco-polar and cross-polar nulls and related parameters frompolarization response descriptorsbased on results generally available in the literature. Since these routines will makeextensiveuse of matrices and complex variables, the library should initially be coded inFortran 77. Consideration should also be given to producing libraries for use with more specialized numericalanalysis and symbolic algebra packages such as Matlab (The Mathworks, Inc.,Natick, MA),Maple (Waterloo Maple Software, Inc., Waterloo, Ont.), Macsyma (Macsyma, Inc., Arlington,MA), and Mathematica (Wolfram Research, Inc., Champaign, IL).Appendix ADESIGN CURVES FOR TOP HAT REFLECTORSA.1 IntroductionThe top hat reflector is a variant of the conventional dihedral corner reflector which has recentlyfound use as a calibration target for airborne imaging radars. It may be regarded as a body ofrevolution formed by rotating a dihedral corner about an axis of symmetry which is perpendicular to both the seam of the dihedral and one of its two reflector panels. As such, its response isindependent of azimuth angle. In practice, a top hat reflector is usually realized by attaching aright circular cylinder to a flat circular plate such that the cylinder axis is perpendicular to andpasses through the center of the plate although other ground plane configurations have beenconsidered. The requirements for the cylinder and base to be flat and mutually orthogonal aresimilar to those for the panels of a conventional dihedral corner reflector.Although the top hat reflector was apparently introduced over twenty-five years ago, it wasnot discussed in the literature until the early 1980’s [1]. It was not treated in standard referencesuntil much later still and, even so, the discussion is brief and contains obvious errors [2], [3].Blejer [4] recently extended the simplified model introduced by Johansen [1] to account for thecontribution of single bounce reflections from the cylinder, the cylinder cap, and the annularbase. The design of top hat reflectors is complicated by the dependence of the scattering crosssection and the angle of maximum response on three parameters: the radius and height ofthe cylinder and the width of the annular base. In addition, the main response lobe becomesincreasingly asymmetrical as the angle of maximum response becomes either much greater ormuch less than 45 degrees to the vertical. Design curves and related material are not generallyavailable.140Appendix A. Design Curves for Top Hat Reflectors 141In this appendix, the problem of designing top hat reflectors with specified response characteristics is considered. In section A.2, an approximate expression for the scattering cross sectionof a top hat reflector is derived as a function of the elevation angleof the incident ray usinga simplified physical optics model which accounts only for the contributionof double-bouncereflections from the cylinder and the ground plane. In section A.3, design curveswhich relatethe physical dimensions of a top hat reflector to the direction of its maximumresponse, itsmaximum scattering cross section, and its elevation beamwidthare presented. In section A.4,the results are used to solve a sample design problem.A.2 AnalysisConsider a top hat reflector which consists of a right circular cylinder of radiusa and heightc attached to a circular ground plane of radius b as shownin Figure A.1. The cylinder axiscoincides with the z-axis and the ground plane lies in the x-y plane. The angleof elevation 8is measured from the positive z-axis. If the plane of incidence contains thecylinder axis andthe opposite panel is of sufficient extent,a top hat reflector will return an incident ray to thesource in the same manner as would a dihedral corner. Since a top hat reflectoris a body ofrevolution, its response pattern is uniform in azimuth.The locus of maximum response fordouble-bounce reflections is an inverted circular cone which intersects the cylinderin the x-yplane.Figure A.1: Problem geometry for scattering by a top hat reflector.Appendix A. Design Curves for Top Hat Reflectors 142The contribution of double-bounce reflections to the scattering cross section of a top hatreflector can be determined using a ray-optical model in which the ground plane is treated as animage plane of finite extent as shown in Figure A.2. If an incident ray intercepts the cylinderat height h, radius a, and angle 9, the reflected ray will intercept the x-y plane at a radiush tan 9+a. If h> (b — a) cot 9, the point of interception will occur beyond the ground planeradius b and the reflected ray will not be returned to the source. The effective height c’ of theportion of the cylinder which contributes to the response is therefore a function of the angle ofincidence and is given byO<07,(A.1)(c,y 9<ir/2,where(b—a’7=tan). (A.2)TOP-HATREFLECTORIMAGE2c 2cNNNFigure A.2: A simplified ray-optical model for scattering by a top hat reflector.Appendix A. Design Curves for Top Hat Reflectors 143If a top hat reflector is sufficiently large in terms of wavelength, the contribution of singlebounce reflections from the cylinder, the cylinder cap, and the ground plane of the reflector tothe backscatter response will be extremely small for other than vertical or horizontal incidenceand can be neglected for the purposes of design. This reduces the problem to one of simplypredicting the forward scattering cross section of a cylinder of radius a and height 2c’ where theplane of incidence contains the cylinder axis. The corresponding problem geometry is shownin Figure A.3. The intersection of the plane of incidence with the surface of the cylinder isreferred to as the specular line. The unit vectors i and . give the directions of the incident andscattered rays. The unit vectors ñ and ñ0 denote the outward normal to the cylinder at anypoint on its surface and along the specular line, respectively, and are given byn = cos+sinq, (A.3)andñ0=cos?SI+sinq. (A.4)specularline\EZJFigure A.3: Problem geometry for forward scattering by a cylinder.Following Knott [5], an approximate expression for the bistatic scattering cross section ofa target can be obtained by evaluating a physical optics integral over the illuminated portionAppendix A. Design Curves for Top Hat Reflectors 144of its surface. In the case of scattering by a cylinder with radius a andheight 2c’, the physicaloptics integral can be expressed as the product of axial and circumferential components= _ji i, (A.5)where the axial component I is given by‘C,I,= Jejkz((i_))dz , (A.6)—C,— 2sin[kc’( (i—A 7Ckc’(2. ( —, ( .and the circumferential component I,!, is given byI,, = I (ñ i) d , (A.8)J—ir/2Although (A.8) will yield an exact solution, the resulting expression contains special functionsand is difficult to handle. It is generally more convenient to approximate IqS by the method ofstationary phase where the stationary phase point is the specular line shown in Figure A.4 andthe stationary phase approximation is of the form [6]2ii-1/2I= Jg(q) e)dq g(0)e’° e’’,(A.9)whereg(q) = (ñ.î),(A.1O)f()= ka (n. (i—(A.11)andg(40)=(n0 i) (A.12)f(,)= ka (ñ0 (1—, (A.13)= ka (n0.(1— ..)). (A.14)Substituting (A.12) — (A.14) into (A.9) gives1/2(n)jka(iio.(f-))[a (n0.—. (A.15)Appendix A. Design Curves for Top Hat Reflectors 145An expression for the bistatic scattering cross section of the cylinder is obtained by substituting (A.15) and (A.7) into (A.5) which gives= —j2c’[2ka]1/2sin kC(Z .(z s)). i)ejka(no())(A.16)n.(i—s) kc (z.(z—s))In the case of forward scattering, the vectors €,,ñ, and are related as follows:(1—.)= —2sin8 ñ,, , (A.17)• (1—.)= —2sin0 , (A.18). (f—= 0, (A.19)= —sin0. (A.20)An expression for the forward scattering cross section of the cylinder is obtained by substituting(A.17) — (A.20) into (A.16) which gives= —j2c’ v”ka sin8&32kasinO• (A.21)Since (A.18) is negative, a phase factorejlr/2is introduced when its square root is extracted.Taking the amplitude of (A.21) and squaring the result gives= acI2sin6.(A.22)Finally, an approximate expression for the scattering cross section of a top-hat reflector isobtained by substituting (A.1) into (A.22) to giveac2sinO,0<87,a(0)= 8ir cos2(A.23)—i- a(b — a)2-y8 < ir/2,where a, b, c, and 8 are defined in Figure A.1 and the parameter 7 is given by (A.2). If theground plane of the reflector is of very large extent, i.e., b>> c, the angle-y—* r/2 and (A.23)reduces toa(8)= ac2sin8,0< 8< ir/2. (A.24)Appendix A. Design Curves for Top Hat Reflectors 146A.3 Design CurvesAngle of Maximum ResponseThe contribution of double-bounce reflections to the response of a top hat reflector will reach itsmaximum value when the angle of incidence is such that any ray incident on the cylinder willintercept the annulus and vice versa. From Figure A.2, this will occur when c = (b — a) cot 6so the angle of maximum response6mis given by6max= tan1(b a)= tan’, (A.25)where is defined as the ratio of the width of the annulus to the height of the cylinder. A graphdepicting8mas a function of is presented in Figure A.4. If the height of the cylinder andthe width of the base are equal, the angle of maximum response is 45 degrees. Doubling theratio will lower the angle of maximum response to approximately 65 degrees while halving itwill raise the angle to approximately 25 degrees.9075600)45301500.1 10.0Figure A.4: Angle of maximum response of a top hat reflector vs. , the ratio of the annuluswidth to the cylinder height.1.0EAppendix A. Design Curves for Top Hat Reflectors147Maximum Scattering Cross SectionSubstitution of (A.25) into (A.23) gives an expression for the maximum scattering crosssectionof a top hat reflector8ir a(b — a)c2Umax= A (b — a)2’(A.26)8ir2_____= —acwhere (sin8max)is given bySflumax=+. (A.27)Consider a top hat reflector with a cylindrical component of fixed radiusa0 and height c0 whichpresents a maximum scattering cross section of o for = 1. As the ratio is altered toobtain the desired angle of maximum response, the maximum scattering cross section of thereflector will vary as shown in Figure A.5. The response will drop off rapidly for < 1 andwillasymptotically approach a value of crc, for > 1.The radius and height of the cylinder can be scaled to compensate for thisvariation andso yield a family of top hat reflectors with different angles of maximum response but thesamemaximum scattering cross section. With the cylinder height c fixed, the scale factor for thecylinder radius a is given by______a/1+2(A.28)a 2With the cylinder radius a fixed, the scale factor for the cylinder height c (and the correspondingannulus width b — a) is given by_____c11+c2—=/. (A.29)c0 2In either case, the corresponding value of b is simply given byb=a+c. (A.30)A graph depicting the scale factors (A.28) and (A.29) as functions of the ratiois presentedin Figure A.6.Appendix A. Design Curves for Top Hat Reflectors 148I04-’0201.50.51.0E0.00.1 10.0Figure A.5: Maximum scattering cross section of a top hat reflector vs. , the ratio of theannulus width to the cylinder height for fixed values of a and c.5.04.03.02.01.00.00.110.0Figure A.6: Scale factors for the cylinder radius a and height c of a top hat reflector vs. , theratio of the annulus width to the cylinder height.1.0CAppendix A. Design Curves for TopHat Reflectors 149Half Power BeamwidthThe half-power elevation beamwidth of a tophat reflector is determined by finding the anglesof incidence above and below theangle of maximum response,0mfor which the scatteringcross section fails to half its maximum valueand taking their difference. From (A.26),A (I.2oir a —a1c(A31)2 — A Vc2+(b_a)2Equating (A.23a) and (A.31) then solving for 6 gives theupper half-power angle,= sin’(2i+2). (A.32)The lower half-power angle3dBis determined in a similar fashion. Equating (A.23b)and(A.28) then cancelling common factorsgives a quadratic equation in sin 6,1 1—sin20A2/1+2 = sineRearranging terms gives an expression instandard form,sin26+(2V-+2)sin0— 1 = 0, (A.34)with roots given by______• _B+v’B2+4sin6= 2(A.35)whereB=1• (A.36)2/1+.2Since is always positive, so too are Band the solution to (A.34) given by• -B+VB2+4sin 0= 2• (A.37)Solving for 6 gives the lower half-power angle,°3dB=2J• (A.38)Appendix A. Design Curves for Top Hat Reflectors 1507560,453015Finally, taking the difference between (A.38) and (A.32) gives the total half-power elevationbeamwidth of the reflectorf—B+/B2+4’\ .I03d]3= sin2— sin2 s/i.+2(A.39)Care must be taken in applying (A.39) since the main response lobe of a top hat reflector isnot symmetric about the angle of maximum response, particularly if> 1. In addition to O3,the total half-power beamwidth, it is convenient to define the quantities0U3dB = 9u3dB —and0e3dB = —as the half-power beamwidths above and below the angle of maximumresponse, respectively. All three quantities are graphed in Figure A.7 as a function of the ratio. Although both the total and upper half-power beamwidths increase with , the lower half-power beamwidth reaches its maximum value when 0.55. The upper and lower half-powerbeamwidths are equal and the top hat presents a symmetrical elevation response pattern when0.51.9000.1 1.0 10.0Figure A.7: Half-power elevation beamwidth of a top hat reflector vs. , the ratio of the annuluswidth to the cylinder height.Appendix A. Design Curves for Top Hat Reflectors151where B’ is given by0)w0(10—0.16UIdB= sin1/1+2)B+%/B2+4”\84dB= Sfl12Similar relations may be derived for the 1 dB beamwidth of the elevation response. Theupper and lower 1 dB angles are given by(A.40)(A.41)111—0.1B’— “ (A42)and the total 1 dB elevation beamwidth is given by0u1dBor-1(—B’ +i/B’2+4-1(10-0.1=2}— sin+2}. (A.43)As in the previous case, it is convenient to define the quantities0uldB = 8uldB — °mand= 9max — 0eldnas the 1 dB beamwidths above and below the angle of maximum response,respectively. All three quantities are graphed in Figure A.8 as a function of the ratio .907545301500.1 1.0 10.0Figure A.8: 1 dB elevation beamwidth of a top hat reflector vs. , the ratio of the annuluswidth to the cylinder height.Appendix A. Design Curves for Top Hat Reflectors 152Since the elevation response pattern of a top hat reflector is generally asymmetricabout thedirection of maximum response, it may be useful to instead specify the median angle of responsegiven by8m=(9,+8)/2. The angle of maximum response and angles of median responsefor1 and 3 dB elevation beamwidths are graphed as a function of theratio in Figure A.9.9075I30w1500.1 10.0Figure A.9: Angle of maximum response and angles of median response for 1and 3 dB elevationbeamwidths vs. , the ratio of the annulus width to the cylinder height.A.4 Design ExampleConsider a requirement for three even-bounce targets to be used in thecalibration of an airborneradar system which operates at a frequency of 10 GHz. The three targets must each presenta maximum scattering cross section of 30 dBsm at elevation angles of 25, 45, and 65 degrees,respectively, over a wide range in azimuth. Top hat reflectorswhich meet these requirementscan be designed using the results presented in theprevious section.The ratio which gives the desired angle of maximum response can be determined usingeither (A.25) or the graph in Figure A.4. It is found that top hat reflectors whichpresent aratio of annulus width to cylinder height of 0.5, 1.0, and 2.0, respectively, willpresent their60451.0Appendix A. Design Curves for Top Hat Reflectors 153maximum response at the specified angles. This information is sufficient to obtain the 1 and 3dB beamwidths of the three reflectors using either (A.39) and (A.43) or the graphs in FiguresA.7 and A.8.The reflector dimensions required to realize a maximum scattering cross section of 30 dBsmat 10 GHz are determined using (A.26). Once the radar wavelength A and the specified angleof maximum response have been specified, the maximum scattering cross section depends onlyon the product of the cylinder radius a and the square of the cylinder height c. The design ofa top hat reflector which presents its maximum response at 45 degrees is considered first. It isconvenient to choose the diameter and height of the cylinder to be approximately equal. From(A.26), this gives a cylinder radius a of 0.75 m and a height c of 1.5 m. Since the ratio is1, the annulus width (b — a) is also 1.5 m which gives a total radius b of 2.25 m. The resultsare then scaled to obtain suitable dimensions for the other two reflectors. Assuming that thecylinder radius is fixed, the dimensions of the other two reflectors can be determined using thescale factors presented in (A.28) and (A.29) and graphed in Figure A.6.The results are summarized in Table A.1. In Figure A.10, the three reflectors are drawnto scale for comparison. In each case, the direction of the maximum response is indicated bya vector. In Figure A.11, the elevation response patterns of the three reflectors are predictedusing Blejer’s [4] formulation which accounts for double-bounce reflections from the interior ofthe reflector and single bounce reflections from cylinder, cylinder cap, and annular base.Appendix A. Design Curves for Top Hat ReflectorsTable A.1: Response Characteristics of Selected Top Hat Reflectors atf= 10 GHz154Fi ureDimensions Response Characteristicsa b Cumax 9max OldBA.10(a) 0.75 m 1.69 m 1.89 m0.5 30 dBsm 25° 27°100A.10(b) 0.75 m 2.25 m 1.50 m 1.030 dBsm 45° 36°15°A.10(c) 0.75 m 3.75 m 1.33 m 2.030 dBsm 65°45021°(c)Figure A.10: Relative size of selected top hat reflectors which present the same maximumscattering cross section. The direction of maximum response is indicated by a vector. (a) 0.5.(b) = 1.0. (c) = 2.0.43NTb(a) (b)TI.. jAppendix A. Design Curves for Top Hat Reflectors 1550•(a)0•(b)(c)a (dB)Figure A.11: Elevation response patterns of selected top hat reflectors which present the samemaximum scattering cross section. (a) = 0.5. (b) = 1.0. (c) = 2.0.-10 -6 -3-10-10 -6-3-100•-10 -6-3-10References[1] E.L. Johansen, “Top hat reflectors cap radar calibration,” Microwaves, vol. 20, no. 12,pp.65—66, Dec. 1981.[2] C.H. Currie and N.C. Currie, “MMW Reflectivity Measurement Techniques,” in Principlesand Applications of Millimetre-Wave Radar, N. C. Currie and C.E. Brown, Eds. Norwood,MA: Artech House, 1987,pp.767—774.[3] R. N. Trebits “Radar Cross Section,” in Radar Reflectivity Measurement: Techniques andApplications, N. C. Currie, Ed. Norwood, MA: Artech House, 1989,pp.44—48.[4] D. Blejer, “Physical optics polarization scattering matrix for a top hat reflector,” IEEETrans. Antennas Propagat., vol. AP-39, pp. 857—859, June 1991.[5] E.F. Knott, J.F. Shaeffer, and M.T. Tuley, Radar Cross Section: Its Prediction, Measurement, and Reduction. Norwood, MA: Artech House, 1985,pp.119—130.[6] N. Bleistein and R.A. Handelsman, Asymptotic Expansions of Integrals. New York: bit,Rinehart, and Winston, 1975,pp.219—220.156Appendix BSCATTERING BY A CONDUCTING GRATINGWITH RECTANGULAR GROOVESB.1 IntroductionIn this appendix, the problem of scattering by a conducting grating with rectangular groovesis solved for the case in which the plane of incidence is perpendicular to the direction of thegrooves [1], [2]. In section B.2, analytical solutions to the scattering problem are derived for thecases of TM- and TE-polarized incidence by mode-matching across the boundary between thefree space and groove regions. Since the grating is periodic in y, the fields in adjacent unit cellsdiffer only by the phase factor exp(jkdsin4) and it is sufficient to consider the fields in a singleunit cell such as the one shown in Figure B.l. Symbols are defined in Table B.l. The fields inthe free space region (x > 0) are represented by an infinite sum of propagating plane waves andevansecent waves while the fields in the groove region (x < 0) are represented by an infinitesum of propagating and evanescent parallel-plate waveguide modes. A doubly infinite set oflinear equations in the complex coefficients of either the free space or groove modes is obtainedby applying the condition of continuity of the tangential electric and magnetic fields over theplanar junction between the free space and groove regions and the property of orthogonalitybetween the normal modes in each set. The set must be truncated or partitioned before it canbe solved numerically. In section B.3, procedures for verifying the correctness and accuracy ofnumerical results obtained by mode-matching are described and the problem of determiningthe minimum number of modes required to accurately represent the fields in the free space andgroove regions and their optimum ratio is considered. In section B.4, an implementation of theanalytical solutions as a pair of subroutines coded in Fortran 77 is presented.157Appendix B. Scattering by a ConductingGrating158yhFigure B.1: A unit cell of a conductinggrating with rectangulargrooves.Table B.1: Definition of SymbolsSymbolDefinitiond period of the gratinga width of the grooveh height of the grooveaid aspect ratio of thegratingA wavelengthçb angleof incidenceqangle of reflection for themth free space mode ordiffracted orderm index for free spacemodesn index for groove modesM number offree space modes in the truncatedset of equationsN number ofgroove modes in the truncated set ofequationsk propagation constantin free space(=2ir/A)k propagationconstant of the nth groove modefmorthonormal basis function for TMand TE free space modesgorthonormal basis function for TMgroove modesg orthonormalbasis function for TE groove modesA2 complexamplitude of the incident TMfree space model3 complex amplitudeof the incident TE free spacemodeAmcomplex amplitude of the mthscattered TM free space modeC complex amplitudeof the nth TM groove mode‘3mcomplex amplitude of the mth scatteredTE free space modeD complex amplitude ofthe nth TE groove modei, characteristic impedance offree space(= f)Fmcomplex reflection coefficient of themth free space modexK aK d >1Appendix B. Scattering by a Conducting Grating159B.2 AnalysisB.2.1 TM PolarizationFree Space RegionIf a TM-polarized plane wave is incident on a gratingof infinite extent such that the planeof incidence is perpendicular to the direction of the grooves, the tangentialcomponents of theelectric and magnetic fields in the free space region are given byH(x, y) = A1ei0stxfo+ Ame_c0mxfm,(B.1)m -_E(x, y) = —A1 cos jejkcos+ Am o cosme_jkCOS mXf, (B.2)m=- oowhereE(x,y)=(B.3)— 1jksinmyB4Jm — e ,k = , (B.5)sinqm = sinq!’j+m, (B.6)cosq = — sii2Qm, I SiflmI< 1, (B.7)= _jsin2m— 1, sinmI> 1.The first term in each of (B.1) and (B.2) represents the incident field while the second termrepresents the infinite sum of propagating and non-propagating waves which comprisethe scattered field. The functionfmis an orthonormal basis function for both TM- and TE-polarizedmodes in the free space region. Equation (B.6) is usually referred toas the grating equation.The remaining quantities and symbols are defined in Table B.1 and Figure B.1.Appendix B. Scattering by a Conducting Grating160Groove RegionThe tangential components of the TM-polarized fields in thegroove region are given byH(x,y)=(B.8)E(x,y)=(B.9)where the functiongnis an orthonormal basis function for TM-polarized modes in thegrooveregion given by=i/cos(nir(x + a/2))(B.10)and,.k =/k2— (!E)2, < k , (B.11)_j(!)2_k2,II>k,e = 1, n=0,(B.12)= 2, n=1,2,3Boundary ConditionsThe y-z plane is the interface between the free space and groove regions. The boundary conditions for the tangential components of the TM-polarized field are given byH(0,y) = H(0,y),I’a/2,(B.13)Hz(O,Y)= —K11(0,y), a/2lvi< d/2,andE(0,y) = E(0,y),lvia/2,(B.14)E(0+,y) = 0, a/2< vi<d/2,where K represents the surface current density along the top surface of the corrugation.Appendix B. Scattering by a Conducting Grating161Mode MatchingThe modal expansions of the tangential magnetic field, (B.1) and(B.8), and the tangentialelectric field, (B.2) and (B.9), areevaluated at x = 0 then substituted into the boundaryconditions (B.13) and (B.14), respectively. Continuity of thetangential magnetic field acrossthe boundary givesAmfm = IyI<a/2,(B.15)m=— n=0while continuity of the tangential electric field across the boundary gives—2Acosqfo+ Acosqf IyI<a/2,(B.16)m=—oc n=0= 0, a/2<IyI<d/2.These expressions have been simplified by combining thecontribution of the incident and specularly reflected components of the fields in the free spaceregion at z = 0 into a single term andredefining A0 asAO=Ai+ARef,(B.17)where .A is the complex amplitude of the incident field and ARe1 isthe complex amplitude ofthe specularly reflected component of the scatteredfield.Inner ProductsThe inner products of the free space and groove modebasis functions are defined asd/2(fm,!’) = Ifmf,ndX,(B.18)J—d/2fa/2(gn,gn’)= Jg,-, dx,(B.19)— a/2a/2(fm,gn)= jfmg,dX,(B.20)—a/2a/2= Jg.fdx. (B.21)Appendix B. Scattering by a Conducting Grating 162Since the basis functions are orthonormal,(fm,fi) = 0, mm’, (B.22)= 1, m=m’, (B.23)(g,gn’)= 0 , ii ii’, (B.24)= 1, n=n’. (B.25)The inner products(fm,gn)and (g,,f> are given byIT 2ksinqrn (kasinq)m.71/—2 2cosj j , n=1,3,5,...v ad (ksincm) — (n2r/a) 2 j(fm,gn)= (B.26)fT 2ksinqm (kasin”Vad (ksinq5m)2— (nir/a)2sin2)n = 0,2,4,...andI(fm,gn), n=1,3,5,...=(B.27)((fm,gn), n=0,2,4,...since (f,g)= (g,f)*,by definition [4, p. 26].Substitution and SolutionThe inner product of (B.15) with g yields an expression for the groove mode coefficients C interms of the free space mode coefficientsAm,Am(fm,gn) = Cn, (B.28)m= -00while the inner product of (B.16) withfmyields—2Acosq+ AmCO5bm= Cntan(knh)(gn,fm) . (B.29)The result is then rearranged to give an expression for the free space mode coefficientAmflterms of the groove mode coefficients C,Am =k,tan(kh)n, fm> +(B.30)n=Ocos4mAppendix B. Scattering by a Conducting Grating163A direct solution for the coefficients of the free space modes can be obtained by substituting(B.28) into (B.30) to giveAm,6mmI=>(B.31)— — —j cosqS,m’——oofl_0__then rearranging the result to yield an infinite set of linearsimultaneous equations in thecoefficients of the free space modesAm,:Am’(fm’,gn)(g,f)tan(kh) —S, = —2A . (B.32)— —j cosq5m’——-oo n—OIt is convenient to express (B.32) as a matrix equation inthe form([B][C]+[G])A=F, (B.33)where the elements of [B], [C],[G], A, and F are given by13 — , :k tan(kh)— \gn,JmJ-—- , .34jkCOSmCnm’= (fm’,gn), (B.35)——Ami = Am’, (B.37)Fm = —2A. (B.38)The number of free space and groove modes included in the solutionmust be reduced froman infinite number to M and N, respectively, before a numerical solution to(B.33) can beobtained. The resulting truncated set of equations is of orderM.Alternatively, the coefficients of the free spacemodes can be determined in an indirect wayby first solving for the coefficients of the groove modes then multiplyingthe result by the matrixwhich relates the two. Since this truncatedset of equations is the order N, the time requiredto compute a solution can be considerably reduced comparedto the direct approach if the ratioof groove modes to free space modes IV/M in the truncated set is muchless than unity. First,Appendix B. Scattering by a Conducting Grating164(B.30) is substituted into (B.28) toyield((g’,f) + 2A6)(fm,gn) =(B.39)m=—-x n’=Omn’Othen the result is rearranged to yield an infinite set of linear simultaneousequations for thecoefficients of the groove modes Ci:Cstan(kh) _icot(knih)) = —2A(f0,g). (B.40)It is convenient to express (B.40) as a matrix equation in theform([C’][B’]+[G’j)X = F’ , (B.41)where the elements of [C’], [B’], [G’], X, and F’ are given byC,cm = (fm,gn>, (B.42)= jkcosqrn(g’,f) ,(B.43)= —6,k1cotk1h, (B.44)X1 = Ctankh, (B.45)F, = —2A (f0,g) . (B.46)On solution, the vector X contains the groove modecoefficients C multiplied by the factortan(kh). From (B.30), the free space mode coefficients Am aredetermined by applying therelationA—[C’]X+2A6.(B.47)From (B.17), the specularly reflected component of the scatteredfield is given byARJ = A0 — A,(B.48)so the specular reflection coefficient for a TM-polarizedincident wave is given bypEE— ARf(B 490— A ‘A0= --;+1.Appendix B. Scattering by a Conducting Grating165B.2.2 TE PolarizationFree Space RegionIf a TE-polarized plane wave is incident on a grating ofinfinite extent such that the planeof incidence perpendicular to the direction of the grooves, thetangential components of theelectric and magnetic fields in the free space region aregiven byE2(x,y)= !3je3kc0txfo+ Bme_c0mxfm,(B.5o)m=-_H(x,y)B2COScSkcostf —BmOS4me_ikc0mxfm, (B.51)m=—oo1owherefm,k,‘q0,sinq5m, and cos4mare defined in (B.3) through (B.7) andj OE(x,y).5“lo ‘‘YGroove RegionThe tangential components of the TE-polarized fields inthe groove region are given byE(x,y)=(B.53)H(x, y) =_jD-!_cos((B.54)where k is given by (B.11) and g is the orthonormalbasis function for TE-polarized modesin the groove region and is given by(B.55)Boundary ConditionsThe y-z plane is the interface between the freespace and groove regions. The boundary conditions for the tangential components of the TE-polarized field aregiven byH(O,y) = H(O,y),IIa/2(B.56)H(O,y) = K(O,y), a/2IId/2,Appendix B. Scattering by a Conducting Grating166where K represents the surface current density along the top surfaceof the corrugation andE(O,y) = E(O,y),yIa/2,(B57)E(O,y) = 0,a/2I’Id/2.Mode MatchingThe modal expansions of the tangential electric field, (B.50) and (B.53),and the tangentialelectric field, (B.51) and (B.54), are evaluated at x = 0 thensubstituted into the boundaryconditions (B.56) and (B.57). Continuity of the magnetic field acrossthe boundary gives2I3 cosçjf0 —>2‘3mCO5mfrn= >2Dcot(kh)g,lxi<a/2, (B.58)m=— n=1while continuity of the electric field givesm3mfm=lxI<a/2,(B.59)= 0, a/d<ixi<d/2.These expressions have been simplified by combining thecontribution of the incident and specularly reflected components of the fields in the free spaceregion into a single term and redefiningas13i +I3Ref,(B.60)where B is the complex amplitude of the incident fieldand!3Refis the complex amplitude ofthe specularly reflected component of the scattered field.Inner ProductsThe inner products of the TE mode basis functionsare defined in a similar way to those forthe TM modes with(fm,g)and (g,,fm) given by2nir/a (kasinqm”iV(ksinm)2_(n/a)2co2)n=1,3,5,...(fm,g)=(B.61)•2n’r/a . (kasinqrn3Vad (ksinbm)2— (nir/a)2 2 )‘ = 2,4,6,...Appendix B. Scattering by a Conducting Grating 167and(fm,g), n=1,3,5,...(g,fm)={, n=2,4,6,...(B.62)Substitution and SolutionThe inner product of (B.59) with g yields an expression for the free space mode coefficients‘3m,Dn(gj,fm)Bm,(B.63)while the inner product of (B.58) withfmyields00Dkcot(kh) 2Bjkcos (f0,g) — I3mjkcosqSrn(fm,g). (B.64)m=— 00The result can be rearranged to yield an expression for the groove mode coefficients D,= 2’3jk cos00jk cos çbrn(f0,gj—kcot(kh)m=-mkcot(kh)(fm,gj. (B.65)A direct solution for the coefficients of the free space modes can be obtained by substituting(B.65) into (B.63) to yield00/00=(jk cosm cos‘fo,g)) (g,frn), (B.66)n=1 \m’=_00kcot(kh)(f,,g)—2k cot(kh)then rearranged to give00 00 00jk cosm______________jk cos— ‘3m‘ = ‘3m’kh)(fm’,gn)fm)_2t3iJt(k h)(f0,g) (g,f)k cot(m’=—oo m’=—c,o 72=1(B.67)The final result is an infinite set of linear simultaneous equations for the coefficientsof the freespace modes‘3m’00 00cosmlCocos‘3m’Ecot(k72h)(fm’,gn) (g,fm) +, =2Bm’=—ook72cot(kh)(fo,g72)(g,fm)(B.68)Appendix B. Scattering by a Conducting Grating 168It is convenient to express (B.68) in the form([B][C] + [G]) X = [B] F, (B.69)where the elements of [B], [C], [G], X, and F are given by:Bmn=(g,f,) , (B.70)— jkcosq5/Bnm’— kcot(kh)fm’,gni, .71—— Vml ,Xmi = Bm’, (B.73)= 2Bi)(f0,g) (B.74)On solution, the vector X contains the free space mode coefficients13mrnAs in the case of a TM-polarized incident wave, it may be advantageous to solve for the coefficients of the groove modes then multiply the result by a matrix which relates the coefficientsof the groove modes to those of the free space modes if N<< M. First, (B.63) is substitutedinto (B.64) to yieldDkcot(kh) = 2B jkcos(f0,g) —(n,(ifm)) jkcosm(fm,g)(B.75)m=— ,i=1which is rearranged to give(Dni(nifm))ikcosbm(fm,g)+ ‘Dikicot(kh) = 2Bjkcosq(f0,g)m=—oon’=l n’=l(B.76)The final result is an infinite set of linear simultaneous equations for the coefficients of thegroove modes Dpi:>jkcos(gn’,fm) (fm,g) +6,kicot(k;h) = 2Bjkcos1(f0,g) . (B.77)n’1 m—It is convenient to express (B.77) in the form([B’J[C’]+ [G’])X’ = F’ , (B.78)Appendix B. Scattering by a ConductingGrating 169where the elements of [B’], [C’], [G’], X’, andF’ are given by:Bm= (fm,g), (B.79)= jkcosq5m(g-1,fm), (B.80)G1 = 6,kcot(kh) , (B.81)X, = ,(B.82)F, = 2Bjkcosc(f0,gj .(B.83)On solution, the vector X contains the groovemode coefficients D. From (B.63), the freespace mode coefficients13mare given by1[B’]x, (B.84)jk cosFrom (B.60), the complex amplitude of the specularlyreflected component of the scattered fieldis given byl3Ref 13o — B,(B.85)so the specular reflection coefficient for a TE-polarizedincident wave is given by:pHH(B.86)L3113iB.3 Verification of Numerical ResultsThe analysis presented in section B.2 is exact but truncationof the mode-matching matrix anduse of a finite word length during computationintroduce unavoidable approximations. Furthermore, errors in analysis and coding sometimes occur. It istherefore desirable to implementprocedures to assess the correctness and accuracy of thenumerical solutions which are obtained.Three numerical checks which are applicable to thisproblem include conservation of energy,reciprocity, and convergence [5].Appendix B. Scattering by a Conducting Grating 170Conservation of Energy and ReciprocityThe physical plausibility of numerical results obtained by mode-matching may be verified bytesting for conservation of energy and/or reciprocity among the propagating modes but neitheris a sufficient condition for a valid solution. Such tests serve mainly to indicate problemsintroduced by errors in coding or excessive round-off during computation [5], [6]. In the contextof scattering by a reflection grating, conservation of energy among the free space modes maybe expressed asIAiI2cos4i Ii4mt2C0m, (B.87)where A and 4j are the complex amplitude and the direction of the incident wave, Am and4’mare the complex amplitude and the direction of the mth free space mode, and the summationincludes only those free space modes which are propagating. Reciprocity is a statement thatthe response of a system will be unchanged when the source and the observer are interchangedand can be expressed as= fm(bm),(B.88)whereFm(q5)is the complex reflection coefficient of the mth diffracted order in response to awave incident at angle q. In the case of the specularly reflected free space mode, çbm = qj and(B.88) reduces to=l(—q,). (B.89)Testing for conservation of energy incurs little computational overhead but testing for reciprocity effectively doubles the length of time required to compute a solution and is rarely doneso routinely.ConvergenceIt is necessary to determine whether a sufficient number of modes have been used to approximate the fields in the free space and groove regions. Although the solution will invariablyconverge to an essentially constant result as the number of modes used in the field expansionAppendix B. Scattering by a Conducting Grating 171are systematically and gradually increased, the manner in which the solution converges and thevalue of the final result will depend on both the number of modes used in each region and theirratio, N/M. The phenomenon, which has been dubbed relative convergence, was originallyidentified in connection with the formulation of the boundary value problem associated witha bifurcated waveguide [7]. It has since been shown to occur in a variety of other problemsincluding diffraction by a strip grating on a dielectric slab [8], [9]. In these cases, Mittra andothers have shown that there usually exists a unique choice for the so-called partitioning ratiowhich will yield the correct solution. Although results obtained using an arbitrary partioningratio and a large number of modes are usually accurate to within engineering tolerances, use ofthe optimum ratio provides both greater accuracy and greater computational efficiency [10].In Figures B.3 through B.6, the phase difference between the TE and TM specular reflectioncoefficients of selected gratings are plotted as a function of the number of groove modes N usedin the solution for fixed values of the corresponding number of free space modes M. In allcases, the gratings have a period of 0.3333) and are illuminated perpendicular to the directionof the grooves at an angle of 45 degrees to the normal. Figures B.3 and B.4 present the resultsfor a pair of twist polarizers with aspect ratios aid of 0.5000 and 0.9999, respectively, whileFigures B.5 and B.6 present the results for a corresponding pair of circular polarizers. Thephase difference converges monotonically until the ratio of groove modes to free space modesreaches a value equal to the aspect ratio of the grating. At that point, the slope of the curveincreases abruptly and the phase difference begins to converge to a new value which depends onthe number of free space modes. In the limit as the number of free space modes becomes verylarge, the difference between the final result and the result obtained at the point of inflectionvanishes which suggests that the result obtained at the point of inflection corresponds to thecorrect solution. This is similar to behaviour reported in connection with other formulations inwhich relative convergence has been observed and implies that the optimum number of groovemodes is given by truncating the product of the aspect ratio of the grating and the number offree space modes to be used in the solution, i.e., N0 = Trunc(a/d. M).Appendix B. Scattering by a Conducting Grating 172a,0Ca,1811790 10 20 30 40 50 60Number of groove modesFigure B.2: Convergence of the phase difference between the TE and TM specular reflectioncoefficients with the number of groove modes for a reflection twist polarizer with d = 0.3333),aid = 0.5000, and h 0.2302A illuminated by a plane wave incident at = 45 degrees.181180179Figure B.3: Convergence of the phase difference between the TE and TM specular reflectioncoefficients with the number of groove modes for a reflection twist polarizer with d = 0.3333 A,aid = 0.9999, and h = 0.3172A illuminated by a plane wave incident atq5= 45 degrees.0 10 20 30 40 50 60Number of groove modesAppendix B. Scattering by a Conducting Grating173C091890 1020 3040 5060Number of groovemodesFigure B.4: Convergence of the phase difference between the TE and TM specular reflectioncoefficients with the number of groove modes for a reflection circular polarizer with d = 0.3333A,aid = 0.5000, and h = 0.1466A illuminated by a plane wave incident at = 45 degrees.9190890 10 2030 40 50 60Number of groove modesFigure B.5: Convergence of the phase difference between the TE and TM specular reflectioncoefficients with the number of groove modes for a reflection circular polarizer with d = 0.3333),aid 0.9999, and h = 0.1641A illuminated by plane wave incident at = 45 degrees.Appendix B. Scattering by a Conducting Grating174B.4 ImplementationSubroutines TMREFL and TEREFL implement the analytical solutions given by(B.41) and (B.78)and are presented in Listings B.1 and B.2. Given theperiod d, groove width a, and groovedepth h of the grating (in metres), the angle of incidence4(in degrees) and wavelength A(in metres) of the incident wave, and the number of free space and groove modes M andNto be used in the solution, they return the complex specular reflection coefficient P0for TM-and TE-polarized incident waves, respectively, and an error code to the calling program.Errorcode 1 indicates that the solution violates conservation of energy among the free spacemodeswhile error code 2 indicates that the mode-matchingmatrix passed to subroutine CDSOLN issingular. Error code 0 indicates that the solution has passed both these tests. If testsfor eitherreciprocity or convergence are required, they must be performed by the callingprogram.Because the solution makes extensive use of complex variables, the subroutines werecodedin Fortran 77. Although the use of language extensions was generally avoided in order to keepthe source code portable between different compilers and platforms, two extensions whicharesupported by virtually all modern Fortran 77 compilers were allowed. First, double precisioncomplex variables (type COMPLEX*16) and the corresponding intrinsic functions were used sothat all floating point calculations could be performed using double precision. Second, variablesof type DOUBLE PRECISION were identified as REAL*8 in type declaration statements for clarity.Otherwise, the code is fully compliant with the ANSI standard and adheres tothe generallyaccepted principles of programming style [11]—[13}.A chart depicting the hierarchy of subprograms which are called by subroutines TMREFL andTEREFL is presented in Figure B.6. Functions TMINNR and TEINNR are presented in ListingsB.3 and B.4. Given the indices of the free space and groove modes rn and n, the propagationconstant k, the grating period d, the groove width a, and a character variable 0 whichspecifiesthe desired sequence of factors, they return the value of the inner products (f,g) and (g,f)as given by (B.26/27) and (B.61/62), respectively. Functions KN, ODD, SEQ, E,and CDCOT arepresented in Listings B.5 through B.9. Complexfunction KN returns the propagation constantAppendix B. Scattering by a Conducting Grating175of the nth groove mode as given by (B.11). Logical function ODD returns true if its argumentis an odd integer. In conjunction with ODD and IF... THEN constructs in functionsTMINNRand TEINNR, real function SEQ is used to select the appropriatesign for the inner products asspecified by (B.27) and (B.55). Complex function CDCOT implements the complexcotangentfunction. The remaining subprograms called by subroutinesTMREFL and TEREFL were suppliedby University Computing Services and perform operationsnot directly supported by Fortran77 such as addition and subtraction of matrices, multiplicationof matrices with vectors, andsolution of systems of linear equations. They are presented in ListingB.10 and are documentedin the publication UBC Matrix [14].TNREFLTEREFLThIINNRKNTEINNRKNICDCOTCDCOTCDADDCDADDODDCDMATV ODDCDMATVSEQCDMULT SEQCDMULTECDSETCDSETCDSOLNCDSOLNFigure B.6: Hierarchy of subprograms called by subroutines TMREFL and TEREFL.The subroutines were compiled on a Sun 4/380 workstation (equipped with 32 Megabytesof RAM) under SunOS 4.1.1 (a variant of 4.2 BSD UNIX) using the standard Sun Fortrancompiler, f77, with all options set to their default values. The UNIX operatingsystem provides a variety of timing and profiling tools which can beused to assist in the evaluation andoptimization of Fortran programs [15]. The total length of time required to execute the twosubroutines, including both system and user time, was determined by runningthe shell utilitytime in conjunction with a test program which simply set the parameters tobe passed, calledsubroutines TMREFL and TEREFL in succession, then terminated.The results are presented inFigures B.7 and B.8 as a function of the number of free space and groove modes usedin thesolution.Appendix B. Scattering by a Conducting Grating 1760ES.0E80 10 20 30 40 506420Number of groove modesFigure B.7: Combined execution time of subroutines TMREFL and TEREFL on a Sun 4/380workstation vs. the number of groove modes used in the solution.0010 20 30 40 50Number of free space modesFigure B.8: Combined execution time of subroutines TMREFL and TEREFL on a Sun 4/380workstation vs. the number of free space modes used in the solution.Appendix B. Scattering by a Conducting Grating 177A call graph analysis of the code was conducted using the profiling utility gprof. Theresults are presented in Figure B.9 for numerical solutions using (a)25 groove modes and 51free space modes and (b) 50 groove modes and 51 free space modes. In each case, calls tosubroutine CDMtJLT (which calculates the productof matrices [B] and [C] from (B.41) and [B’]and [C’] from (B.78)) account for approximately two-thirds ofthe total execution time. Theloading increases quadratically with the number of groove modesand linearly with the numberof free space modes. In contrast, calls to subroutine CDSOLN (which solvesfor the amplitudeof the groove modes) account for only one-fifth of the total execution time where the loadingalso increases quadratically with the number of groove modes but is independentof the numberof free space modes. The remainder of the execution time is spent calculating the elementsof the various vectors and matrices in (B.41) and (B.78). Since the elements of matrices[B]and [C] (and [B’] and [C’]) are independent on the groove depth h, these results suggest thatexecution time could be reduced considerably by not recalculating their product if h is the onlyparameter which has changed since the previous call. Such a modification could be implementedby inserting a SAVE statement into each subroutine in order to preserve the values of the localvariables between calls and applying an appropriate branch on entry.Figure B.9: Execution profile of subroutines TMREFL and TEREFL for numerical solutions using(a) 25 groove modes and 51 free space modes and (b) 50 groove modes and 51 free space modes.(a) (b)Appendix B. Scattering by a Conducting Grating 178ListingB.1: SUBROUTINE TMREFL(D, A, H, PHI, L, MFS, NWG, REFL, ERR)SUBROUTINE TJIREFL( D, A, H, PHI, L, NFS, luG, REFL, ERR)* Given:*D - grating period (metres)*A — groove width (metres)*H - grating depth (metres)*PHI — angle of incidence (degrees)*L — wavelength (metres)* KFS— number of free space modes* lUG — numberof groove modes* Result:* REFL — complexreflection coefficient of zeroth free space mode* ERR— error return 0 = no errors,*1 = energr in free space modes not conserved,*2 = matrix to be inverted is singular.** Allfloating point calculations are done in double precisionPARAJIETER (MAXFS=51, MAXUG51)* KAXFS, MAXUG- maximum array dimensions*(see also SUBROUTINE TMIINR)EXTERNAL CDCOTINTEGER IFS, JUG, ERR, I, 3, P1, N, ZEROREAL*8 D, A, H, PHI, L, K, PRRAD, P1, S, SINPH, SUMCOMPLEX*16 REFL, B, C, •E, F, G, RCOMPLEX*16 CDCOT, CDET, COSPH, KI, ThuRDIMENSION KOLAXFS)DIMENSION COSPH(-MAXFS/2 :IIAXFS/2), SINPH(-MAXFS/2 :IIAXFS/2)DIMENSION B(MAXFS,MAXWG), C(MAXWG,MAXFS)DIMENSION GOIAXUG HAXUG), EOIAXUG,MAXUG)DIMENSION FOIAXUG5COMMON P1, SIIPHZERO = (MFS+1)/2P1 = 4D0*DATAN(1DO)** convert angle to radiansPHRAD = PHI * PI/180D0* calculate free spacepropagation constantK = 2D0*PI/L*tabulate values of sin[phi(m)] and cos[phi(m)) for later useDO 50, I = -MFS/2, IIFS/2S = DSIN(PHRAD) + DBLE(H)*L/DSINPH(M) = SIF ((S.S).LE.1D0) THENCOSPE(H) = DCMPLX(DSQRT(1DO-S*S), ODO)ELSECOSPH(M) = DCJIPLX(ODO, -DSQRT(S*S-1DO))EIDIFSO COITINUEAppendix B. Scattering by a Conducting Grating 179Listing B.1: (continued)* set up the matrices* matrixBDO 100, I = 1, IFSDO 100, 3 = 1, lUGI = I - ZERO1=3 -1B(I,i) = TIIIINR(M, I, K, D, A, ‘gf’) *$Kl(K,1,A)/(DCMPLX(ODO,K)*COSPH(M))100 COITIIUE* matrixCDO 200, 3 = 1, lUGDO 200, I = 1, IFS1=3 -111 = I - ZEROC(3,I) = TNIN1R(N, 1, K, D, A, ‘fg’)200 CONTIJUE** [C]*[B)—> [E]CALL CDNULT(C,B,E,NWG,IIFS,NWG,MAXUG,MAIFS,MIZWG)* vector FDO 300, 3 = 1, lUG13—110F(3) = -2.0*TNIIIR(M, 1, K, D, A, ‘fg’)300 COITIIUE* matrixG400 CALL CDSET(G, lUG, lUG, ?IAXVG, (ODO,000))DO 500, 3 1, lUGl3-1G(J,J) = —CDCOT(K1(K,1,A)*H)500 CONTINUE*• [E]+[G] —> [E]CALL CDADD(E,G,E,IUG,IUG,NAXUG,KAXUG,MAIWG)• [E]*[X] = [F]; solve for [I]; [X]—>[F]CALL CDSOL1(E ,F,1WG,IIAXVG ,CDET)IF (ABS(CDET).LT.1E-20) THENREFL = (000,000)ERR2RETURNEIDIF* [B]*[F] —> [R]CALL CDMATV(B ,F,R,MFS,NWG,JIAXFS)R(ZERO) = R(ZERO) + (100,000)*check for conservation of energy in propagating free space modes*(a necessary but not a sufficient condition for a valid solution)SUM = ODOERR = 0DO 1000, I = 1, IFSI = I - ZEROIF (ABS(SINPHOI)) .LE.1DO)$ SUM = SUM+ (CDABs(R(I))**2) COSPH(H)1000 CONTINUEIF ( ABS(SUI-DCOS(PHRAD)).GT.0.000IDO ) ERR = 1* calculation complete!REFL = R(ZERO)*(-1DO,000)RETURNENDAppendix B. Scattering by a Conducting Grating 180ListingB.2: SUBROUTINE TEREFL(D, A, H, PHI, L, MFS, NWG, REFL, ERR)SUBROUTINE TEREFL( D, A, H, PHI, L, PIFS, IVG, REFL, ERR)* Given:* D — grating period(metres)* A — groove width(metres)* H — grating depth(metres)*PHI — angle of incidence (degrees)* L — wavelength(metres)*MFS — number of free space modes• JUG — number of groove modes* Result:*REFL — complex reflection coefficient of zeroth free space mode* ERR — error return 0 = no errors,• 1 energy in free space modes not conserved,* 2 = matrix to be invertedis singular.** All floating point calculations are done in double precisionPARAMETER (MAXFS=51, MAXUG=51)* MAXFS, MAXUG - maximum array dimensions* (see also SUBROUTINE TEINIR)EXTERNAL CDCOTINTEGER MFS, JUG, ERR, I, J, P1, N, ZEROREAL*8 D, A, H, PHI, L, K, PERAD, P1, S, SINPH, SUMCOMPLEX*16 REFL, B, C, E, F, G, RCOMPLEX*16 CDCOT, CDET, COSPH, KR, TEIJJRDIMENSION R(MAXFS)DIMENSION COSPH(—MAXFS/2:JIAXFS/2), SIJPH(-MAXFS/2:MAIFS/2)DIMENSION BOIAXFS,MAIUG), C(MAXVG,MAXFS)DIMENSION G(KAXWG KAXUG), E(MAXUG,MAZUG)DIMENSION F(MAXWGSCOMMON P1, SINPHZERO = (MFS+1)/2P1 = 4D0*DATAN(1DO)** convert angle to radiansPURAD = PHI * PI/180D0* calculate free space propagation constantK = 2D0*PI/L* tabulate values of sin[phi(m)J and cos[phi(m)]for later useDO 50, P1 = —HFS/2, PIFS/2S DSIN(PHRAD) + DBLE(M)*L/DSIJPH(M) SIF ((S.S).LE.1DO) THENCOSPH(M) = DCMPLX(DSQRT(1DO-S*S), ODO)ELSECOSPHOI) = DCMPLX(ODO, —DSQRT(S*S-1DO))ENDIF50 CONTINUE*Appendix B. Scattering by a Conducting Grating 181Listing B.2: (Continued)* set up the matrices* matrixEDO 100, I = 1, MFSDO 100, J 1, lUGN = I - ZERO1=1B(I,3) = TEIIIROI, 1, K, D, A, ‘gf’) *$ DCMPLI(ODO,K)*COSPE(M)100 CONTIIUE*matrixCDO 200, 3 = 1, JUGDO 200, I = 1, MFS‘=3N = I - ZEROC(J,I) = TEIIIR(M, 1, K, D, A, ‘fg’)200 COITIJUE*• [C)*tBJ —> [E]CALL CDMULT(C,B,E,IVG,MFS,IUG,MAXUG,MAIFS,MAIWG)*vectorFDO 300, 3 1, lUG1=3M0F(3) = 2.0*TEIIIR(M, 1, K, D, A, ‘fg’) •$ DCMPLZ(ODO,K)* DCIIPLX(DCOS(PERAD))300 COITIJUE*matrix G400 CALL CDSET(G, lUG, JUG, NAXUG, (ODO,ODO))DO 500, 3 = 1, lUG‘=3G(J,J) = Kl(K,l,A)*CDCOT(Kl(R,l,A)*E)500 CONTINUE -* [E)+[G]—> [E)CALL CDADD(E,G,E,IVG,IWG,NAIUG,MAXWG,KAIWG)*[E]*[X] = [F]; solve for [I]; [X]—>[F]CALL CDSOLI(E ,F,IWG,NAXUG,CDET)IF (ABS(CDET).LT.1E-20) THENREFL = (000,OD0)ERR=2RETURNEIDIF* [B]*[F) -> [R]CALL CDMATV(B,F,R,MFS,IUG,NAXFS)R(ZERO) R(ZERO)-DCMPLI(ODO,K) * DCIIPLX(DCOS(PHRAD))* check for conservation of energylit propagating free space modesSUM = ODOERR = 0DO 1000, I = 1, MFSN = I — ZEROIF (ABS(SIIPH(M)) .LE.1DO)$ SUM = SUN + C CDABS( R(I)/( DCMPLI(ODO,K) •$ COSPE(M) ) ) )*.2 COSPH(M)1000 CONTINUEIF C ABS(SUM-DCOS(PHRAD)).GT.0.000100 ) ERR = 1* calculation complete!REFL R(ZERO)/(DCMPLI(ODO,K)*DCMPLZ(DCOS(PHRAD)))RETURNENDAppendix B. Scattering by a ConductingGrating 182ListingB.3: COMPLEX*16 FUNCTIONTMINNR(M, N, K, D, A, 0)COKPLEX*16 FUNCTIONTMIIIRO(, I, K, D, A, 0)PARAMETER (KAIFS=51)LOGICAL ODDINTEGER M, IREAL.8 K, D, A, DEN, E, NUN, P1, SEQ,SINPHCEARACTER*2 0DIMENSION SIJPH(—MAXFS/2 :MAIFS/2)COMMON P1, SIJPEDEN = K.SIIPR(M)*K*SINPH(M) - (DBLE(N)*pI/A).(DBLE(I)*pI/A)IF (ODD(I)) THENNUN = DCOS (A.KeSINPH (N) /2D0) *2D0.K*SINPH(K)TMINNR = DCMPLI(SEQ(O))*DCMPLI(ODO,—IUM/(DEI4’DSQRT(A*D)))ELSENUN = SIN(A*K.SIIPH(M)/2D0) *2D0*K*SIIPROI)TMIIIR = DCMPLI(DSQRT(E(N)/(A.D)) * IUJI/DEN,ODO)EIDIFRETURNENDListing B.4: COMPLEX*16 FUNCTIONTEINNR(M, N, K, D, A, 0)COMPLEI*16 FUNCTIONTEINNR(M, I, K, D, 1, 0)PARAMETER (JIAXFS=51)LOGICAL ODDINTEGER N, NREAL*8 K, D, A, DEN, NUN,P1, SEQ, SIIPECRARACTER*2 0DIMENSION SIIPH(-KAIFS/2 :KLIFS/2)COMMON P1, SliPliDEN K*SIIPH(M)*K*SIIPHOI) - (DBLE(I)*PI/A)*(DBLE(I)*PI/A)IF (ODD(I)) THENNUN = DCOS(A*K*SIIPH(M)/2D0)*2D0*DBLE(N).PI/ATEIIIR = DCMPLX(-(IUN/DEI)*DSQRT(2D0/(A*D)),ODO)ELSENUN = SIN(A*KSSINPH(M)/2D0)*2D0*DBLE(N)*PI/ATEIINR = DCWPLI(SEQ(O))*$ DCMPLI(ODO,-(JUMIDEI)*DSQRT(2D0/(A*D)))EIDIFRETURNENDAppendix B. Scattering by a Conducting Grating 183ListingB.5: COMPLEX*16 FUNCTION KN(M, N, K, D, A, 0)CONPLEI*16 FUNCTION KI(K,I,A)REAL*8 K, A, P1,QINTEGER NP1 = 4D0*DATAN(1DO)Q= DBLE(N)*PI/AIF (ABS(Q).LE.K) THENKN DCNPLX(DSQRT(K*K-Q.Q), 0Db)ELSEIN = DCMPLI(ODO, —DSQRT(Q*Q-K*K))ENDIFRETURIENDListing B.6: LOGICAL FUNCTION ODD(N)LOGICAL FUNCTION ODD(N)* true if N is odd, false if N is evenINTEGER NIF (KOD(N,2).EQ.o) THENODD = .FALSE.ELSEODD = .TRUE.ENDIFRETURNENDListing B.7: REAL*8 FUNCTION SEQ(O)REAL*8 FUNCTION SEQ(O)* if n is odd, <fm,gn> = <gn,fm> the order is unimportant*if n is even, <fm,gn> -<gn,fm> the order must be accounted for’CHAR.ACTER*2 0IF (O.EQ.’fg’) THENSEQ = 1DbELSESEQ = -1DbENDIFRETURNENDListing B.8: REAL*8 FUNCTION E(N)REAL*8 FUNCTION E(N)*calculates the Neumann numberINTEGER NIF (N.EQ.b) THENE = 1DbELSEE = 2DbENDIFRETURNENDAppendix B. Scattering by a Conducting Grating 184Listing B.9: COMPLEX*16 FUNCTION CDCOT(ARG)COMPLEX*16 FUICTIOJ CDCOT(LRG)* find the cotangent of a complex numberCWIPLEX*16 CDSIJ, CDCOS, ARGREALs8 DIMAGIF (DIIIAG(ARG).GT.100DO) TEEJCDCOT = (ODO,—1DO)ELSE IF (DIIIAG(ARG).LT.—100DO) TBEkCDCOT = (0DO,1DO)ELSE IF (ABs(cDsII(ARG)).LT.1E-1o) TRUICDCOT = (1E1O,O.O)ELSECDCOT = CDCOS(ARG)/CDSIJ(AKG)EWDIFRETURIENDAppendix B. Scattering by a Conducting Grating185Listing B.1O: UBC Complex Matrix SubroutinesSUBROUTIIE CDADD(1,B,C,M,1,1DIMA,1Dh1,1DIMC)*[C) <— [A)+[B] Ref: UBC Matrix, p. 15COMPLEX*16 A(NDIMA,*),B(IDIKB,*) ,C(NDIMC,*)DO 1 3=1,1DO 1 I=1,M1 C(I,J) = 1(1,3) + B(I,3)RETURNENDSUBROUTINE CDMATV(A,V,V,M,1,NDINI)**w <- [A].v Ref: USC Matrix, p. 17COMPLEX*16 A(NDIMA,e) ,V(*) ,w(*)DO 1 1=1,1U(I) = (0.00,0.00)DO 1 3=1,11 11(I) = w(I) + 1(1,3) • V(i)RETURNENDSUBROUTINE CDMULT(A,B,C,M,J,L,NDIMA,NDIMB,NDIMC)• [C] <— [A].[B] Ref: USC Matrix,p.18COMPLEX*16 1(IDIMA,*),B(IDIMB,.) ,C(NDIIIC..)DO 1 J=1,LDO 1 11,MC(I,J) = (0.DO,O.D0)DO 1 1=1,11 C(I,J) = C(I,J) + 1(I,K) * B(K,J)RETURNENDSUBROUTIIE CDSET(1,M,1,IDIMA,X)• [C] — [A]+[B] Ref: USC Matrix, p. 21COMPLEX*16 *(IDIMA,*),IDO 1 3=1,1DO 1 I=1,M1 1(1,3) = IRETURNENDSUBROUTINE CDSUB(A,B,C,M,N,NDIMA,NDIMB,NDIMC)*[C] <— [1]-[B] Ref: UBC Matrix, p. 22COMPLEI*16 A(NDIMA,*) ,B(IDIMB,.),C(IDIMC,*)DO 1 J11DO 1 11,M1 C(I,J) = 1(1,3) — B(I,J)RETURNENDSUBROUTINE CDVMAT(V,A,W,M,N,NDIMA)** a <- v.[A] Ref:USC Matrix, p. 24COMPLEX*16 1(1DIMI,s) ,V(*) ,V(e)DO 1 3=1,NW(J) = (0.D0.0.D0)DO 2 11,M2 W(J) = V(J) + v(I) • 1(1,3)1 CONTINUERETURNENDAppendix B. Scattering by a Conducting Grating 186Listing B.1O: (continued)SUBROUTINE CDSOLI (A, B, I ,MM ,DET)**[A).x = b; x—>b Ref: UBC Matrix, p. 55C This routine finds the solution of a system of equations AIB, and takesC paging into consideration; i.e., operations are done by columns,C A— matrix of coefficintsC B— vector of right hand sides; on exit array B wiil contain the solution.C I— order of matrix.C MM- first dimension of a.C DET- determinant of matrix is DET*10**JEXCOMMON /NATEXP/ JEXDIMENSION AOIfI,N),B(*)COMPLEX*16 A,B,DET,T,TB,DCMPLXREAL*8 RMAI , QDET , CDABSJEXODETDCMPLX(1 .D0,0.D0)C decompose a into upper triangular matrixDO 6 K1,IM=KRIIAX=CDABS (A (K , K) )IF(K.EQ.I) GO TO 14KP1K+1DO 1 IKP1,NIF(CDABS(A(I,K)).LE.RMAX) GO TO 1M1RAAXCDABS(A(I ,K))1 CONTINUE14 T=A(M,K)IF(RAAX.LT.1.D-20) GO TO 20IF(K.EQ.N) GO TO 5A (11, K) =A (K , K)A (K , K)TBB (K)BOO =B(K)B(K)=TBDO 2 I=KP1,NA(I X)=-A(I K)/T2 B(I5B(I)+AI,K)*TBDO 4 JKP1,NTBA(M,J)A(M,J)=A(X,J)A(K,J)TBIF(CDABS(TB).EQ.0.D0) GO TO 4DC 3 IKP1,N3 A(I,J)A(I,J)+A(I,K)*TB4 CONTINUEC now get determinant5 DETDET*TQDET=CDABS (DET)IF(QDET.LT.1.D15) GO TO 25DET=DET*1.D—15JEXJEX+1525 IF(QDET.GT.1.D-15) GO TO 30DETDET*1 .015JEXJEX-1 530 IF(M.NE.K) DET—DET6 CONTINUEC now do back substitutionIF(N.EQ.1) GO TO 9IM1N-1DO 8 KB1,NM1KM1N-KBKKM1+1B(K)B(K)/A(K,K)T=-B(K)DO 8 1=1 XMl8 B(I)B(I+A(I,K)*T9 B(1)=B(1)/A(1,1)RETURN20 DETDCMPLX(0.D0,0.D0)JEXORETURNENDReferences[1] A. Hessel, J. Shmoys, and D. Y. Tseng, “Bragg-angle blazing of diffraction gratings,”J. Opt. Soc. Am., vol. 65,pp.380—384, Apr. 1975.[2] J. W. Heath, Scattering by a Conducting Periodic Surface with a Rectangular GrooveProfile, M.A.Sc. Thesis, Univ. of British Columbia, Vancouver, B.C., 1977.[3] J.J. Van Zyl and F.T. Ulaby, “Scattering matrix representation for simple targets,” inRadar Polarimetry for Geoscience Applications. F.T. Ulaby and C. Elachi, Eds. Norwood,MA: Artech House, 1990,Pp.17—19.[4] J. W. Dettman, Mathematical Methods in Physics and Engineering. New York: McGraw-Hill, 1969,p.26.[5] C. P. Wu, “Methods of Solution,” in Theory and Analysis of Phased Array Antennas,N. Amitay, V. Galindo, and C. P. Wu, Eds. New York: Wiley, 1972,pp.75—119.[6] N. Amitay and V. Galindo, “On energy conservation and the method of moments in scattering problems,” IEEE Trans. Antennas Propagat., vol. AP-17,pp.747—751, Nov. 1969.[7] R. Mittra, “Relative convergence of the solution of a doubly infinite set of equations,” J.Res. Nat. Bur. Stand., vol. 67D, Mar-Apr. 1963,pp.245—254.[8] T. Itoh and R. Mittra, “Relative convergence phenomenon arising in the solution of diffraction from strip grating on a dielectric slab,” Proc. IEEE, vol. 59,pp.1363—1365, Sept. 1971.[9] R. Mittra, T. Itoh, and T-S Li, “Analytical and numerical studies of the relative convergence phenomenon arising in the solution of an integral equation,” IEEE Trans. MicrowaveTheory Tech., vol. MTT-20,pp.96—104, Feb. 1972.187References 188[101C.P. Wu, “Convergence test and the relative convergence problem,” in Computer Techniques for Electromagnetics. Oxford: Pergamon, 1973, pp. 300—304.[11] B.W. Kernighan and P.J. Plauger, The Elements of Programming Style. New York:McGraw-Hill, 1974.[12] M. Metcalf, Effective Fortran 77. Oxford: Clarendon Press, 1985.[13] B. Kruger, Efficient Fortran Programming. New York: Wiley, 1990.[14] T. Nicol, Ed. UBC Matrix: A Guide to Solving Matrix Problems. Vancouver, B.C.: Univ.of British Columbia Computing Centre, Mar. 1982.[15] M. Loukides, UNIX for Fortran Programmers. Sebastopol, CA: O’Reilley & Associates,1990.Appendix CCIRCULAR POLARIZATION SELECTIVE REFLECTORSC.1 IntroductionBecause the polarization scattering matrix which corresponds to acircular polarization selectiveresponse cannot be diagonalized when expressed with respect to alinearly polarized basis, atrihedral corner reflector cannot be modified to present such aresponse simply by placing finsor corrugations of appropriate dimensionsand orientation along one of its interior surfaces oracross its aperture. Alternative methods for obtaining such a response based onthe additionof a suitable transmission polarizer to either a linear polarizationselective or a twist polarizingtrihedral corner reflector are proposed.C.2 ConceptA circular polarization selective reflector will return anincident wave of the chosen sense backto the source but will either absorb an incident wave ofthe opposite sense or scatter it in adifferent direction. With respect to a circularly polarized basis,the normalized polarizationscattering matrices of left and right circular polarizationselective reflectors are given bySLL SLR1 0[SJ = =(C.1)SRL SRR0 0andSLL SLR 0 0= . (C.2)SRL SRR0 1The co-polar and cross-polar response of circular polarizationselective reflectors are plottedas a function of the polarization state of theincident wave in Figures C.1 and C.2. The189Appendix C. Circular PolarizationSelective Reflectors190parameters E and T refer to the ellipticityand tilt angles of the corresponding polarizationellipse.Although the amplitudes of theco-polar and cross-polar responses of a circular polarizationreflector are invariant under rotationof the target about the line-of-sight, the relativephase ofthe response varies linearly withthe angle of rotation.In Chapter 4, it is shown that a conventionaltrihedral corner reflector can be modified topresent a given polarizationresponse along its boresight simply by placing conductingfins orcorrugations of appropriate dimensions andorientation along one of its three interiorsurfacesor across its aperture if thecorresponding polarization scattering matrix [S]can be transformedthrough rotation about the line-of-sight into a diagonalmatrix [S’] of the form[S’]= SH0, (C.3)0S(7A polarization scattering matrixexpressed with respect to a circularly polarizedbasis [Sj canbe transformed to the equivalent polarizationscattering matrix expressed with respect toalinearly polarized basis [S] byapplying the unitary change of basis transformation[SJ= {U*J[S]{U*]_l, (C.4)where [U] is the transformation matrix forchange of basis from linear to circular polarizationgiven by11 1—j. (C.5)Applying (C.4) to the matrices of (C.1) and (C.2)gives the normalized linear polarizationscattering matrices of left and right circularpolarization selective reflectors,[Se]= SHH SHy =—1(C.6)SVH Svv21andSHH SHy= 1—j(C.7)SVH Svv —j1Appendix C. Circular Polarization SelectiveReflectors 191(b)Figure C.1: Normalized response of a leftcircular polarization selective reflector as a functionof the polarization state of the incident wave.(a) Co-polar response. (b) Cross-polar response.(a)b.(b)Figure C.2: Normalized response of a right circularpolarization selective reflector as a functionof the polarization state of the incident wave. (a)Co-polar response. (b) Cross-polar response.(a)Appendix C. Circular Polarization Selective Reflectors 192If a depolarizing reflector with the polarization scattering matrix of (C.6) or (C.7) is rotatedabout the line-of-sight through an angle a, the off diagonal elements of its new polarizationscattering matrix [S’] are given byS1 = 2SHVcos2 a— SHy + (Svv — SHH)sin a cos a, (C.8)where the off diagonal elementsSHyandSVHare identical in all coordinate frames for themonostatic case [1]. Since the diagonal elements of the original matrix are real and the offdiagonal elements are imaginary, S, is complex. In order to diagonalize [5’], the angle ofrotation a must be chosen such that both the real and imaginary components of S, are setto zero, i.e.,Re(Sy)= (Svv—SHH)sinacosa = 0, (C.9)Im(Sy) = 2SHV cos2 a— SHy = 0.(C.10)Since SHHSvv,(C.9) is satisfied only if a = n r/2 where n is an integer. However, thesevalues of a will not satisfy (C.10) unlessSHy = 0.Therefore, the polarization scatteringmatrices of (C.6) and (C.7) cannot be diagonalized and it is not possible to realize a trihedralcorner reflector with a circular polarization selective response simply by placing conducting finsor corrugations of appropriate dimensions and orientation along one of its interior surfaces oracross its aperture.C.3 Proposed ImplementationsMethod IIt may be possible to realize a trihedral corner reflector with a circular polarization selectiveresponse by incorporating both transmission and reflection polarizers in the modified reflector.One such scheme is shown in Figure C.3. Here a transmission circular polarizer of the typeshown in Figure C.4 is placed across the aperture of a linear polarization selective reflector whichhas been oriented to return horizontally polarized incident waves but reject vertically polarizedincident waves. Subject to mechanical constraints, any of the several types of transmissioncircular polarizers that have been described in the literature (e.g., [2]) would be satisfactory.Appendix C. Circular Polarization Selective Reflectors193(a) (b)Figure C.3: A proposed implementation of a circularpolarization selective reflector using atransmission circular polarizer and a linear polarization selective reflector. (a)Left circularpolarization selective reflector. (b) Right circular polarization selective reflector.In Figure C.3(a), the axis of the transmission circular polarizer isoriented at -45 degreesto the vertical. Since left circularly polarized incidentwaves are converted to horizontal polarization by the transmission polarizer, they are returned by the linear polarizationselectivereflector then converted back to left circular polarization as they passthrough the transmissionpolarizer a second time. However, right circularly polarized incidentwaves are converted tovertical polarization by the transmission polarizer and are either scatteredaway from the sourceor absorbed by the linear polarization selective reflector. Thus, a left circularpolarization selective reflector has been realized. The polarization scattering matrix of the reflector isgivenby[Sc]= [T][S][T},— 11+j 1+j —1 011+j 1—j—21—j 1+j 0 021—j 1+j (C.11)— j—1j— jlwhere [T] and [5] are the polarization matrices of the transmission circular polarizerorientedas shown and a vertical polarization selective reflector, respectively, according to the BSAAppendix C. Circular Polarization Selective Reflectors 194LINEAR POLARIZATION4a 0.671 >>‘qfZX0 FREE SCEWAVELENGThCIRCULAR POLARIZATIONFigure C.4: A parallel plate transmission circular polarizer. (from [2, sec. 23-5])convention. If the transmission circular polarizer is oriented with its axis at 45 degrees to thevertical instead, as shown in Figure C.4(b), a right circular polarization selective reflector isrealized. If the response of the polarizer is now given by [T’], the polarization scattering matrixof the reflector will be given by[Sn] = [T’][S][T’j,11+j —1+j —1 011+j —1+j= 2—1+j 1+j 0 021+j‘ (C.12)——1—j— 2—j1Method IIA second scheme for implementing a circular polarization selective trihedral corner reflector isshown in Figure C.5. Here a circular polarization selective surface (CPSS) [4, 5] is placed acrossthe aperture of a trihedral twist reflector. Over most aspects, incident waves of the chosen sensepass through the CPSS and are returned back to the source by the trihedral twist reflector whileincident waves of the opposite sense are specularly reflected away from the source as suggestedby Figure C.6(a).Appendix C. Circular Polarization Selective Reflectors 195Figure C.5: A proposed implementation of a circular polarization selective reflector using a circular polarization selective surface and a trihedral twist reflector. (a) Left circular polarizationselective reflector. (b) Right circular polarization selective reflector.If the angle of incidence is normal to the reflector aperture, however, both incidentpolarizations will be returned to the source and the response will no longer appear to be polarizationselective, as shown in Figure C.6(b). For a reflector with corners of equal length, this condition occurs along the boresight or direction of maximum response which is clearly undesirable.This problem also limits the performance of the gridded trihedral, a linear polarization selectivereflector which is realized by mounting a closely spaced parallel grid of thin wiresacross thereflector aperture as shown in Figure C.7.(.R.L(R.L R(a)(b)Figure C.6: Scattering by a right circular polarization selective trihedral corner reflector.(a) Oblique incidence. (b) Normal incidence.(a) (b)Appendix C. Circular Polarization Selective Reflectors19650Figure C.7: A linear polarization selective gridded trihedral and itsco-polar and cross-polarazimuthal response patterns. (from [3])Backscatter returns from either the circular polarizationselective surface or the closelyspaced grid of parallel wires can be avoided by mounting thereflector in such a way that itsboresight is pointed away from the radar over most aspects.However, this is undesirable if thecorners of the reflector are of equal length because thescattering cross section of the reflectorwill be substantially reduced compared to its maximumvalue and will become much moresensitive to small changes in the orientation of the reflector. Analternative solution is to forcethe boresight and the normal to the aperture to pointin different directions by incorporatingbilateral symmetry into the reflector geometry. Considera trihedral corner reflector with onecorner of length c and two of length a as shown inFigure C.8. If the ratio of the center andside corner lengths of the reflector (c/a) is increasedslightly from unity, the normal to thereflector aperture will rise in elevation relative to the the reflector boresightangle as describedin Chapter 3. The elevation beamwidth of the reflectorwill also increase slightly at the expenseof its azimuthal beamwidth. Decreasing the ratio (c/a)instead would have the opposite effectalthough the angle of separation will be less pronounced.In either case, the angle of separation,8 between the normal to the aperture and the reflectorboresight will increase and the desiredresult will be achieved as suggested by Figure C.8(b).Trihedral Responseto Cross ComponentGrid Response toParallel ComponentC050 25 25Aspect AngleAppendix C. Circular Polarization Selective Reflectors 197aC(a)a(b)Figure C.8: A bilaterally symmetric trihedral corner reflector with triangular panels.0is theangle between the boresight and the normal to the aperture. (a) Front view. (b) Side view.The polarization scattering matrix of a hybrid reflector which combines a left CPSS withpolarization matrix [Tj and a trihedral twist reflector with polarization matrix [S] is given by[Se] = [T] {S} {T]1—1j—1 01—1j=j1 0 1• j1— 1—1j— jl(C.13)By substituting a right CPSS with polarization transmission matrix [T’] for the original, aright circular polarization selective reflector can be realized. Its polarization scattering matrixis given by[Sj = [T’]{S][T’]— 1—1.—j—1 01—1—j——j1 0 1—j1—1—1—j——j1(C.14)Appendix C. Circular Polarization Selective Reflectors198A method for realizing a circular polarization selective surface was first described byTilstonet a!. [4]. An alternative implementation which is easier and less expensive to fabricate thanthe original was subsequently devised by Morin [5]. The improved surface consists of a slab oflow dielectric constant material of thickness A/4 which has been divided into a square patternof cells of dimension A/2 by A/2. A resonant element which consists of a single piece of wirewhich has been bent into three sections is inserted into each cell such that the two end sectionsare flush against the front and back faces of the dielectric slab. Left and right CPSS elementsare shown in Figure C.9 where the z-y plane is the plane of the aperture.(b)Figure C.9: Elements of a circular polarization selective surface (CPSS). (a) Left CPSS element.(b) Right CPSS element.C.4 DiscussionTwo methods for realizing circular polarization selective trihedral corner reflectors have beenproposed in this appendix. The performance of such reflectors and their usefulness in practice will depend on a number of factors which have not been considered here including (1)degradation of the polarization response of either of the polarizers for incidence off the reflector boresight, (2) possible degradation of the polarization response due to multiple reflectionsbetween the transmission and reflection polarizers, and (3) the ease with which the modifiedreflector can be assembled and the mechanical ruggedness of the finished product.x3/8(a)References[1] S.H. Bickel, “Some invariant properties of the polarizationscattering matrix,” Proc. IEEE,vol. 53,pp.1070—1072, Aug. 1965.[2] W.B. Offut and L.K. DeSize, “Methods ofPolarization Synthesis,” in Antenna EngineeringHandbook, 2nd ed., R. C. Johnson and H. Jasik, Eds.New York: McGraw-Hill, 1984,chap. 23.[3] .J.A. Scheer, “Radar ReflectivityCalibration Procedures,” in Radar Reflectivity Measureinent: Techniques and Applications, N.C. Currie, Ed. Norwood, MA: Artech House,1989,p.109.[4] W.V. Tilston, T. Tralman, and S.M. Khanna, “Apolarization selective surface for circularpolarization,” in IEEE/A P-S Symp. Dig. (Syracuse, NY), June 1988, pp. 762—765.[5] G.A. Morin, “A simple circularpolarization selective surface (CPSS),” in IEEE/A P-SSymp. Dig. (Dallas, TX), May 1990,pp.100—103.199Appendix DEXPERIMENTAL ARRANGEMENTD.1 IntroductionDuring the course of this study, the microwave antenna rangelocated on the roof of the Electrical Engineering building at the University of British Columbia was upgradedand used tomeasure the response patterns of prototype depolarizing trihedral corner reflectors as describedin Chapter 4. The principles of radar cross section measurementhave been widely discussedin the literature, e.g., [1]—[5]. In this appendix, the modifications andimprovements that weremade to the UBC antenna range in support of the prototype reflectormeasurement programare briefly described. In section D.2, the general layout of the antenna range is described. Insections D.3 and D.4, respectively, the design and implementation of the CW radarapparatusand digital pattern recorder that are discussed. In section D.5, theresults of tests performed toverify the suitability of the antenna range for use in the measurement program are presented.Recommendations for future modifications and improvements are offered.D.2 OverviewA block diagram of the UBC microwave antenna range as itwas configured for the prototypereflector measurement program is shown in Figure D.1. The CW radar apparatus is designedto measure the co-polar and cross-polar response of targets in the range from 8—12 GHz. Thedigital pattern recorder is designed to calibrate the CW radar apparatus and record the response patterns of antennas or targets under test. It replaces the Scientific Atlantaseries 1520mechanical chart recorder used previously.200—E-B00 A A U’It1I CD CDm<:.:*Lm-...D2,•:..•:O0CD•:•.:•:••.IL[I1I4$?..G)S:-ai.0:•._z—‘ —I1G)WiTitImI\/_Xr°0—IoG) j%.)).3Cl,T.0 IN r’) 3I.Appendix D. Experimental Arrangement202The outdoor portion of the antenna range is shown in Figure D.2. The ScientificAtlantamodel 5851 model tower travels on a carriage along 15 metres of track down the center of theroof of the east wing of the Electrical Engineering building. The model tower supportstheantenna or target under test at a height of 2.2 metres and can be configured either to rotatethe device about a vertical axis for conventional azimuthal pattern measurements or toroll thedevice about a horizontal axis for polarization response measurements. Detailed mechanicalspecifications are presented in the model tower operating manual [6].D.3 CW Radar ApparatusThe CW radar apparatus consists of a CW transmitter, a microwave receiver equippedwith anexternal crystal mixer, and two standard gain horns mounted side by side on a custom-builtfeedthrough mounting adapter. A block diagram of the transmitter and a photograph of thetransmitter shelf are shown in Figures D.3 and D.4, respectively. The signal source is a MarconiFigure D.2: Photograph of the radar cross section measurement range.Appendix D. Experimental Arrangement 2036052B microwave signal generator which can provide at least 25 mW (14 dBm) of power to amatched load over the range from 8—12 GHz. The transmitted signal is sampled using an HPX752C 20 dB directional coupler and monitored using a Marconi 6593A VSWR meter equippedwith a HP X485B detector mount. The exact frequency of the transmitted signal is determinedby tuning the HP X532B frequency wavemeter until a dip is observed on the VSWR meter.The HP X382A precision attenuator is used to adjust the transmitter output power duringcalibration. Detailed specifications for each component are presented in references [7]—[9].The receiver section of the radar apparatus consists of a Scientific Atlanta model 171OAPportable microwave receiver equipped with a model M8.2 external crystal mixer. Detailedspecifications are given in the receiver operating manual [10]. The receiver local oscillator (LO)is tunable from 0.985—2.5 GHz. During operation, the LO output is fed through a RF pad tothe LO arm of a frequency selective tee. This tee couples the LO signal through a coaxial cableto the external mixer where harmonic mixing takes place. The resulting 45 MHz IF signal isconducted back through the same coaxial cable and frequency selective tee to a 45 MHz IFpreamplifier and subsequent stages. The receiver provides a signal to the monitor meter onits front panel which is proportional to the received signal strength. This signal was tappedand passed through a signal conditioning unit containing an op-amp based current-to-voltageconverter, a low pass filter, and an adjustable gain block in order to provide the 0—10 VDCoutput signal required by the digital pattern recorder.A pair of Scientific Atlanta model 12-8.2 standard gain horn antennas are used as thetransmitting and receiving antennas. The antennas present a boresight gain of 22.10 ± 0.05dB at a wavelength of 3.2 cm with E-plane and H-plane half-power beamwidths of 12.5 and13.5 degrees, respectively [11]. They are mounted side by side on a custom-built feedthroughmounting adapter which is attached to the side of the building penthouse at a height of 2.2metres as shown in Figures D.5 and D.6. The transmitting horn is mounted so that it isvertically polarized. The mounting adapter permits the receiving horn to be rotated by 90degrees in order to permit either the co-polar (0vv) and cross-polar (CHV) response of theAppendix D. Experimental Arrangement204target to be measured. Arrangements for aligning thehorns are shown in Figure D.7. Coarsealignment is performed using bubble levels which are mounted on thehorns and waveguideas shown. The alignment is completedby optical sighting an alignment target mounted onthe model tower using cross hairs mounted on thefront of the horns and the flanges of thewaveguide.Figure D.3: Block diagram of the CW radar transmitter.Standard Gain HornScientific AtlantaModel 12-8.2FrequencyMeterHPXS32BPrecision 10dBVariable DirectionalAttenuator CouplerHP X382A HP X752CDetectorMountHPX485BFigure D.4: Photograph of the CW radar transmitter.CD 0••ti CD CD CD CD+-CDCl)—CD—. CDCD Cl) CIaI.01Appendix D. Experimental Arrangement206Figure D.6: Profile view of the RCS measurement range.— a2mOm(a)Standard Gain HornFront viewCross-hairs25cm(b)Figure D.7:Arrangement for mechanically aligning the transmitting and receiving antennas.(a) Alignmentaids mounted on the transmitter shelf. (b) Alignment target mounted on themodel tower.Horn Antennas-3dBOm 5m lOm 15mCircular Levels90° twist section(as required)Side\1iewFeedthrough Mounting AdapterAppendix D. Experimental Arrangement 207A graphical representation of the link budget for the radar cross section measurement rangeas configured for measurement of prototype trihedral corner reflectors with a corner length of60 cm is presented in Figure D.8. The results correspond to the ideal case in which the poweravailable to the radar receiver is given by the radar range equation,p— PtGtGr(/A2)— (47r)3(R/A)4whereG = G,.= 22.1 dB are the gains of the transmitting and receiving antennas, respectively,u = 27.3 dBsm is the radar cross section of the target, ..\ = 3.18 cm is the radar wavelength,and R = 11 m is the range to the target. In order to achieve the required dynamic range, theminimum return from the target must be greater than either the background reflectivity of thefacility or the receiver thermal noise level.Ea-19ci)-J-49ci)Facility Background Level (Crosspolar)Receiver Thermal Noise Level140Transmitted Power LevelMaximum RCS ReturnDynamic RangeMinimum RCS Return0dB-10dB-20dB-30dBFacility ckgroundJçvel(opplar)Figure D.8: Radar cross section measurement range link budget.Appendix D. ExperimentalArrangement208D.4 Digital Pattern RecorderThe digital pattern recorderconsists of an industry standard personal computerequipped witha colour VGA display, an 80386SX-25 CPUwith a numerical coprocessor, 2 MB RAM,an80 MB hard disk drive, andcustom-designed interface hardware anddata acquisition software.A PC-26 analog-to-digital converter (ADC)expansion card (Boston Technologies, Boston, MA)equipped with the custom-designed signalconditioning unit described in the previous section is used to measure theoutput from the portable microwave receiver. Acustom-built12-bit synchro-to-digital converter (SDC) expansioncard based on a Control Sciences Inc.(Chatsworth, CA) 168F309 integrated12-bit SDC module is used to measure theangle ofrotation of the antenna or target undertest. Design considerations for PC-basedlaboratoryinstrumentation have been widely discussedin the literature, e.g., [12], [13]. A photographofthe digital pattern recorder, the positionercontrol unit, and the portable microwavereceiver isshown in Figure D.9.Figure D.9: Photograph of the digital pattern recorder,positioner control unit, and portablemicrowave receiver.Appendix D. Experimental Arrangement209A set of five independent program modules perform calibration, data acquisition, and diagnostic functions. The program modules store calibration data, configurationnotes, and responsepattern data in ASCII format and employ a consistent user interface.They were coded usingTurbo Pascal, version 5.5 (Borland International, Scotts Valley, CA) andmake use of graphicsroutines from the Science and Engineering Tools for Turbo Pascal subroutinelibrary, version 6.1(Quinn-Curtis, Needham, MA). The function of each module is summarized in Table D.1.Table D.1: Digital Pattern Recorder Program ModulesModule Name FunctionSynchro Test DPRSYN • display position of model tower in real time andverify correct operation of the synchro interfaceReceiver Calibration DPR_CAL • generate relative calibration curve for receiver andverify correct operation of the receiver interfacePattern Recorder DPR..PAT • record pattern data and configuration notesfor analysis and presentationData View DPR_VIEW • view and compare previously recorded dataReceiver Stability DPRSTAB • sample receiver output signal and generateamplitude distributions and frequency spectraThe synchro test module, DPLSYN, is used to verify that the model tower synchro, thepositioner control unit synchro repeater, the synchro-to-digital converter and the interconnectingcables are functioning correctly. It can also be used to calibrate the speed control on the positioner control unit. The program displays the current angle of rotation and rate of rotationofthe antenna or target under test while it generates a plot of the angle of rotation versustime.A sample screen display is presented in Figure D.10.The receiver calibration module, DPR_CAL, is used to measure the transfer characteristic ofthe microwave receiver. During a measurement sequence, the pattern recorder module usesthisdata to translate the binary code read from the ADC into a relative measure of the receivedsignal strength in decibels. Relative calibration of the CW radar apparatus is performed withthe equipment configured as shown in Figure D.11(a). After the operatorhas set the receiverinput signal to the highest level likely to be encountered duringthe measurement, the programAppendix D. Experimental Arrangement 210prompts the operator to decrement the received signal strength in 5 dB steps using either theprecision microwave attenuator in the transmitter or the IF step attenuator in the receiver andsample the receiver output. A sample screen display is presented in Figure D.10. After thecalibration sequence is complete, the operator is prompted to either save the relative calibrationdata to a ifie or begin the calibration sequence again. Absolute calibration of the CW radarapparatus can be performed either directly by measuring the response of a calibration target ofsimilar size to the target under test or indirectly by measuring the effective path loss betweenthe transmitter and receiver.An arrangement for performing absolute calibration using an indirect method is shown inFigure D.12. First, the response of the target is measured using the configuration of Figure D.12(a). Next, the crystal mixer is removed from the receiving horn and mounted on the10 dB directional coupler in the transmitter as shown in Figure D.12(b). Finally, the precisionmicrowave attenuator is adjusted until an identical response is observed at the receiver output.The total insertion loss includes the contributions of both the attenuator and the directionalcoupler. Since the distance to the target, the radar wavelength, and the gain of transmittingand receiving antennas are known, the absolute radar cross section of the target can be determined simply by equating the geometric path loss predicted by the radar equation to the totalinsertion loss. Closure is obtained when the results of direct and indirect calibration agree towithin a prescribed tolerance.The pattern recorder module, DPRPAT, is used to record both the parameters of the testconfiguration and the response pattern of the antenna or target under test. A sample parameter entry screen is shown in Figure D.13. When a new measurement sequence is begun, theparameter entry screen is replaced by a data acquisition screen which displays the responsepattern as it is being measured. The data acquisition screen can be configured to present theresults on either a rectangular or polar chart spanning either + 90 or +180 degrees. A sample± 180 degree polar display with a dynamic range of 30 dB is shown in Figure D.14. After datacollection is complete, the operator is returned to the parameter entry screen and prompted toeither save the response pattern data to a file or begin the measurement sequence again.Appendix D. Experimental Arrangement211Figure D.10: Digital pattern recorder: synclirotest screen.Use arrow keys for nenu selectionFigure D.11: Digital pattern recorder: receiver calibrationscreen.Appendix D. Experimental Arrangement212(a)TransmittingChannel (V)ReceivingChannel (VocH)(b)Figure D.12: Equipment configuration for performing relative and absolute calibration of theCW radar apparatus.90 twist section(required)Digital Pattern RecorderDigital Pattern RecorderAppendix D. Experimental Arrangement213X Axis Label:X Axis Scale:Sa,ipling Interval:Use arrou ke!Js for eenu eiect ionFigure D.13: Digital pattern recorder: parameter entry screen.Output Folder;Frequency:Polarization:9.4 GUzJert 1Output Filenai,e: stiCalibration Status-Calibration File:Plot TeAveraging:Receiver iplitude: Standard GainiFix4.... .w3Test Description:V Axis Label:Top Of Scale:Dmanic Range: 3BEGIN EX! TIFigure D.14: Digital pattern recorder: data acquisition screen.Appendix D. Experimental Arrangement214The last two program modulesare test articles used to evaluate functions whichwill beincorporated into future versions of the receiver calibrationand pattern recorder modules. Thereceiver stability module, DPR..STAB, is used to assess the shortand long-term stability of theportable microwave receiver while a constant amplitudesignal is applied to its input. After theoutput of the receiver has been sampled over a period oftime ranging from several minutes toseveral hours, the results are either viewed directly astime series or processed to yield amplitudedistribution functions and/or frequency spectra.The data viewing module, DPR_VIEW, is amodified version of the pattern recorder modulewhich is used to view and compare previouslyrecorded data.D.5 Facility EvaluationRadar cross section measurements are affected bya combination of random and systematicerrors. It is convenient to depict these errors bythe signal flow graph presented in FigureD.15 where S represents the actual response of the target, Rand T represent errors due to(1) deviation of the incident field from a planewave and (2) multipath reflections from thesurrounding facility, and I represents the contribution of generalbackground due to (3) returnsfrom the target support structure and surroundings and (4) directtransmission between thetransmitting and receiving antennas. In order to determine the suitabilityof a RCS facility foruse in a measurement program, a series of tests must beperformed to assess the magnitude ofthese errors and their potential impact on the accuracy of the resultsobtained.TransmitterIReceiverRFigure D.15: Error model for radar cross section measurement.Appendix D. Experimental Arrangement215Ideally, the target is located sufficiently far away from the transmitting antenna that theincident field presents a planar wavefront. For a small outdoor antenna range of the typeconsidered here, such an arrangement is rarely practical since the contribution of the multipathray to the response increases rapidly with range. It is often suggested thatsatisfactory resultswill be obtained if the phase deviation over the largest dimension of the target is less than 7r/8radians. This leads to the so-called far-field criterion which requires that the range to the targetbe greater than2 d2where d is the largest dimension of the target and A is the radar wavelength. In practice,physicalconstraints often make it necessary to measure radar cross section of targets at ranges aslow asone-quarter of the recommended value. The consequences are similar to those encountered whenconducting antenna measurements under similar circumstances. In particular, the apparentvalue of the radar cross section obtained by direct solution of theradar range equation will beless than the actual value that would have been obtained had the incident field been aplanewave.A trihedral corner reflector of the type used in the measurementprogram has triangularpanels of equal length and presents a maximum radar cross section of4ir £umax =where £ is the corner length of the reflector and A is the radar wavelength.Although the largestphysical dimension of a reflector of this type is the distance from the tip ofone corner to themid-point of the opposite panel, the maximum dimension of the effectiveaperture is muchsmaller and is given byas suggested by Figure D.16. The maximum aperture dimension, maximum radar crosssection,and far-field range of trihedral corner reflectors with various corner lengths aresummarized inTable D.2.Appendix D. Experimental Arrangement 216(a)Figure D.16: The effective aperture of a trihedral corner reflector with triangular panelsforincidence (a) along the boresight and (b) at an azimuth angleof 30 degrees.Table D.2: Response of Trihedral Corner Reflectors withTriangular Panels at9.445 GHzCorner length Max. aperture MaximumRCS Far-field range(cm) dimension (cm) (m2) (dBsm) (m)15 14 2 3.2 230 28 33 15.3 545 42 169 22.3 1160 57 537 27.3 2075 71 1310 31.2 31Four tests were conducted in order to assess the suitabilityof the antenna range for use inthe prototype reflector measurement program. All the tests were performedat the standardmarine radar frequency of 9.445 GHz.In the first test, the polarization response of the CWradar apparatus was evaluated. Thereceiving horn was removed from the feedthrough mounting adapter andattached to the modeltower. The response of the receiving horn was measured ata range of 11 metres as the hornwas rotated about its boresight. The results generallyagree with the predicted values andare presented in Figure D.17. The isolation between the transmittingand receiving horns isconsiderably greater than 30 dB when the horns areorthogonally polarized.In the next two tests, the contributions of near-field and multipatheffects to the responseof a target under test were evaluated.A conventional trihedrai corner reflector with a corner(b)Appendix D. Experimental Arrangement217length of 60 cm was attached to the model tower. First, the response of the target was measuredas its range was increased from 8 to 16 metres in 0.25 m increments.The results are presentedin Figure D.18. At ranges greater than about 12 metres, rapid variationsin the response areobserved which suggest that multipath effects are becoming important.At ranges less thanabout 9 metres, evidence of a reduction in the response whichis apparently due to near-fieldeffects is observed. At ranges between9 and 12 metres, the response varies as hR4 with onlya few perturbations. A second test was performed in order to assess the effectof small changesin range on the fine structure of the response pattern. The azimuthal response pattern of thereflector was measured at ranges of 10 and 12 metres and compared. The results arepresentedin Figure D.19 where it can be seen that the response patterns are nearly identical. On this basisof these results, it was concluded that the optimum range at which to measure the prototypereflectors is 11 metres.In the last test, the contributions of direct transmission between the horns and unwantedreturns from the model tower and surroundings were evaluated. The conventional trihedralcorner reflector was removed from the model tower and the positioner head was coveredby asmall section of microwave absorber. The co-polar and cross-polar responses of the backgroundwas measured in turn as the model tower was rotated in azimuth ata range of 11 metres.The results are presented in Figure D.20. The co-polar response rises above-30 dB at oniy afew angles and never rises above -28 dB relative to the maximum responseof the prototypereflector. The cross-polar response never rises above -30 dB.During the evaluation of the facility, several problems were noted. The gears that rotatethepositioner head about a horizontal axis exhibit a small but noticeable backlash. The effectisparticularly obvious when a large target such as a prototype reflector is mounted on themodeltower. Since procedures for reducing the backlash require the services of a skilled machinist[6],the problem could not be corrected simply by conducting routine maintenance during the courseof the measurement program. It is recommended that the problem be corrected in the nearfuture, however.Appendix D. Experimental Arrangement 2180.0-2.0I&-10.0-12.08.0 10.0 12.0 14.0 16.0Range (m)Figure D.18: Boresight response of a conventional trihedral corner reflector vs. range.0•0dB-10dB-20 dB-30 dB-Measured•...Theoretical.09kFigure D.17: Polarization response ofthe receiving horn at a range of 11 m.— MeasuredIdealAppendix D. Experimental Arrangement2190.0-10.01-20.0— MeasuredTheoretical-30.0-45.0 -30.0 -15.0 0.0 15.0 30.045.0Azimuth Angle-(deg)Figure D.19: Azimuthal response pattern of a conventional trihedral cornerreflector at rangesof 10 and 12 m.0.0_vvU:2-bC;)10.0a,4-a,a,-20.0-30.0 .. it\. . . . fl ,i. .A.-45.0 -30.0 -15.0 0.0 15.0 30.0 45.0Azimuth Angle-,(deg)Azimuthal response pattern of the model tower at a range of 11 m.Figure D.20:Appendix D. Experimental Arrangement 220The transfer characteristic of the portable microwave receiver tends to drift noticeably afteronly an hour. As a result, it was necessary to recalibrate the receiver each time a new seriesof measurements was conducted. Since the receiver is over twenty years old and is nearing theend of its useful life, it is recommended that consideration be given to replacing it.As noted earlier, multipath reflections were observed when the range to the target wasgreater than about 12 metres. It might be possible to significantly reduce multipath effects(and extend the model tower’s useful range of travel) by employing either a berm or a seriesof radar fences to scatter the multipath ray as described by Knott [2, pp. 369—370]. It isrecommended that an experimental program be conducted to assess the effectiveness of suchmethods in suppressing multipath effects at this antenna rangeReferences[1] C. G. Bachman, Radar Targets. Lexington, MA: D. C. Heath, 1982,pp.109—158.[2] E.F. Knott, J.F. Shaeffer, and M.T. Tuley, Radar Cross Section: ItsPrediction, Measurement, and Reduction. Norwood, MA: Artech House, 1985,pp.315—382.[3] Radar Cross-Section Measurements with the HP 8510Network Analyzer. Product NoteNo. 8510-2, Santa Rosa, CA: Hewlett-Packard, Apr. 1985.[4] R.B. Dybdal, “Radar cross section measurements,” Proc. IEEE,vol. 75,pp.498—516,Apr. 1987.[5] N.C. Currie, Ed., Radar Reflectivity Measurement:Techniques and Applications, Norwood,MA: Artech House, 1989.[6] Series 5800 Model Towers, 3rd ed. Atlanta, GA:Scientific-Atlanta, Dec. 1970.[7] Instruction Manual for Signal SourceType 6058B. Stevenage, UK: Marconi Instruments(Sanders Div.), 1976.[8] Coaxial Waveguide Catalog andMicrowave Measurement Handbook. Palo Alto, CA:Hewlett-Packard, 1979.[9] Instruction Manual for VSWR Indicator and Selective AmplifierType 6593A. Stevenage,UK: Marconi Instruments (Sanders Div.), 1976.[10] Series 1710 Portable Microwave Receiver, 2nd ed.Atlanta, GA: Scientific-Atlanta,May 1970.[11] W. T. Slayton, Design and Calibration of Microwave Antenna GainStandards. NRL Report. No. 4433, U.S. Naval Research Lab, Washington, DC, 1954.221References222[12] S.C. Gates and J. Becker, Laboratory Automation usingthe IBM PC. Englewood Cliffs,NJ: Prentice-Hall, 1989.[13] B.G. Thompson and A.F. Kuckes, IBM-PC in the Laboratory.Cambridge: CambridgeUniv. Press, 1989.
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Radar cross section enhancement for radar navigation...
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Radar cross section enhancement for radar navigation and remote sensing Michelson, David George 1993
pdf
Page Metadata
Item Metadata
Title | Radar cross section enhancement for radar navigation and remote sensing |
Creator |
Michelson, David George |
Date | 1993 |
Date Issued | 2009-04-14T18:43:59Z |
Description | Recent developments in radar navigation and remote sensing have led to a requirement for rugged yet inexpensive location markers and calibration targets which present both a very large scattering cross section and a specified polarization response over a wide angular range. This study considers several problems related to the analysis and design of passive radar targets derived from corner reflectors. Transformation of the polarization response of a target between global and local coordinate frames is shown to correspond to rotation of the polarization basis by a prescribed angle which is a function of both the coordinate transformation matrix and the direction of propagation. Once the angle of rotation has been determined using either spherical trigonometry or vector algebra, any polarization descriptor can be transformed between coordinate frames by application of a suitable rotation operator. The scattering cross section and angular coverage of a conventional trihedral corner reflector can be altered by modifying the size and shape of its reflecting panels. A numerical algorithm based on physical optics is used to predict the contribution of triple-bounce reflections to the response of a reflector with polygonal panels of arbitrary shape. If three-fold symmetry is broken and the reflector is simply required to present bilateral symmetry, it is found that the scattering cross section, elevation beamwidth, and azimuthal beamwidth of the reflector can be chosen independently of each other. A method for altering the polarization response of a conventional trihedral corner reflector by adding conducting fins or corrugations to one its interior surfaces is proposed. In calculating design curves for twist-polarizing or circularly polarizing reflectors by mode-matching, optimum accuracy and efficiency are obtained by setting the ratio of free space to groove modes equal to the ratio of groove width to the period. Methods for obtaining linear and circular polarization selective responses are considered. The contribution of triple-bounce reflections to the response of such reflectors is a function of the direction of incidence, the orientation of the reflector, the dimensions of the corrugations, and the size and shape of the reflecting panels. Experimental results show that prototype twist-polarizing and circularly polarizing reflectors respond essentially as predicted. |
Extent | 4473892 bytes |
Genre |
Thesis/Dissertation |
Type |
Text |
File Format | application/pdf |
Language | eng |
Collection |
Retrospective Theses and Dissertations, 1919-2007 |
Series | UBC Retrospective Theses Digitization Project |
Date Available | 2009-04-14 |
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.0065024 |
URI | http://hdl.handle.net/2429/7038 |
Degree |
Doctor of Philosophy - PhD |
Program |
Electrical and Computer Engineering |
Affiliation |
Applied Science, Faculty of Electrical and Computer Engineering, Department of |
Degree Grantor | University of British Columbia |
Graduation Date | 1994-05 |
Campus |
UBCV |
Scholarly Level | Graduate |
Aggregated Source Repository | DSpace |
Download
- Media
- [if-you-see-this-DO-NOT-CLICK]
- ubc_1994-894504.pdf [ 4.27MB ]
- Metadata
- JSON: 1.0065024.json
- JSON-LD: 1.0065024+ld.json
- RDF/XML (Pretty): 1.0065024.xml
- RDF/JSON: 1.0065024+rdf.json
- Turtle: 1.0065024+rdf-turtle.txt
- N-Triples: 1.0065024+rdf-ntriples.txt
- Original Record: 1.0065024 +original-record.json
- Full Text
- 1.0065024.txt
- Citation
- 1.0065024.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Country | Views | Downloads |
---|---|---|
United States | 23 | 1 |
China | 20 | 5 |
Singapore | 19 | 0 |
India | 16 | 0 |
Russia | 11 | 0 |
France | 10 | 2 |
Canada | 9 | 0 |
Germany | 6 | 125 |
United Kingdom | 5 | 4 |
Japan | 4 | 0 |
Unknown | 3 | 0 |
New Zealand | 2 | 0 |
Brazil | 2 | 1 |
City | Views | Downloads |
---|---|---|
Unknown | 58 | 132 |
Beijing | 11 | 3 |
Hangzhou | 6 | 0 |
Mumbai | 5 | 0 |
Ashburn | 5 | 0 |
Penza | 4 | 0 |
Montreal | 4 | 0 |
Tokyo | 4 | 0 |
Indore | 3 | 0 |
Wilmington | 3 | 0 |
Cheyenne | 2 | 0 |
Uberlandia | 2 | 0 |
Saint Petersburg | 2 | 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-0065024/manifest