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

Problema

O mapeamento por Radar de Abertura Sintética (SAR) de alvos móveis sobre terrenos tridimensionais (ou no ar) apresenta um desafio duplo de estimação: determinar com precisão tanto o deslocamento no cross-range induzido pelo movimento quanto a altura do alvo. Embora os sistemas convencionais de Velocity SAR (VSAR) que utilizam arranjos de antenas lineares possam estimar deslocamentos de cross-range via Transformada de Fourier Discreta (DFT), eles falham em resolver a altura do alvo porque amostram espacialmente apenas ao longo da direção azimutal. Além disso, quando as velocidades ou alturas dos alvos são grandes, o envoltório de fase (phase wrapping) inerente ao processamento por DFT leva à ambiguidade, restringindo o intervalo inequívoco de estimação. Os métodos existentes para resolver a ambiguidade, como arranjos lineares de frequências múltiplas ou não uniformes, são projetados primariamente para cenários unidimensionais (1D) e não abordam diretamente a natureza acoplada em 2D do deslocamento de cross-range e da estimação de altura em ambientes 3D.

Metodologia

O artigo propõe um framework geral utilizando arranjos de antenas planares e o Teorema Chinês dos Restos Multidimensional (MD-CRT) para estimar conjuntamente o deslocamento de cross-range e a altura do alvo.

1. Modelo de Sinal e Formulação 2D-DFT:

   • Os autores modelam a geometria SAR com uma plataforma movendo-se a uma velocidade constante e um arranjo de antena planar (situado no plano $x-z$) com elementos dispostos em uma rede gerada por uma matriz $M$.
   • Assume-se que o alvo se move com um vetor de velocidade 3D.
   • Após o processamento padrão de SAR (incluindo compensação de fase quadrática), a imagem complexa através do arranjo planar exibe uma estrutura exponencial complexa 2D.
   • O deslocamento de cross-range induzido pelo movimento e a altura do alvo são unificados em um único vetor de parâmetros 2D $\mathbf{g}$.
   • A aplicação de uma Transformada de Fourier Discreta 2D (2D-DFT) através dos índices da antena revela que o vetor de parâmetros $\mathbf{g}$ está relacionado ao pico detectado via uma formulação de resto vetorial módulo uma matriz de inteiro específica $N$ (o módulo da matriz 2D-DFT). Especificamente, $\mathbf{g} \equiv \mathbf{r} \pmod{N}$, onde $\mathbf{r}$ é o vetor resto observado.

2. Construção do Módulo de Matriz Ótimo:

   • Ao contrário dos arranjos lineares 1D onde o módulo é um escalar, o cenário 2D envolve um módulo de matriz $N$. A estrutura de $N$ afeta como os erros no vetor resto se propagam para o vetor de parâmetros estimado.
   • O artigo deriva uma construção ótima para $N$ sob uma geometria de arranjo fixa ($M$) e um número fixo de elementos de antena (fixo $|\det(N)|$).
   • A escolha ótima é $N = \alpha \text{adj}(M^\top)$, que minimiza o fator de amplificação de erro $\|M^{-\top}N^{-1}\|_\infty$, garantindo que os erros de resto sejam equilibrados entre ambas as dimensões em vez de serem amplificados em uma única direção.

3. Framework de Multi-subarranjos para Resolução de Ambiguidade:

   • Um único arranjo planar produz apenas uma estimativa "dobrada" (o vetor resto) se o verdadeiro vetor de parâmetros estiver fora do intervalo inequívoco definido pelo Paralelepípedo Fundamental (FPD) de $N$.
   • Para resolver isso, os autores propõem um framework de multi-subarranjos utilizando $J$ subarranjos, cada um com um gerador de rede diferente $M_j$ e o correspondente módulo $N_j$.
   • Isso gera um sistema de congruências de matriz: $\mathbf{g} \equiv \beta_j \mathbf{r}_j \pmod{\beta_j N_j}$.
   • O verdadeiro vetor de parâmetros $\mathbf{g}$ é recuperado resolvendo este sistema usando o MD-CRT, que efetivamente amplia o intervalo inequívoco para o FPD do Mínimo Múltiplo Comum à Direita (LCRM) dos módulos.

