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

다음은 "Moving Target SAR Imaging Using Planar Arrays and Multidimensional Chinese Remainder Theorem (MD-CRT)–Part I: A General Framework" 논문의 상세 기술 요약입니다.

문제 정의

3차원 지형(또는 공중) 상의 이동 표적에 대한 합성 개구 레이더(SAR) 이미징은 두 가지 추정 과제를 동시에 안고 있습니다: 운동에 의한 크로스 레인지 시프트(cross-range shift)와 표적 높이를 모두 정확하게 결정하는 것입니다. 선형 안테나 배열을 사용하는 기존의 벨로시티 SAR(VSAR) 시스템은 이산 푸리에 변变换(DFT)를 통해 크로스 레인지 시프트를 추정할 수 있지만, 공간 샘플링이 방위각 방향으로만 이루어지기 때문에 표적 높이를 분해할 수 없습니다. 또한, 표적의 속도나 높이가 클 경우 DFT 처리 과정에서 발생하는 위상 모호성(phase wrapping)으로 인해 추정 가능한 무모호 범위(unambiguous range)가 제한됩니다. 다중 주파수 또는 비균일 선형 배열과 같은 기존의 모호성 해결 방법들은 주로 1차원(1D) 시나리오를 위해 설계되었으며, 3차원 환경에서의 크로스 레인지 시프트와 높이 추정 사이의 결합된 2차원적 특성을 직접적으로 다루지 못합니다.

방법론

본 논문은 평면 안테나 배열과 다차원 중국인 나머지 정리(MD-CRT)를 활용하여 크로스 레인지 시프트와 표적 높이를 공동으로 추정하는 일반적인 프레임워크를 제안합니다.

1. 신호 모델 및 2D-DFT 정식화:

   • 저자들은 일정한 속도로 이동하는 플랫폼과 행렬 $M$에 의해 생성된 격자 위에 배치된 평면 안테나 배열($x-z$ 평면에 위치)을 갖는 SAR 기하 구조를 모델링합니다.
   • 표적은 3차원 속도 벡터를 가지고 이동한다고 가정합니다.
   • 표준 SAR 처리(이차 위상 보상 포함)를 거친 후, 평면 배열 전체의 복소 이미지는 2차원 복소 지수 구조를 나타냅니다.
   • 운동 유도 크로스 레인지 시프트와 표적 높이는 하나의 2차원 파라미터 벡터 $\mathbf{g}$로 통합됩니다.
   • 안테나 인덱스에 대해 2D-DFT를 적용하면, 파라미터 벡터 $\mathbf{g}$는 특정 정수 행렬 $N$(2D-DFT 행렬 모듈러스)에 대한 벡터 나머지 공식 $\mathbf{g} \equiv \mathbf{r} \pmod{N}$을 통해 관측된 벡터 나머지 $\mathbf{r}$과 연관됨을 알 수 있습니다.

2. 최적 행렬 모듈러스 구축:

   • 1차원 선형 배열에서 모듈러스가 스칼라인 것과 달리, 2차원 설정에서는 행렬 모듈러스 $N$이 관여합니다. $N$의 구조는 벡터 나머지의 오차가 추정된 파라미터 벡터로 어떻게 전파되는지에 영향을 미칩니다.
   • 본 논문은 고정된 배열 기하 구조($M$)와 고정된 안테나 소자 수(고정된 $|\det(N)|$) 하에서 $N$에 대한 최적 구축 방식을 도출합니다.
   • 최적의 선택은 $N = \alpha \text{adj}(M^\top)$이며, 이는 오차 증폭 계수 $\|M^{-\top}N^{-1}\|_\infty$를 최소화하여 나머지 오차가 한 방향으로 증폭되지 않고 두 차원에 걸쳐 균형 있게 배분되도록 보장합니다.

3. 모호성 해결을 위한 다중 서브어레이(Multi-Subarray) 프레임워크:

   • 단일 평면 배열은 실제 파라미터 벡터가 $N$의 기본 평행육면체(Fundamental Parallelepiped, FPD) 외부의 무모호 범위에 있을 경우, "접힌(folded)" 추정치(벡터 나머지)만을 제공합니다.
   • 이를 해결하기 위해, 서로 다른 격자 생성기 $M_j$와 그에 대응하는 모듈러스 $N_j$를 가진 $J$개의 서브어레이를 사용하는 다중 서브어레이 프레임워크를 제안합니다.
   • 이는 $\mathbf{g} \equiv \beta_j \mathbf{r}_j \pmod{\beta_j N_j}$ 형태의 행렬 합동식 체계를 생성합니다.
   • 실제 파라미터 벡터 $\mathbf{g}$는 MD-CRT를 사용하여 이 체계를 해결함으로써 복구되며, 이는 모듈러스들의 최소 공배 우측 다중(Least Common Right Multiple, LCRM)의 FPD로 무모호 범위를 효과적으로 확장합니다.

