Moving Target SAR Imaging Using Planar Arrays And Multidimensional Chinese Remainder Theorem (MD-CRT)--Part I: A General Framework

Probleemstelling

Synthetische Apertuur Radar (SAR) beeldvorming van bewegende doelen boven driedimensionaal terrein (of in de lucht) presenteert een dubbele schattingsuitdaging: het nauwkeurig bepalen van zowel de door beweging veroorzaakte cross-range verschuiving als de hoogte van het doel. Terwijl conventionele Velocity SAR (VSAR)-systemen die gebruikmaken van lineaire antenne-arrays cross-range verschuivingen kunnen schatten via Discrete Fourier Transform (DFT), falen zij in het oplossen van de doelhoogte omdat zij slechts ruimtelijk langs de azimuth-richting bemonsteren. Bovendien leidt bij grote doelvelociteiten of -hoogten de fase-wrapping die inherent is aan DFT-verwerking tot ambiguïteit, wat het onambigue schattingsbereik beperkt. Bestaande methoden voor het oplossen van ambiguïteit, zoals multi-frequentie of niet-uniforme lineaire arrays, zijn primair ontworpen voor scenario's in één dimensie (1D) en adresseren niet direct de gekoppelde 2D-aard van cross-range verschuiving en hoogte-schatting in 3D-omgevingen.

Methodologie

Het artikel stelt een algemeen framework voor dat gebruikmaakt van planaire antenne-arrays en de Multidimensionale Chinese Reststelling (MD-CRT) om de cross-range verschuiving en de hoogte van het doel gezamenlijk te schatten.

1. Signaalmodel en 2D-DFT Formulering:

   • De auteurs modelleren de SAR-geometrie met een platform dat beweegt met een constante snelheid en een planaire antenne-array (liggend in het $x-z$ vlak) met elementen gerangschikt op een rooster gegenereerd door een matrix $M$.
   • Er wordt vanuit gegaan dat het doel beweegt met een 3D-snelheidsvector.
   • Na standaard SAR-verwerking (inclusclusief kwadratische fasecompensatie) vertoont het complexe beeld over de planaire array een 2D complexe exponentiële structuur.
   • De door beweging veroorzaakte cross-range verschuiving en de hoogte van het doel worden verenigd in een enkele 2D-parametervector $\mathbf{g}$.
   • Het toepassen van een 2D Discrete Fourier Transform (2D-DFT) over de antenne-indices laat zien dat de parametervector $\mathbf{g}$ gerelateerd is aan de gedetecteerde piek via een vectorrestformulering modulo een specifieke integer matrix $N$ (de 2D-DFT matrix modulus). Specifiek geldt $\mathbf{g} \equiv \mathbf{r} \pmod{N}$, waarbij $\mathbf{r}$ de geobserveerde vectorrest is.

2. Optimale Constructie van de Matrix Modulus:

   • In tegen af van 1D lineaire arrays waar de modulus een scalar is, betreft de 2D-setting een matrix modulus $N$. De structuur van $N$ beïnvloedt hoe fouten in de vectorrest voortplanten naar de geschatte parametervector.
   • Het artikel leidt een optimale constructie voor $N$ af onder een vaste array-geometrie ($M$) en een vast aantal antenne-elementen (vaste $|\det(N)|$).
   • De optimale keuze is $N = \alpha \text{adj}(M^\top)$, wat de foutversterkingsfactor $\|M^{-\top}N^{-1}\|_\infty$ minimaliseert, waardoor wordt gewaarborgd dat fouten in de restwaarde evenredig verdeeld zijn over beide dimensies in plaats van versterkt te worden in één richting.

3. Multi-Subarray Framework voor Ambiguïteitsresolutie:

   • Een enkele planaire array levert slechts een "gevouwen" schatting (de vectorrest) op als de ware parametervector buiten het onambigue bereik valt dat wordt gedefinieerd door het Fundamentele Parallelepiped (FPD) van $N$.
   • Om dit op te lossen, stellen de auteurs een multi-subarray framework voor met $J$ subarrays, elk met een verschillende roostergenerator $M_j$ en een bijbehorende modulus $N_j$.
   • Dit genereert een systeem van matrix-congruenties: $\mathbf{g} \equiv \beta_j \mathbf{r}_j \pmod{\beta_j N_j}$.
   • De ware parametervector $\mathbf{g}$ wordt hersteld door dit systeem op te lossen met behulp van de MD-CRT, wat effectief het onambigue bereik vergroot naar het FPD van de Kleinste Gemeenschappelijke Rechter Veelvoudige (LCRM) van de moduli.

