Ⅰ. 서 론
FMCW(Frequency Modualted Continuous Wave) 신호기반의 레이다 시스템은 기존의 펄스(chirp-pulse 신호기반) 레이다 시스템보다 소형/경량화에 여러 장점이 있어 최근 다양한 연구목적을 위해 활용되고 있다. 특히, 레이다 영상 획득을 위한 SAR(Synthetic Aperture Radar) 시스템으로의 활용도는 시스템적 측면의 장/단점을 고려해 다양한 레이다 영상화 신호처리 기법들이 연구되고 있다[1]~[4].
특히, FMCW 신호기반의 SAR 시스템은 동일한 변조주기에서 레이다 신호의 송·수신이 연속적으로 이루어지는 특징이 있어, 시스템 구성에 따라 상대적으로 송·수신 안테나의 연속이동효과(continuous motion effect)가 두드러지는 문제가 있다[5]. 이는 항공기 기반 시스템의 SAR 영상 복원을 위해 활용되는 back-projection 알고리즘을 이용한 영상 복원 시 영상 품질을 저하시키는 요인으로 작용하여, FMCW 신호 모델을 이용한 연속이동효과의 분석과 최적화된 영상 복원 기법 적용이 필요하다.
이를 위해 본 논문에서는 연속이동효과에 의한 레이다 수신 신호의 위상변화를 FMCW-SAR 신호 모델을 통해 분석하였다. 연속이동효과에 의한 위상변화는 시간-도플러(t, ku) 영역의 신호처리 기법에 의해 제거될 수 있으며, 영상 복원을 위한 거리압축신호에 적용된다. 그리고 SAR 영상 복원을 위해 적용된 back-projection 알고리즘에 연속이동효과 보정을 위한 수정된 정합필터를 적용하였으며, 그 결과를 실제 항공기 기반 SAR 시스템의 관측 데이터에 적용하여 연속이동효과 보정 결과를 검증하였다.
본 논문에서는 FMCW-SAR 신호 모델과 이를 이용한 연속이동효과 및 보정 결과를 모의 실험하고, 실제 영상 복원을 위해 back-projection 알고리즘의 수정된 정합필터 적용 사례를 II장 FMCW-SAR 신호 모델과 III장 연속이동효과 모의실험, 그리고 IV장 실험 결과에서 각각 설명한다.
Ⅱ. FMCW-SAR 신호 모델
FMCW-SAR 신호 모델은 시스템 특성을 분석하고, 레이다 원시신호 생성과 영상 복원 과정을 모의실험하기 위한 유용한 도구로써 활용될 수 있다. 특히, FMCW-SAR 시스템의 연속적인 안테나 위치 변화 특성을 그림 1과 같이 도식화할 수 있으며, 안테나 위치와 목표물 함수는 u와 f(x, y)로 각각 표현될 수 있다[6].
FMCW 신호 모델은 식 (1)과 같이 운용 주파수(f0)와 주파수 변조율(Kr = 대역폭/변조주기)을 이용해 톱니파형태의 주파수 변조 신호로 나타낼 수 있으며, 수신 신호는 목표물로부터 반사되는 지연 시간(τ)을 이용해 식 (2)와 같이 표현할 수 있다. 이때, 목표물 지연 시간은 식 (4)와 같이 목표물 위치(xi, yj)와 안테나 위치(u)의 거리(R)로부터 계산된다[7].
FMCW 신호기반의 레이다 시스템은 송·수신 신호를 주파수 혼합기를 통해 주파수 하향 변환하고, 두 신호의 차에 해당하는 비트 주파수(beat frequency) 성분을 샘플링하여 영상 복원 과정에 이용한다[7]. 이와 같은 과정을 식(3)으로 정리하여, 다음 과정에서 FMCW-SAR 시스템의 수신 신호로 이용된다(단, c는 광속 3×108 m/s이다.)[8].
위 식 (1)~(3)은 표현의 편의상 신호세기 성분이 생략되었으며, 모의실험을 위해 안테나 빔패턴, 목표물 거리와 반사도 등이 고려되어야 한다[9].
FMCW-SAR 시스템에서 고려되어야 할 연속이동효과는 레이다 신호 송·수신 과정에도 안테나의 위치가 지속적으로 이동하는 것으로, 기존에 안테나 위치를 정지된 위치로 근사화시킨 ‘start-stop’ 근사기법만으로 그 처리가 부족한 측면이 있다[5],[10]. 이를 보완하기 위해서 시간에 따른 안테나 위치 변화를 항공기 이동 속도(vr)와 측정 시간(sweep time)을 고려해 식 (5)와 같이 수정하였다.
식 (6)은 각각의 송신신호마다 정지된 안테나 위치로 표현된 기존 거리 함수 R(u)를 항공기 이동 속도와 시간 변화에 따른 안테나 위치 변화의 세부 정보를 추가한 2차원 배열의 거리함수 R(t, u)로 나타낼 수 있다.
식 (7)은 연속이동효과를 포함한 수정된 거리함수 R(t, u)를 이용해 거리성분에 대한 공간 주파수(kr = 2π/λ) 영역의 신호로 표현된 것이다. 이때, 안테나 위치 변화를 포함한 공간 주파수 영역의 수신 신호는 비트 주파수 스펙트럼으로 표현되는 거리 정보의 퓨리에 역변환 관계에 놓이며, 이는 시간영역과 동일한 영역의 신호임을 알 수 있다. 따라서 앞선 식 (3)은 시간 영역 신호에 연속이동효과를 더해 식 (7)과 같이 정리될 수 있다. 이를 시간-도플러 S(t, ku) 영역의 신호로 변환하는 과정은 식 (8)~(10)과 같다.
식 (8)은 도플러(ku) 영역으로의 퓨리에 변환 과정을 나타낸 것으로 변수치환(u‘ = u − vrt) 과정과 정지 위상 근사기법(stationary phase approximation)을 적용한 과정을 나타낸 것이며, 식 (10)의 관계식을 이용해 식 (8)의 결과식을 다음과 같이 얻을 수 있다.
위의 최종식은 항공기 이동 속도에 따른 위상 변화 부분(식 (11) 우변 첫 번째 항)과 도플러 영역의 목표물 위상함수로 구분된다.
Ⅲ. 연속이동효과 모의실험
항공기 기반 FMCW-SAR 시스템의 연속이동효과 특징을 분석하고, 이를 보정하기 위해 앞서 설명한 FMCW 신호 모델과 시간-도플러 S(t, ku) 영역 신호처리 기법을 적용한 모의 실험을 수행하였다. 영상 복원을 위해 back-projection 알고리즘을 이용하였으며, 복원 영상 내 연속이동효과를 효과적으로 표현하기 위해 실제보다 넓은 안테나 방사패턴을 갖도록 설정하였다(표 1 참고).
설정 항목 | 설정값 | 비고 |
---|---|---|
주파수 대역 | 10.0~10.5 GHz | 해상도 30 cm |
변조율(Kr) | 5e11 Hz/s | chirp-rate |
샘플링 주파수 | 1.2 MHz | 1,024 points |
측정 시간 주기 | 1 ms | PRF = 1,000 Hz |
안테나 방사 패턴 | 30° | 반전력빔폭(HPBW) |
FMCW-SAR 시스템에서 발생되는 연속이동효과는 목표물과 안테나 간의 거리변화를 나타낸 식 (6)으로써 그 특성을 확인할 수 있다.
그림 2는 목표물 지연 시간(τ(t, u)=2/c·R(t, u)) 특성을 시간에 따른 안테나 위치 변화를 모두 고려해 작성된 것으로 1 ms의 샘플링 주기 내 위치 변화를 모두 계산한 것으로 각 샘플링 시간과 그에 따른 안테나의 연속적인 위치 변화(1,024개)를 모두 그린 것으로 방위각 방향으로 중심점이 이동하는 것을 확인 할 수 있다.
그림 3은 앞선 측정 시간과 안테나 위치의 함수로 모델링된 목표물 지연 시간을 이용해 생성된 레이다 원시신호의 거리 압축 신호 s(w, u)를 모의실험한 것으로 연속이동효과 유·무를 비교한 결과이다.
식 (3)~(6)을 이용한 모의 실험을 통해 목표물 궤적이 방위각 방향으로 이동하는 FMCW-SAR 시스템의 연속이동효과를 확인하였다. 이와 같은 연속이동효과가 복원된 SAR 영상에 미치는 영향을 분석하기 위해서 back-projection 알고리즘을 이용해 레이다 원시 신호를 SAR 영상으로 복원하였다.
그림 4는 FMCW-SAR 시스템의 연속이동효과를 보정하기 위한 신호처리과정을 비교 설명한 것으로 ‘start-stop’ 근사기법(연속이동효과 미적용)으로 생성된 레이다 원시신호를 복원한 기준 SAR 영상과 연속이동효과를 적용한 레이다 원시 신호를 복원한 SAR 영상, 그리고 연속이동 효과 보정을 위한 시간-도플러 S(t, ku) 영역 신호처리기법이 적용된 SAR 영상을 각각 그림 4(a)∼(c)에서 나타낸다. 이때, 목표물의 분산 전력 특성 변화를 정밀하게 나타내기 위해 모의 실험 설정 해상도 30 cm보다 세밀한 해상도로 복원된 영상을 목표물 중심점을 기준으로 확대하였다.
연속 이동 효과를 적용한 그림 4(b) 레이다 원시 신호를 이용한 복원 영상은 기준 영상 대비 목표물의 분산 전력 패턴이 거리 방향으로 퍼지는 특성을 보이며, 방위각 방향으로는 항공기 이동 속도에 비례한 왜곡된 위상 성분으로 목표물 분산 전력 패턴이 왜곡되는 특성을 나타낸다.
항공기 이동 속도와 샘플링 시간에 비례하는 연속이동효과의 위상 변화 성분을 제거하기 위해서 앞서 설명한 식 (12)와 식 (13)의 시간-도플러 영역 신호처리 기법을 back-projection 알고리즘 입력 신호에 적용하여 그림 4(c)와 같은 연속이동효과 보정 결과를 얻었다. 시간-도플러 영역 신호처리 기법을 이용한 연속이동효과 보정 결과는 그림 4(c) 좌우 패턴과 같이 거리 방향의 목표물 분산 전력 패턴이 기준 영상의 패턴 특성과 매우 유사하게 보정됨을 확인 할 수 있다. 그러나 목표물의 중심점 부분의 영상 복원 개선 효과가 충분치 못해 이를 개선하기 위해 추가적으로 영상 복원을 위해 사용된 back-projection 알고리즘에 수정된 정합필터를 적용하였다.
식 (14)는 기존 back-projection 알고리즘의 정합필터 과정을 나타낸 것으로 복원 영상의 목표물 지연 시간(td)은 안테나 위치(u)와 복원 영상 내 목표물 즉 각 픽셀과의 거리를 이용해 계산될 수 있으며, u는 M개의 정지된 안테나 위치로 처리된다. 반면에, 연속이동효과 보정을 위해 수정된 정합 필터 과정은 식 (15) 및 식 (16)과 같이 시간(t)에 따른 변화를 추가한 목표물 지연 시간 특성 함수로 변형 될 수 있으며, 이는 기존의 back-projection 알고리즘의 안테나 위치 입력 정보(um)를 시간 샘플(tn) 개수(N)만큼 세분화하고, 추가된 안테나 위치정보 만큼 목표물의 분산전력이 누적되는 효과가 있다.
반면에 식(15)의 수정된 정합필터 과정은 시간샘플 개수 N에 비례한 안테나 세부 위치를 이용한 영상 복원 과정으로 기존 back-projection 알고리즘의 연산량을 N배 만큼 증가시키는 단점이 있으나, 이는 일정 비율로 시간 샘플 개수를 낮춰 계산함으로써 보정 정확도 유지와 계산 효율 저하를 보완 할 수 있다. 다시 말해, 연속이동효과 보정을 위해 수정된 정합필터는 기존 안테나 위치 입력 정보(um)에 항공기 이동 속도와 시간 변화를 고려한 세부위치 입력 정보(vr tn)를 더한 것으로 기존 back-projection 알고리즘을 이용한 영상 복원 과정에 약 5∼10배(예, N/10∼N/100 비율 모의 실험 적용) 증가된 입력 정보를 추가해 개선된 목표물 분산 전력 특성을 얻을 수 있다.
그림 5는 시간-도플러 s(t, ku) 영역에서 연속이동효과 위상성분을 보정한 레이다 거리압축 신호를 입력으로 한 back-projection 알고리즘의 영상 복원 결과이다. 수정된 정합 필터 적용으로 목표물의 분산전력 특성은 중심점에 더욱 집중된 결과를 확인할 수 있으며, 이를 거리방향과 방위각방향으로 각각 추출하여 기준 SAR 영상(REF)과 연속이동효과(CME) 보정 전/후의 결과로 비교한 것이다. 그 결과 시간-도플러 영역에서의 위상 변화 보정 기법(CMC)과 수정된 정합필터(MF)를 적용한 경우, 거리/방위각 방향 해상도 각각 보정 전 37.2 cm/6.5 cm, 보정 후 30.1 cm/8.9 cm로 분석되었으며, 거리방향 해상도 개선과 목표물 중심점 부근의 분산 전력 특성 개선효과가 있음을 확인하였다. 단, 방위각 방향의 해상도는 안테나 세부위치 정보를 추가한 수정된 정합필터로 인해 기준 신호 대비(8 cm) 약 0.9 cm 증가한 수치를 확인하였다.
Ⅳ. 실험 결과
FMCW-SAR 시스템 연속이동효과 보정 기법 성능을 비교/검증하기 위해 항공기 기반 FMCW-SAR 시스템을 이 용한 현장관측 자료를 활용하였다. 본 실험연구에 적용된 FMCW-SAR 시스템은 저고도/근거리 비행 실험에 최적화된 ImSAR社(모델명 NanoSAR-B, 최대출력 1 W)에서 제작된 시스템으로 운용 효율을 높이기 위해 delayed-dechirp 방식으로 운용된다[10]. 이는 상대적으로 낮은 샘플링 주파수(fs=1.2 MHz)와 적은 샘플개수(1,024 points)를 보다 효율적으로 사용하기 위한 방식으로 약 480 m의 상쇄 거리(range offset)를 포함하여 운용되었다(표 2 참고).
레이다 항목 | 제원특성 | 비고 |
---|---|---|
운용 주파수 | 10~10.5 GHz | 2채널 |
송신출력 | 최대 1 W | (최대 30 dBm) |
측정시간 | 약 30초 | PRF = 1,000 Hz |
샘플링 주파수 | 1.2 MHz | 샘플개수 1,024개 |
안테나 방사 패턴 | 12° | HH-pol.(수평편파) |
비행 고도/속도 | 425 m / 320 km/h | − |
연속 이동 효과 보정 기법 검증을 위해 사용된 관측자료는 2016년 1월 대부도 인근 매립지 주변에서 측정된 것으로 넓은 평지로 이루어진 지형에 시험용 목표물(삼각수동전파반사기, 크기 50 cm)을 설치하고, 그 특성을 복원 영상으로부터 분석하였다.
그림 6은 FMCW-SAR 시스템의 연속이동효과 보정성 능을 비교분석한 결과로 보정 전/후의 SAR 영상이며, 시험용 목표물 설치 지역 내 확대 영상으로 보정 성능을 분석하였다. 시험용 목표물을 이용한 보정 성능 분석 결과는 표 3과 같으며, 복원 영상 내 시험용 목표물의 분산 전력 특성은 그림 7과 같이 거리 방향 성분의 개선 효과를 확인 할 수 있었다. 이때, 이론적인 경사거리(slant-range) 해상도는 약 31 cm로 목표물 전파 입사각을 고려할 때 지면거리 해상도는 약 44 cm 정도가 되며, 윈도우 함수를 적용한 일반적인 PSLR(Peak-Sidelobe Ratio) 성능 –20 dB보다 낮은 결과를 확인하였다[9].
항 목 | 보정 전 | 보정 후 | 비고 |
---|---|---|---|
해상도 | 0.63 m | 0.51 m | 거리 방향 |
0.46 m | 0.47 m | 방위각 방향 | |
PSLR(peak sidelobe ratio) | −22.6 dB | −25.9 dB | 거리 방향 |
−20.4 dB | −20.1 dB | 방위각 방향 |
Ⅴ. 결 론
본 연구에서는 FMCW-SAR 시스템의 연속이동효과를 보정하기 위해 시간-도플러 s(t, ku) 영역 신호처리기법과 SAR 영상 복원을 위한 back-projection 알고리즘의 수정된 정합 필터 기법을 제안하였으며, 연속이동효과 보정검증을 위해 FMCW 신호 모델을 이용한 모의실험을 수행하였다. 또한, 제안된 연속이동효과 보정 기법의 성능검증을 위해 실제 항공기 기반 FMCW-SAR 원시자료를 이용하였으며, 분석결과 거리방향 해상도 개선(측정값 51 cm, 이론값 44 cm)과 PSLR 2∼3 dB 개선 효과를 확인하였다.