PDF Download
PubReader
Export Citation
Email the Author
Share Facebook
Share Twitter
Cited by 0 Article
Metrics (View:4,704)
Ⅰ. 서 론
SAR의 대표적인 운용 기법으로 spotlight, stripmap 모드가 있다. Spotlight 모드는 영상 획득 동안 SAR 안테나가 영상 중심점을 지향하도록 안테나 빔을 조향하며 획득하기에 합성 개구면 각도(synthetic aperture angle)를 방위 빔 폭에 상관없이 증가시킬 수 있으므로 방위 방향으로 고해상도 영상 획득이 가능한 반면, 영상의 방위 관측폭이 빔 폭에 의해서 제한된다[1]. 이와 달리 stripmap 모드는 안테나 방향은 비행체 진행 방향과 직교하는 측면 방향(side-looking)으로 고정시켜 놓고 영상을 획득하게 되므로 방위 방향으로 관측폭 제한 없이 영상 획득이 가능하나, 방위 해상도가 안테나 유효 빔 폭에 의해서 제한된다[2]. 따라서 고해상도 모드로 집중 감시가 필요한 지역의 영상 획득 시에는 spotlight 모드의 운용을 고려하며, 넓은 지역에 대해 정찰 목적의 영상 획득 시에는 상대적으로 해상도를 낮춘 stripmap 모드 운용을 고려한다. 항공기 SAR의 경우, 임무 중 바람, 진동 등에 의한 요동을 경험하게 되며, GPS/IMU와 같은 항법 장치에 의한 요동 보상 이후 잔여 오차를 추정하여 제거하는 데이터 기반의 자동 초점(autofocus, AF) 기법을 반드시 적용해야 한다.
Stripmap 모드 데이터를 spotlight 모드로 처리하기 위한 연구가 수행[3]～[6]되었으나, 어떻게 AF를 적용해야 최종 성능 개선이 가능한지까지 포함하지는 않았다. 본 논문에서는 항공기 탑재 SAR에서 stripmap 모드로 획득한 데이터에 대해서 spotlight SAR 모드 형태의 데이터 처리 및 AF 적용 기법을 제시한다. Stripmap 모드의 요구 방위 해상도 대비 충분히 넓은 빔 폭 가정이 타당할 경우[5],[6], stripmap 모드 알고리듬인 RDA를 이용한 spotlight SAR 처리 형태의 영상 형성 기법을 설명한다. 그리고 AF 전처리 함수 적용을 통해 기존의 AF 기법을 수정 없이 쉽게 적용할 수 있음을 보이며, 방위 방향 영상 품질 개선에 필수적임을 보인다. 그리고 spotlight SAR 처리한 영상 블록들의 방위 방향 병합을 통해 최종적으로 방위 방향으로 stripmap 모드의 연속 영상 형성이 가능함을 실제 비행시험을 통해 획득한 원시데이터를 처리하여 보이도록 한다.
Ⅱ. 본 론
Stripmap 모드 원시 데이터를 spotlight 모드로 처리하기 위해서는 방위 해상도와 관측폭 사이의 관계 및 기하학적 제한 사항을 반드시 고려해야 한다. 이론적으로 stripmap 모드에서 점표적의 방위 해상도는 그림 1처럼 방위 방향 최대 유효 빔 폭인 θeff으로만 제한되지만[2], 그림 2와 같이 요구 해상도 대비 실제 안테나의 유효 빔 폭이 충분히 큰 경우라면, 요구 방위 해상도 ρaz는 식 (1)처럼 θeff 대신 θsyn의 함수로 표현할 수 있다[1]～[3].
여기서 λ는 송신 주파수의 파장, θsyn는 합성 개구면 각도를 의미한다. θsyn는 θeff보다 항상 작게 되며, spotlight 모드로 처리 가능한 방위 관측폭 SWaz는 그림 2에서처럼 나타나며, 식 (2)로 표현할 수 있다[3].
작은 각도 근사(small angle approximation)로 위의 수식은 식 (3)과 같이 근사화해서 표현할 수 있다.
식 (3)에서 θeff과 θsyn이 같다면, 일반적인 stripmap 모드 형태가 되어 spotlight 모드로 처리 가능한 관측폭은 이론적으로 전혀 발생하지 않으므로 본 논문에서 고려하지 않는 경우가 된다. 즉, θsyn이 작아질수록, 즉 요구 해상도를 낮출수록, SWaz은 커지게 된다. 요구 해상도 만족을 위한 합성 개구면 길이(synthetic aperture length, SAL) 동안 비행체가 이동하면서 stripmap 모드로 운용했을 때, SWaz이내의 모든 표적들에 대해서는 마치 spotlight 모드로 획득한 데이터로 간주할 수 있다. SAL 대비 SWaz가 너무 작을 경우, 최종 방위 방향 영상 병합에 필요한 영상들이 많아져 연산량에 부담이 되므로, 요구 방위 해상도에 의해서 결정되는 SAL(≈SRc · θsyn) 대비 SWaz가 너무 작아지지 않도록 해상도와 관측폭 사이 trade-off를 체계 설계에 반드시 고려해야 한다.
RDA는 stripmap 모드의 대표적인 영상 형성 알고리듬으로서 기본적으로 거리 압축, 거리-도플러 영역에서 거리 셀 이주 보상 RCMC(range cell migration correction), 마지막으로 방위 압축을 통해 방위 시간 영역에서 영상을 형성하게 된다. 일반적인 stripmap 모드의 RDA 방위 압축 시에는 유효 빔 폭에 의한 방위 대역폭을 모두 고려하므로, 방위 위치가 다른 표적들이 경험하는 도플러 대역폭이 하나의 중심 주파수(스퀸트가 없을 때 0°)를 기준으로 동일하게 나타난다[2]. 다만 각 표적 별로 필요한 SAL의 시작 위치와 종료 위치가 방위 방향 표적 위치 별로 모두 달라지므로 각 표적 별로 다른 SAL 동안 경험하는 요동의 형태는 모두 달라지며, spotlight SAR와 달리 AF에 의한 성능 개선이 크게 저하되는 주요 원인이 된다[7]. 이와 달리 그림 2의 관측폭 SWaz 내의 모든 표적들은 spotlight 모드처럼 SAL 동안 동일한 요동을 경험한다. 즉, 관심 영역 SWaz에 대해서 SAL동안 RDA에 의한 영상 형성 시 고려 사항은 다음 2가지로 요약할 수 있다.
RDA의 방위 정합 필터는 다음 식 (4)로 모델링된다.
여기서 η는 방위 방향 시간, ηc는 방위 방향의 합성 개구면 시간(synthetic aperture time, SAT)의 중심, v는 비행체의 비행 속도, 은 경사거리 벡터를 나타낸다. 전체 유효 빔 폭 θeff을 고려해야 하는 stripmap RDA에서 SAL은 θeff에 의해서 결정되나, spotlight 처리 시에는 θsyn에 의한 SAL만을 고려해서 생성해야 한다.
RDA에 의한 SAR 영상 형성 시 RCMC 이후 항법 장치에 의한 요동 보상을 수행한다[8]. 이를 통해 요동에 의한 위치 오차는 거리 셀 이내의 크기로 감소시킬 수 있기에 무시가 가능하나, 방위 방향 잔여 위상 오차 성분은 영상의 해상도, 부엽 품질에 민감한 영향을 미치며, 큰 열화 요인이 된다. 품질을 개선하기 위한 AF 기법의 적용은 반드시 필요하며, 잘 알려진 PGA(phase gradient autofocus)의 경우, PFA(polar format algorithm)와 같은 방위 dechirp에 의한 영상 형성 알고리듬에 최초 제안된 알고리듬[1],[9]이기에 방위 펄스 압축을 통해 시간 영역에 영상을 형성하는 RDA에 바로 적용할 수 없다. 이를 극복하기 위해 RDA 처리 과정 중 방위 dechirp 신호를 만들어 위상 오차를 추정 및 보상하고, 방위 압축 시 적용하는 형태의 알고리듬이 제시되었으나, 매우 복잡하고 연산량이 많다는 단점이 있다[7]. 방위 시간 영역에 영상을 형성하는 2차원 주파수 영역 알고리듬인 RMA(range migration algorithm), 2차원 시간 영역 알고리듬인 BPA(back projection algorithm) 등에 의한 spotlight SAR 영상 형성 결과에 PGA를 적용하기 위한 연구가 이전에 수행되었다[10],[11]. Spotlight 모드에서 방위 방향으로 서로 떨어진 표적들의 도플러 대역은 중심 표적 위치 대비 각 표적의 방위 위치 차이에 의한 도플러 주파수 차이만큼 이격되어 나타난다. 즉, 각 표적의 방위 도플러 대역폭이 영상 중심 도플러 주파수를 중심으로 정렬이 되어야 AF에 의한 위상 오차 추정 및 보상이 가능하므로 AF 적용 전 도플러 전처리 연산 반드시 필요함을 입증하였다. 이를 통해 spotlight RDA 영상에서도 AF를 적용하기 위해서는 도플러 전처리 연산이 반드시 필요함을 추정할 수 있다. 스퀸트가 없다고 할 때, 도플러 위치 정렬함수 수식은 다음과 같이 나타난다[10],[11].
여기서 은 1×M 경사거리 벡터 , 는 N×1 방위 방향 영상의 위치 벡터이다. 도플러 위치 정렬 함수 ϕ는 N×M 행렬이 되며, RDA 영상 영역에 적용할 경우, 방위 주파수 영역에서 도플러 정렬이 수행된 결과를 얻을 수 있으며, 이를 그림 3에 나타내었다. 거리 빈 별 모든 표적들이 각 표적들의 방위 위치에 따른 도플러 주파수 대역위치에 나타나던 신호들이 도플러 위치 정렬 함수에 의해 방위 방향 도플러 중심 주파수(영상 중심 위치)를 기준으로 정렬되었음을 알 수 있다. 다만, 도플러 대역폭이 경사거리가 멀어질수록 감소하는 이유는 SAL 동안 거리 별로 다른 위치에 있는 표적들이 경험하는 도플러 대역폭이 거리에 따라 달라지기 때문이며, 도플러 대역폭은 식 (6)과 같다.
즉, 표적의 도플러 대역폭은 경사거리에 반비례하므로 경사 거리별 도플러 대역폭을 스케일링(scaling)하여야 서로 다른 거리에 위치한 표적들이 SAL 동안 경험한 동일 위상 오차 성분들을 정렬시킬 수 있으며, 이후 위상 오차를 추정해야 AF 적용 시 정확한 성능 개선을 이룰 수 있다. 거리 관측폭이 작을 경우, 경사 거리별 도플러 대역폭의 차이가 크지 않아 무시할 수 있는 정도라면 스케일링 없이 AF 적용을 고려할 수 있다. 그러나 거리 방향 관측폭이 커진다면 반드시 도플러 대역폭 스케일링 후 AF를 수행해야 정확한 위상 오차 추정이 가능하다. 그림 4는 그림 3(b)의 도플러 정렬 함수 적용 이후 결과에 대역폭 스케일링을 수행한 결과이다. 마치 모든 경사거리의 표적이 동일한 도플러 대역폭을 갖는 것처럼 스케일링되었음을 확인할 수 있다.
다음 그림 5는 도플러 위치 정렬함수와 대역폭 스케일링 함수를 적용한 SAR 영상이다. 영상에서 경사거리가 멀어질수록 방위 방향 영상 폭이 비례해서 감소한 것은 경사 거리가 멀어질수록 영상의 해상도가 낮아지므로(도플러 대역폭이 작아지므로) 영상 영역 zero-padding에 의한 도플러 영역 내삽(interpolation)으로 대역폭 스케일링을 수행하였기 때문이다. 전처리 연산이 끝난 SAR 영상에 대해서 기존 AF 기법을 수정 없이 적용할 경우, 완전한 성능 개선을 이룰 수 있게 된다. AF 적용 이후에는 도플러 영역에서 원래 표적들의 대역폭이 되도록 영상의 zero-padding 영역을 제거하는 리스케일링(re-scaling)을 반드시 적용해야 이후 영상 블록에 대한 매칭 기법 적용이 가능하다. 즉, SAL 동안 거리 방향으로 넓은 지역에 대한 영상 형성 시 근거리 지역보다 원거리 지역의 해상도가 점진적으로 저하되므로, 반드시 영상의 최대 거리에서 해상도 만족이 가능하도록 체계 설계를 수행해야 한다.
Spotlight RDA 영상에 AF를 적용하여 영상 품질을 개선하는 방식의 영상 형성을 방위 방향으로 지속 수행할 경우, 광역 stripmap 영상을 형성할 수 있으며, 이를 아래 그림 6에 나타내었다. 다만 영상 블록 별 중첩을 반드시 고려해서 블록을 나누어야 중첩 부분 영상의 상호 상관(cross-correlation) 이용한 방위 이동 편차 보상 매칭 기법을 적용[12]할 수 있다. 이를 통해 stripmap 모드 운용의 장점인 방위 방향 연속 영상 형성이 가능해진다. 본 논문에서 제안한 항공용 SAR에서의 spotlight RDA처리 기법을 이용한 stripmap SAR 영상 형성 기법에 대한 블록 다이어그램을 기존(conventional) RDA 처리 과정을 포함하여 그림 7에 자세하게 나타내었다.
Ⅲ. 실험 결과 분석
본 논문에서 제안한 처리 기법을 자체 비행 시험을 통해 확보한 SAR 원시데이터를 활용[13]하여 검증하기로 한다. Stripmap 모드의 미터 급 요구 해상도를 위한 θsyn ≈ 1° 대비 방위 유효 빔 폭 θeff는 약 3° 이상으로 충분히 크기에 stripmap 데이터의 spotlight RDA 영상 형성 및 AF 적용이 가능하도록 체계 설계를 수행하였다. Stripmap SAR 영상 획득의 주요 운용 변수는 비행 고도 3.04km, 경사거리 약 20.59km, 스퀸트 각도 0°, 전체 영상의 SAT는 약 57초, 각 영상 블록의 SAT는 약 9초이다. 스퀸트가 발생하지 않도록 제어하였고, 도플러 중심 주파수가 0 Hz에 가깝게 되었기에 제로 도플러 변위를 위한 처리는 필요하지 않았다. 방위 방향으로 8개의 spotlight SAR 영상 블록이 되도록 중첩하여 처리하고, 방위 영상 매칭 기법을 적용하여 최종적으로 stripmap 모드 형태의 SAR 영상을 형성하였고, 전체 SAR 영상의 크기는 거리 방향 약 6.236 km, 방위 방향 약 2.25 km이다. RDA에 의한 spotlight SAR 영상 형성 방식은 앞서 언급한 2가지 고려 사항을 반영하여 처리한다면 그 외에는 일반적으로 잘 알려진 RDA 기법과 동일하다. 2가지 고려 사항에 대한 수행 과정은 다음과 같으며, 그림 8에 자세히 나타내었다. SAR 영상에서 θeff 빔 폭에 의한 aliasing을 회피하려면 방위 압축 수행 전 영상 블록의 방위 방향으로 충분한 zero- padding 수행이 필요함을 강조하였다. 여기서는 θeff에 의한 방위 방향 샘플 수에서 θsyn에 의한 방위 방향 샘플수를 빼고 난 샘플 수만큼 zero-padding을 수행한 결과를 그림 8(a) 나타내었다. 그리고 식 (4)의 정합필터를 생성하여 블록 데이터 방위 압축을 수행하여 SAR 영상을 형성한 결과를 그림 8(b)에 나타내었다. 실제 획득한 원시 데이터에서는 θeff보다 큰 범위에서도 aliasing에 영향을 주는 신호가 유입될 수도 있으며, 그림 8(b)에서 보면 방위 방향 위 아래로 aliasing이 발생하였음을 확인할 수 있다. 그러나 방위 압축이 끝난 SAR 영상에서는 그림 8(c)에서처럼 오직 θsyn에 의한 중심 부근 SAR 영상 부분만 필요하며, 이 부분 그림 8(d)만 aliasing 열화를 회피하며 되기에 논문에서 제안한 처리 기법 적용 이후 열화가 없음을 확인하였다.
본 논문의 다음 실험 결과는 spotlight RDA에 의한 영상 형성 이후 AF를 적용하는 데 필수적인 도플러 전처리 함수 및 도플러 대역폭 스케일링 함수의 적용 전/후 결과를 단계적으로 비교 분석한다. 이를 통해 논문에서 제안한 spotlight RDA에 의한 영상 형성 및 AF 적용 기법의 유효성을 확인하는데 집중하기로 한다.
먼저 영상 획득 동안 비행체가 경험한 요동을 거리, 고도 방향으로 분석한 결과를 그림 9에 나타내었다.
첫 번째 영상 블록에 대해서 AF를 적용하지 않은 SAR 영상 형성 결과(a), 도플러 위치 정렬 함수/대역폭 스케일링 함수 적용 없이 PGA를 수행한 결과(b), 도플러 위치 정렬 함수만을 적용하고 대역폭 스케일링 함수는 적용하지 않고 AF를 수행한 결과(c), 도플러 위치 정렬 함수/대역폭 스케일링 함수를 모두 적용하고, AF를 수행한 결과(d)를 그림 10에 나타내었다. (a)를 보면 GPS/IMU에 의한 요동 보상을 수행한 결과이지만, 그림 9의 요동을 경험하면서 보상되지 않은 잔여 위상 오차에 의한 영상 품질 열화가 크게 발생한 것을 알 수 있다. (b) 영상을 통해 도플러 위치 정렬 함수 및 대역폭 스케일링 함수 적용 없이 PGA 바로 적용 시에, 영상의 품질이 (a) 영상 대비 일부 개선된 것처럼 보이나, 실제 SAR 영상의 활용성 측면에서 본다면 유의미한 개선이 있다고 보기 어렵다. (c) 영상은 (b) 영상 대비 개선이 더 이루어졌음을 확인할 수 있으나, 여전히 방위 해상도 및 부엽 품질 측면에서 크게 저하되어 있음을 확인할 수 있다. 논문에서 제안한 전처리 기법들을 모두 적용하고, PGA를 수행한 (d) 영상에서 성능 개선이 크게 이루어졌음을 쉽게 할 수 있다. 정량적인 분석을 위해 영상 내 서로 다른 위치 존재하는 임의의
밝은 표적 두 군데(1번 표적 : 좌측 상단, 2번 표적: 영상 중심)에 대한 성능 평가를 수행하였다. 코너 초점기(corner reflector, CR)가 없는 영상이라 해상도를 포함한 엄밀한 성능 분석은 불가하였고, PSLR(peak side lobe ratio), ISLr(integrated side lobe ratio)등의 방위 방향 부엽 성능 측정 결과를 표 1에 나타내었다.
(a) 영상의 품질은 심각하게 저하되어 1, 2번 표적을 식별할 수 없었기에 비교에서 제외하였다. 먼저 (b), (c) 영상 대비 (d) 영상의 부엽 성능이 크게 개선된 결과로부터 제안 기법의 유효함을 쉽게 확인할 수 있었다. 특히 (b) 대비 (c) 영상에서 부엽 성능이 상대적으로 큰 개선이 이뤄졌음을 확인할 수 있다. (c)의 경우 (d) 대비 해상도 측면에서 큰 열화가 발생하였음을 SAR 영상을 통해 확인할 수 있는데, 즉 해상도가 저하로 인해 상대적으로 주엽이 넓어지다 보니 부엽 성능이 영상의 품질 대비 상대적으로 양호하게 측정되었다고 볼 수 있다.
최종적으로 8개의 영상 블록에 대해서 방위 매칭 기법을 통해 stripmap SAR 영상 형성한 결과를 그림 11에 나타내었다.
Ⅳ. 결 론
Stripmap 모드에서 획득에서 SAR 원시데이터에 대해서 spotlight RDA의한 영상 형성 및 AF 적용 기법에 대한 연구를 수행하였다. 대표적인 stripmap 모드 영상 형성 알고리듬인 RDA을 적용하여 spotlight 모드 영상 형성 시 반드시 필요한 고려 사항을 제시하였다. 그리고 RDA에 의한 spotlight 영상 형성 결과에 AF를 적용하기 위해서는 도플러 위치 정렬 함수 및 대역폭 스케일링 처리가 필요함을 식별하고, 실제 비행시험을 데이터를 통해 제안 처리 기법의 유효성을 입증하였다. 광역의 방위 빔 폭과 방위 요구 해상도의 trade-off를 통해 stripmap 모드 영상 형성 시에도 AF에 의한 성능 개선에 매우 유리한 spotlight RDA 영상의 블록 매칭 기법 적용을 통해 광역의 영상 형성이 가능함을 확인하였다.