4. Robuust Herstel met Praktische Fouten:

   • Het artikel behandelt praktische kwesties waarbij de frequentievector geen geheel getal is, wat leidt tot kwantiseringsfouten in de 2D-DFT piekdetectie.
   • Een benaderd 2D-DFT piekmodel wordt geïntroduceerd, waarbij niet-gehele frequentievectoren worden vervangen door nabijgelegen gehele vectoren om het probleem terug te mappen naar het standaard integer vector-restframework.
   • Robuuste MD-CRT wordt geïntegreerd om zowel kwantiseringsfouten (ruisloos geval) als additieve signaalruis te behandelen.
   • De auteurs stellen voldoende voorwaarden vast voor robuust herstel en leiden expliciete reconstructiefoutgrenzen af voor zowel de ruisloze als de ruizige setting. Deze grenzen zijn afhankelijk van de schalingfactoren $\alpha_j$ en de ruisniveaus.

Belangrijkste Bijdragen

• Verenigde 2D-Formulering: Het artikel vestigt een verenigde 2D-vectorrestformulering voor de gezamenlijke schatting van cross-range verschuiving en doelhoogte met behulp van planaire arrays, waarmee eerdere 1D VSAR-concepten worden uitgebreid.
• Optimalisatie van de Matrix Modulus: Het biedt een theoretische afleiding voor de optimale 2D-DFT matrix modulus die foutvoortplanting minimaliseert, een structurele ontwerpvrijheid die uniek is voor de 2D-setting.
• MD-CRT Framework voor SAR: Het past de MD-CRT aan voor SAR-beeldvorming en demonstreert hoe meerdere planaire subarrays met verschillende roostergeneratoren gezamenlijk ambiguïteiten kunnen oplossen en het onambigue bereik kunnen vergroten buiten wat mogelijk is met een enkele array.
• Robuustheidsanalyse: Het werk introduceert een benaderd piekmodel voor niet-gehele frequenties en biedt rigoureuze voldoende voorwaarden en foutgrenzen voor robuust herstel in aanwezigheid van kwantisatie en additieve ruis.

Resultaten

Het artikel presenteert numerieke simulaties om het theoretische framework te valideren:

• Vergroting van het Onambigue Bereik: Een simulatie met twee co-prime planaire subarrays demonstreert dat het gezamenlijke onambigue bereik (gedefinieerd door de LCRM van de subarray-moduli) strikt groter is dan dat van een enkele subarray.
• Lokalisatienauwkeurigheid: In een ruisloos scenario herstelt het twee-subarray framework de locatie van het doel op de weg succesvol. In contrast hiermee resulteert het gebruik van een enkele subarray in een gereconstrueerde doelpositie die verschoven is ten opzichte van het werkelijke pad vanwege ambiguïteit.
• Foutgrenzen: De afgeleide foutgrenzen bevestigen dat het vergroten van de schalingfactor $\alpha$ (effectief het vergroten van de array-omvang of het gebruiken van zero-padding) de reconstructiefout vermindert, wat de nauwkeurigheid verbetert.

Betekenis en Claims

De auteurs stellen dat dit werk de theoretische basis vormt voor het ontwerp en de prestatieanalyse van planaire arrays in moving target SAR-beeldvorming.

• Het framework maakt de gezamenlijke schatting van 3D-bewegingsparameters (verschuiving en hoogte) mogelijk, wat niet haalbaar is met conventionele 1D lineaire arrays.
• Het voorgestelde multi-subarray aanpak biedt een systematische methode om het onambigue bereik te vergroten zonder dat complexe hardwarewijzigingen nodig zijn, afgezien van de array-geometrie.
• Het artikel vermeldt expliciet dat Deel II van deze tweeledige serie onderzoek zal doen naar specifieke planaire subarray-ontwerpen en meer gedetailleerde simulaties over robuustheid zal bieden, wat aangeeft dat Deel I zich richt op het algemene framework, de wiskundige formulering en de theoretische grenzen.

De betekenis ligt in het overbruggen van de kloof tussen 1D VSAR-technieken en de complexere vereisten van 3D-terreinbeeldvorming, door een mathematisch rigoureuze aanpak te bieden voor ambiguïteitsresolutie en foutbeheer in SAR-systemen.