4. Recuperação Robusta com Erros Práticos:

   • O artigo aborda questões práticas onde o vetor de frequência não é inteiro, levando a erros de quantização na detecção de pico da 2D-DFT.
   • Um modelo de pico 2D aproximado é introduzido, substituindo vetores de frequência não inteiros por vetores inteiros próximos para mapear o problema de volta ao framework padrão de vetor resto de inteiros.
   • O MD-CRT Robusto é incorporado para lidar tanto com erros de quantização (caso sem ruído) quanto com ruído de sinal aditivo.
   • Os autores estabelecem condições suficientes para recuperação robusta e derivam limites explícitos de erro de reconstrução para os cenários sem ruído e com ruído. Esses limites dependem dos fatores de escala $\alpha_j$ e dos níveis de ruído.

Principais Contribuições

• Formulação 2D Unificada: O artigo estabelece uma formulação de vetor resto 2D unificada para a estimação conjunta do deslocamento de cross-range e altura do alvo usando arranjos planares, estendendo os conceitos anteriores de VSAR 1D.
• Otimização do Módulo de Matriz: Fornece uma derivação teórica para o módulo de 2D-DFT ótimo que minimiza a propagação de erro, uma liberdade de design estrutural única ao cenário 2D.
• Framework MD-CRT para SAR: Adapta o MD-CRT para o mapeamento SAR, demonstrando como múltiplos subarranjos planares com diferentes geradores de rede podem resolver conjuntamente ambiguidades e ampliar o intervalo inequívoco além do que é possível com um único arranjo.
• Análise de Robustez: O trabalho introduz um modelo de pico aproximado para frequências não inteiras e fornece condições suficientes rigorosas e limites de erro para recuperação robusta na presença de quantização e ruído aditivo.

Resultados

O artigo apresenta simulações numéricas para validar o framework teórico:

• Ampliação do Intervalo Inequívoco: Uma simulação com dois subarranjos planares coprimos demonstra que o intervalo inequívoco conjunto (definido pelo LCRM dos módulos dos subarranjos) é estritamente maior do que o de qualquer subarranjo individual.
• Precisão de Localização: Em um cenário sem ruído, o framework de dois subarranjos recupera com sucesso a localização do alvo na estrada. Em contraste, o uso de um único subarranjo resulta em uma localização de alvo reconstruída deslocada do caminho real devido à ambiguidade.
• Limites de Erro: Os limites de erro derivados confirmam que o aumento do fator de escala $\alpha$ (efetivamente aumentando o tamanho do arranjo ou usando zero-padding) reduz o erro de reconstrução, melhorando a precisão.

Significância e Alegações

Os autores afirmam que este trabalho fornece a fundação teórica para o design de arranjos planares e análise de desempenho no mapeamento SAR de alvos móveis.

• O framework permite a estimação conjunta de parâmetros de movimento 3D (deslocamento e altura), o que não é alcançável com arranjos lineares 1D convencionais.
• O método de multi-subarranjos proposto oferece um método sistemático para ampliar o intervalo inequívoco sem exigir mudanças complexas de hardware além da geometria do arranjo.
• O artigo afirma explicitamente que a Parte II desta série de duas partes investigará designs específicos de subarranjos planares e fornecerá simulações mais detalhadas sobre robustez, indicando que a Parte I foca no framework geral, formulação matemática e limites teóricos.

A significância reside em preencher a lacuna entre as técnicas de VSAR 1D e os requisitos mais complexos de mapeamento de terreno 3D, oferecendo uma abordagem matematicamente rigorosa para a resolução de ambiguidades e gestão de erros em sistemas SAR.