4. 실제적 오차를 고려한 강건한 복구:

   • 주파수 벡터가 정수가 아닐 경우 발생하는 양자화 오차 문제를 다룹니다.
   • 비정수 주파수 벡터를 인접한 정수 벡터로 대체하여 문제를 표준 정수 벡터 나머지 프레임워크로 매핑하는 근사 2D-DFT 피크 모델이 도입되었습니다.
   • 양자화 오차(무잡음 케이스)와 가법 신호 잡음을 모두 처리하기 위해 **강건한 MD-CRT(Robust MD-CRT)**가 통합되었습니다.
   • 저자들은 강건한 복구를 위한 충분 조건을 확립하고, 무잡음 및 잡음이 있는 설정 모두에 대해 명시적인 재구성 오차 경계를 도출했습니다. 이 경계는 스케일링 인자 $\alpha_j$와 잡음 수준에 따라 달라집니다.

주요 기여

• 통합 2차원 정식화: 본 논문은 평면 배열을 사용하여 크로스 레인지 시프트와 표적 높이를 공동으로 추정하기 위한 통합된 2차원 벡터 나머지 정식화를 확립하여, 기존의 1차원 VSAR 개념을 확장했습니다.
• 행렬 모듈러스 최적화: 오차 전파를 최소화하는 최적의 2D-DFT 행렬 모듈러스에 대한 이론적 도출을 제공하며, 이는 2차원 설정에서 고유한 구조적 설계 자유도를 제공합니다.
• SAR를 위한 MD-CRT 프레임워크: 서로 다른 격자 생성기를 가진 다수의 평면 서브어레이가 어떻게 공동으로 모호성을 해결하고 단일 배열보다 넓은 무모호 범위를 확보할 수 있는지 보여줌으로써 MD-CRT를 SAR 이미징에 적응시켰습니다.
• 강건성 분석: 근사 피크 모델을 도입하여 비정수 주파수에 대응하였으며, 양자화 및 가법 잡음이 존재하는 환경에서 강건한 복구를 위한 엄격한 충분 조건과 오차 경계를 제시하였습니다.

결과

논문은 이론적 프레임워크를 검증하기 위해 수치 시뮬레이션을 제시합니다:

• 무모호 범위 확장: 두 개의 서로 소(co-prime)인 평면 서브어레이를 사용한 시뮬레이션 결과, 결합된 무모와 범위(서브어레이 모듈러스의 LCRM으로 정의됨)가 단일 서브어레이의 범위보다 엄격히 더 넓음을 입증했습니다.
• 국지화 정확도: 무잡음 시나리오에서, 두 서브어레이 프레임워크는 도로 위의 표적 위치를 성공적으로 복구했습니다. 반면, 단일 서브어레이를 사용하면 모호성으로 인해 재구성된 표적 위치가 실제 경로에서 벗어나게 됩니다.
• 오차 경계: 도출된 오차 경계는 스케일링 인자 $\alpha$를 증가시키는 것(효과적으로 배열 크기를 키우거나 제로 패딩을 사용하는 것)이 재구성 오차를 줄여 정확도를 향상시킨다는 점을 확인시켜 줍니다.

의의 및 주장

저자들은 본 연구가 이동 표적 SAR 이미징을 위한 평면 배열 설계 및 성능 분석의 이론적 토대를 제공한다고 주장합니다.

• 이 프레임워크는 기존의 1차원 선형 배열으로는 달성할 수 없는 3차원 운동 파라미터(시프트 및 높이)의 공동 추정을 가능하게 합니다.
• 제안된 다중 서브어레이 접근 방식은 배열의 기하학적 구조를 넘어 복잡한 하드웨어 변경 없이도 무모한 범위를 확장할 수 있는 체계적인 방법을 제공합니다.
• 본 논문은 이 2부작 시리즈의 Part II에서 구체적인 평면 서브어레이 설계와 강건성에 대한 더 상세한 시뮬레이션을 조사할 것임을 명시하며, 이는 Part I이 일반적인 프레임워크, 수학적 정식화 및 이론적 경계에 집중하고 있음을 나타냅니다.

본 연구의 의의는 1차원 VSAT 기술과 3차원 지형 이미징의 더 복잡한 요구 사항 사이의 간극을 메우며, SAR 시스템을 위한 수학적으로 엄밀한 모호성 해결 및 오차 관리 접근 방식을 제공한다는 데 있습니다.