논문/REGULAR PAPERS

SAR-GMTI에서 지상이동표적의 속도 추정 기법

배창식1, 전현무1, 양동혁1, 양훈기1,*
Chang-Sik Bae1, Hyeon-Mu Jeon1, Dong-Hyeuk Yang1, Hoon-Gee Yang1,*
Author Information & Copyright
1광운대학교 전파공학과
1Department of Radio Sciences and Engineering, Kwangwoon University
*Corresponding Author: Hoon-Gee Yang (e-mail: reply1988@kw.ac.kr)

© Copyright 2019 The Korean Institute of Electromagnetic Engineering and Science. This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: Dec 12, 2016; Revised: Jan 12, 2017; Accepted: Jan 25, 2017

Published Online: Feb 28, 2017

요약

본 논문에서는 2채널 displaced phase center antenna(DPCA) 기반의 SAR-GMTI 시스템에서 지상이동 표적의 속도 정보를 추정하는 알고리즘을 제안한다. 제시된 알고리즘에서는 이동표적의 across-track 속도는 기존의 along-track interferometry(ATI) 기법을 이용하여 추정이 가능하다고 가정한 후 표적의 along-track 속도를 추정하는 방법을 제시한다. 이를 위해 기존 이동표적 구조를 변형하여 이동표적 속도를 영으로 만들고, 이를 레이다 속도에 반영하여 레이다 속도가 변화된 새로운 기하학적 구조를 얻는다. 이후 합성 개구면(synthetic aperture) 내의 subaperture 신호에 대한 푸리에 변환을 통해 공간 주파수 중심점 위치를 이용하여 표적의 along-track 속도를 추정한다. 시뮬레이션을 통해 제안한 알고리즘의 속도 추정 성능을 검증하고, 또한 추정된 속도가 보상되어 이동표적 영상의 해상도 및 SINR이 개선됨을 보인다.

Abstract

A ground moving target’s velocity estimation algorithm applicable for a SAR-GMTI system using 2 channel displaced phase center antenna(DPCA) is proposed. In this algorithm, we assume target’s across-track velocity can be estimated by along-track interferometry (ATI) and present a method to estimate target’s along-track velocity. To accomplish this method, we first transform a radar-target geometry in which a moving target has zero velocity via altering a radar velocity such that target’s velocity is reflected into it and next manipulate the spectral centers of the subapertures within the synthetic aperture. The validity of the proposed algorithm is demonstrated through simulation results showing the performance of the target’s velocity estimation and the enhancement of reconstructed target image quality in terms of resolution and SINR.

Keywords: SAR; GMTI; SAR-GMTI; Along-Track Velocity Estimation; Subaperture

Ⅰ. 서 론

현재 민수와 방위사업의 다양한 분야에서 지상이동표적을 관측하고자 하는 수요가 늘어나면서 SAR 시스템에 ground moving target indication (GMTI) 기능을 접목시킨 SAR-GMTI 시스템이 주목을 받고 있다[1],[2].

SAR-GMTI 시스템에서는 정지된 지형 이미징 외에 지상이동 표적의 속도로 인한 SAR 영상의 왜곡을 보상하기 위해 표적의 속도를 정확히 추정하는 것이 중요한 과제로 여겨지고 있다. 표적의 across-track 속도의 경우, 각 채널에 수신된 신호의 위상을 이용하는 along-track interferometry(ATI) 기법[3]으로 간단히 추정이 가능하다.

Along-track 속도 추정은 linear frequency modulated (LFM) 형태의 azimuth 신호를 추정하는 과정과 동일하기 때문에 정합필터뱅크(matched filter bank), fractional Fourier transform(FrFT), 그리고 time frequency analysis와 같이 기존에 사용되는 LFM 신호 추정 기법들이 사용되고 있다. 대부분의 기법들은 높은 계산량과 구현상의 문제들로 실제 SAR-GMTI 시스템에 적용되기 힘든 점이 있다. 따라서 알고리즘을 개선시키거나 새로운 접근 방법을 통한 속도 추정 알고리즘 연구[4][9]가 최근에 많이 진행되고 있다.

계산량을 줄이기 위한 방법으로 실제 도로 환경에 대한 사전 지식(a priori knowledge)을 이용한 연구[6]와 second order Wigner-Ville distributuion(WVD)을 통한 time-frequency analysis 기법[7]을 이용한 연구 등이 진행되어왔다. 또한, multiple signal classification(MUSIC), Capon 기법을 이용하여 속도를 추정[8]하거나 segmented keystone transform(SKT)과 Doppler Lv’s transform을 결합한 SKT-DLVT[9]와 같은 새로운 방법에 대한 연구들도 계속해서 보고되고 있다.

본 논문에서는 합성 개구면(synthetic aperture)내의 subaperture 신호에 대한 푸리에 변환을 통해 공간 주파수 도메인에서 spectral center 위치 정보를 이용하여 표적의 along-track 속도를 추정하는 새로운 알고리즘을 제안한다. 이를 위해 2장에서는 SAR-GMTI 시스템 동작 개요를 설명하고, 3장에서는 표적의 속도 추정 알고리즘 및 이를 적용한 이동표적 SAR 이미징 방법을 제시하며, 4장에서는 시뮬레이션을 통해 제안한 알고리즘을 조사한다.

Ⅱ. 2채널 SAR-GMTI 시스템

지상이동표적에 대한 SAR 기하학적 구조는 그림 1과 같이 나타낼 수 있다. 그림 1에서 플랫폼은 y축 방향(azimuth)으로 vp = (0, vpy)의 속도로 이동하고 있고, 표적은 vt = (vtx, vty)의 속도로 이동하고 있다. Azimuth 방향의 시간, 즉 slow-time t가 0일 때 플랫폼은 (0, 0)에 위치해 있고, 표적은 (r0, 0)에 위치해 있다. 플랫폼이 이동하면서 표적으로부터 수신한 azimuth 신호의 위상 ϕa는 다음과 같이 표현할 수 있다.

ϕ a ( t ) 4 π λ ( r 0 + v t x t + 1 2 r 0 ( v t y v p ) 2 t 2 )
(1)
jkiees-28-2-139-g1
그림 1. | Fig. 1. 지상이동표적에 대한 SAR 기하학적 구조 | SAR geometry for ground moving target.
Download Original Figure

여기서 λ는 파장이다. 식 (1)을 시간에 대해 미분하면 도플러 주파수를 얻을 수 있으며, 이는 정지표적 신호의 것과 비교하여 DC 성분의 주파수가 존재하고, 변형된 기울기의 일차 주파수(linear frequency) 성분을 가짐을 확인할 수 있다. 이러한 성분들로 인하여 신호가 azimuth 방향의 정합필터를 통과함에 있어서 부정합(mismatch)이 발생하고, 그 결과 SAR 영상의 왜곡(퍼짐 및 이격)이 발생한다. 따라서 SAR 영상의 왜곡 현상을 보상하기 위한 표적의 속도 추정이 필요하다.

그림 2는 DPCA를 위한 2채널 SAR-GMTI 시스템 구조를 나타낸 것이다. DPCA는 이동표적의 속도를 추정하기 앞서 클러터 환경에서 이동표적을 탐지하기 위해 수행되며, DPCA 조건을 만족하기 위해 플랫폼은 t1에서 채널 1의 effective phase center와 t2에서 채널 2의 effective phase center가 일치하도록 이동한다. 이후 pulse repetition interval(PRI)이 Δt = t2t1이 되도록 설정하면 시스템은 DPCA 조건을 만족하게 된다.

jkiees-28-2-139-g2
그림 2. | Fig. 2. 2채널 SAR-GMTI 시스템 | 2-channel SAR-GMTI system.
Download Original Figure

Ⅲ. 표적 속도 추정

3-1 Across-Track 속도 추정

표적의 across-track(x축) 속도는 기존에 널리 사용되는 방법인 ATI 기법을 이용하여 추정할 수 있다고 가정하였다. 플랫폼이 2채널 구조로 DPCA 조건을 만족하는 경우, 각 채널에 수신된 신호간의 ATI를 적용시킬 수 있다. ATI를 수행한 신호의 위상, 즉 ATI 위상 Δϕ식 (2)와 같이 나타낼 수 있으며, 이때의 위상을 통해 표적의 across-track 속도 vtx를 추정할 수 있다.

Δ ϕ 4 π λ v t x Δ t
(2)
3-2 Along-Track 속도 추정

그림 3그림 1의 이동표적을 고정시킨 후 플랫폼의 상대적 움직임을 이용하여 기하학적 구조를 나타낸 것이다. 이때 레이다에 수신된 신호는 그림 1의 환경과 동일하다. 각각 vt = (vtx, vty), vp = (0, vpy)의 속도로 이동하던 표적과 플랫폼은 새롭게 변형된 기하학적 구조(그림 3)에서 v′t = (0, 0), v′p = (v′px, v′py)의 새로운 상대속도로 이동하는 것으로 생각할 수 있다. 이때 v′p은 이동표적의 속도와 플랫폼의 속도의 벡터적 합으로, 다음과 같이 표현할 수 있다.

v p = ( v p x , v p y ) = ( v t x , v p y v t y )
(3)
jkiees-28-2-139-g3
그림 3. | Fig. 3. 변형된 지상이동표적에 대한 SAR 기하학적 구조 | Transformed SAR geometry for ground moving target.
Download Original Figure

즉, 표적을 고정시킨 변형된 기하학적 구조는 플랫폼과 빔의 각도가 θ만큼 틀어지고(squinted), 플랫폼과 표적간의 가장 가까운 거리(표적의 위치에서 플랫폼의 이동궤적으로 내린 수선의 발)에 해당하는 지점이 Δx만큼 이격된 형태의 새로운 기하학적 구조를 가지는 것으로 생각할 수 있다. 이때 틀어진 각도 θ와 이격된 거리 Δx는 다음과 같이 구할 수 있다.

θ = tan 1 ( x p x v p y ) Δ x = r 0 sin θ
(4)

여기서 r0t=0일 때의 플랫폼과 표적과의 거리이다. 플랫폼의 이동궤적의 중심으로부터 Δx만큼 이격된 지점은 표적과의 가장 가까운 거리임과 동시에 영 도플러(zero Doppler) 주파수를 갖는 지점이기도 하다. 따라서 Δx는 다음과 같은 관계식을 통해 나타낼 수 있다.

Δ x = v p N T
(5)

여기서 Nt = 0일 때 플랫폼과 표적에 의해 발생하는 도플러 주파수로부터 제로 도플러가 되는 지점인 t = t0까지의 azimuth 샘플(sample) 개수로, 이는 다음(그림 4)과 같은 방법으로 구할 수 있다.

jkiees-28-2-139-g4
그림 4. | Fig. 4. Subaperture 푸리에 변환을 이용한 도플러 주파수 및 파라미터 N계산 과정 | Subaperture Fourier transform processing to acquire Doppler frequency and parameter N.
Download Original Figure

먼저 t = 0을 중심으로 유한한 크기의 subaperture 내의 azimuth 샘플에 대해 푸리에 변환(Fourier transform)을 수행하여 표적의 도플러 주파수, 즉 spectral center를 구한다. 이때 spectral center는 spectrum의 peak에 해당하고, 표적이 다수의 scatterer로 구성되는 환경에서는 peak 근처에서 1차 moment를 이용하여 구할 수 있다. 이후 subaperture의 중심점(t)을 이동해가며 각 subaperture의 spectral center를 조사하여 그 값이 영 도플러 주파수를 가질 때까지 이 과정을 반복하고, 이때의 중심점에 해당하는 인덱스와 최초 중심점 인덱스의 샘플 간격으로 N을 구한다. 각 샘플은 PRI 간격으로 샘플링이 되어 있으므로 NT는 플랫폼이 Δx만큼 이동하는데 걸린 시간 t0가 된다. 지상의 고정 클러터 성분은 DPCA 과정으로 거의 제거가 가능하여 spectral center를 추정하는데 큰 어려움이 없지만 nonstationary 클러터가 존재하는 환경에서는 spectral center의 추정이 어려울 수 있다. 이 경우에는 subaperture의 크기를 줄임으로써 spectral center에 위치한 피크를 sharp하게 할 수 있으나, 이로 인해 공간 주파수의 해상도가 떨어지는 단점이 있다. 또한, subaperture의 크기는 푸리에 변환의 계산량과 직접적인 관계가 있으므로 알고리즘의 계산량은 subaperture의 크기에 따라 달라진다.

다음으로 식 (4)(5)v′p에 대한 식으로 정리하면 다음과 같다.

v p = r 0 N T sin θ = r 0 N T sin ( tan 1 ( | v p x | v p y ) ) = ( v p x ) 2 + ( v p y ) 2
(6)

여기서 |v′px|은 Δx > 0인 조건을 만족하기 위해 사용된 것으로 v′px의 값을 알고 있는 환경이므로 가능하다. a = |v′px|, b = r0/NT, x = v′py로 치환하여 식을 정리하면 다음과 같은 x에 대한 4차 방정식을 얻을 수 있다.

x 4 + 2 a 2 x 2 + ( a 4 a 2 b 2 ) = 0
(7)

이때 식 (7)에 대한 해 x는 다음과 같다.

x 1 , 2 = ± 2 a 2 + 4 a 4 4 ( a 4 a 2 b 2 ) 2 x 3 , 4 = ± 2 a 2 4 a 4 4 ( a 4 a 2 b 2 ) 2
(8)

지상이동표적의 along-track 속도 vtyvpy보다 작다고 가정한다면, 즉 v′py = vpy-vty > 0이라면, 이때 x는 다음과 같은 단일 해를 가지게 된다.

x = ± 2 a 2 + ( 2 a 2 ) 2 4 ( a 4 a 2 b 2 ) 2
(9)

최종적으로 식 (9)를 풀어쓰면 v′py에 대한 식 (10)을 얻을 수 있다. 따라서 t = 0일 때의 플랫폼과 표적과의 거리 r0, PRI T, 영 도플러까지의 인덱스 개수 N, 그리고 표적의 across-track 속도 vtx(=−v′px)를 이용하여 식 (10)을 계산하고, 이를 식 (3)에 대입하면 표적의 along-track 속도를 추정할 수 있게 된다.

v p y = 2 | v p x | 2 + 4 | v p x | 4 4 ( | v p x | 4 | v p x | 2 ( r 0 N T ) 2 ) 2
(10)
3-3 SAR 영상 보상

추정한 이동표적의 속도가 (vtx^,vty^)라면 이를 이용하여 생성한 이동표적에 대한 azimuth 정합필터는 다음과 같다.

exp [ j ( 2 π 2 v t x ^ λ t + π 2 ( v t y ^ v p y ) 2 λ t 2 ) ]
(11)

식 (11)을 적용하여 이동표적에 대한 위치 왜곡 및 해상도, signal to interference plus noise ratio(SINR)이 개선된 복원 영상 결과는 Ⅳ장에서 제시한다.

Ⅳ. 시뮬레이션

본 절에서는 속도 보정이 이루어지지 않은 SAR-GMTI의 성능을 보이고, 이어서 표적 속도 추정 및 이에 따른 SAR-GMTI 성능을 비교하였다.

본 논문에서 표적의 속도 추정 시뮬레이션에서 사용한 파라미터는 표 1과 같다.

표 1. | Table 1. 시뮬레이션 파라미터 | Simulation parameters.
Parameter Value
Operating frequency 5.3 GHz
Pulse repetition interval 2.5 ms
Range chirp slope(kr) 250 GHz/s
Pulse duration 25 us
Synthetic aperture length(L) 600 m
Platform velocity(vp) (0, 150 m/s)
Target velocity(vt) (2 m/s, 5 m/s)
Target location @ t = 0 (20,300 m, 0)
Download Excel Table

그림 5는 DPCA 방법의 GMTI를 통한 클러터 억제 및 이동표적의 탐지 성능 개선을 보인 것이다. 클러터의 경우, 지형 환경에 따라 다양한 분포로 모델링이 가능하므로 본 논문에서는 모든 셀에 걸쳐 가우시안 분포로 존재한다고 가정하였다. 또한, 시스템의 물리적 오차 및 바람과 같은 미세한 변화로 인한 동일한 장면간의 decorrelation을 나타내기 위해 1이 아닌 0.8의 상관인수(correlation factor)를 가지도록 생성하였다. 그림 5(a)는 GMTI를 수행하지 않은 이동표적에 대한 SAR 영상으로, 클러터로 인해 이동표적의 식별이 힘든 것을 확인할 수 있다. 그림 5(b)는 GMTI를 수행한 후의 결과로 클러터 레벨이 상당히 낮아지고, 그에 따라 이동표적이 뚜렷이 나타남을 확인할 수 있다. 하지만 그림 5(b)는 정지표적 파라미터를 이용하여 생성한 azimuth 정합필터를 통과한 결과로, 이동표적의 t = 0에서의 위치(20,300 m, 0)와 다른 곳에 이미지가 형성되었고, 또한 정합필터의 부정합으로 인하여 표적의 영상이 퍼지고, 최대 피크가 낮아짐을 확인할 수 있다.

jkiees-28-2-139-g5
그림 5. | Fig. 5. 지상이동표적 SAR 영상 | SAR image for ground moving target.
Download Original Figure

이동표적의 속도 추정을 위해서 50개의 azimuth 샘플을 갖는 subaperture를 이용하여 spectral center를 조사하였다. 그림 6(a)t = 0일 때 공간 주파수 축에서 표적 스펙트럼 분포를 나타내며, 그림 6(b)t = t0에서의 표적 스펙트럼으로 영을 중심으로 위치해 있는 것을 확인할 수 있다. 이때 N은 772개로 이 N식 (10)에 대입시켜 변형된 기하학적 구조에서 플랫폼의 along-track 속도 v^py를 구하였고, 이를 식 (3)에 대입하여 표적의 along-track 속도 vty^를 구하여 표 3의 결과를 얻었다.

jkiees-28-2-139-g6
그림 6. | Fig. 6. 도플러 주파수 | Doppler frequency.
Download Original Figure
표 2. | Table 2. 변형된 geometry에서 속도 파라미터 | Velocity parameters on transformed geometry.
Parameter Value
Platform velocity(vp) (−2 m/s, 145 m/s)
Target velocity(vt) (0, 0)
Download Excel Table
표 3. | Table 3. 시뮬레이션 결과 - 파라미터 | Simulation results - parameters.
Parameter Value
Estimated platform along-track velocity(v^py) 145.11 m/s
Estimated target along-track velocity(vty^) 4.89 m/s
Download Excel Table

그림 7은 환경 변화에 따른 제안한 알고리즘의 추정성능을 확인하고자 상관인수를 바꾸어가며 평균 속도 추정 오차를 분석한 결과이다. 상관인수가 낮은 환경, 즉 두 장면간의 decorrelation이 강한 환경에서는 속도 추정 오차가 상당히 큰 것을 확인할 수 있다. 이는 상관인수가 낮을수록 DPCA를 통한 클러터 제거 성능열화가 발생하고, 이로 인해 표적의 spectral center를 정확히 추정할 수 없기 때문이다.

jkiees-28-2-139-g7
그림 7. | Fig. 7. 상관인수에 따른 속도 추정 오차 | Error of estimated velocity with respect to correlation factor.
Download Original Figure

그림 8표 3의 추정된 속도를 바탕으로 생성한 azimuth 정합필터를 통과한 SAR 영상의 결과를 나타낸 것이다. 이격된 표적이 원하는 위치에 나타나는 것을 확인할 수 있고, 그림 5(b)와 비교하였을 때 표적의 해상도가 향상됨을 확인할 수 있다. 또한, 표적의 최대 피크가 증가함에 따라 셀의 SINR이 그림 5(b)의 27.1 dB에서 42.6 dB로 15.5 dB 증가함을 확인할 수 있다. 보상된 SAR 영상은 추후 이동표적의 식별 및 시각화를 위한 용도로 사용될 수 있다.

jkiees-28-2-139-g8
그림 8. | Fig. 8. 추정된 속도로 보상한 지상이동표적 SAR 영상 | Compensated SAR image for ground moving target using estimated velocities after GMTI.
Download Original Figure

Ⅴ. 결 론

본 논문에서는 2채널 SAR-GMTI 시스템에서 표적의 along-track 속도를 추정하는 알고리즘을 제안하였다. 이동하는 표적을 고정시킨 기하학적 구조에서 플랫폼의 움직임을 관측하고, 이로 인한 이동표적 SAR 영상의 왜곡 특성을 분석하였다. Azimuth 샘플들의 subaperture 푸리에 변환을 통하여 공간 주파수를 획득하고, 이와 관련된 파라미터들을 이용하여 이동표적의 along-track 속도를 구할 수 있는 관계식을 유도하였다. 추정된 속도(4.89 m/s)는 실제 속도(5 m/s)와 비교하여 약 2 %의 오차를 가지며 잘 추정된 것을 확인할 수 있다. 추정된 속도를 바탕으로 azimuth 정합필터를 설계하여 온전한 이동표적 SAR 영상을 복원한 결과 영상의 해상도가 향상되고 SINR이 증가함을 확인할 수 있었다.

Acknowledgements

이 연구는 2016학년도 광운대학교 교내 연구비의 지원으로 연구되었음.

References

[1].

D. Cerutti-Maori, I. Sikaneta, and C. H. Gierull, “Optimum SAR/GMTI processing and its application to the radar satellite RADARSAT-2 for traffic monitoring”, IEEE Trans. Geosci. Remote Sens., vol. 50, pp. 3868-3881, Oct. 2012.

[2].

S. Suchandt, H. Runge, H. Breit, U. Steinbrecher, A. Kotenkov, and U. Balss, “Automatic extraction of traffic flows using TerraSAR-X along-track interferometry”, IEEE Trans. Geosci. Remote Sens., vol. 48, no. 2, pp. 807-819, Feb. 2010.

[3].

S. V. Baumgartner, G. Krieger, “Multi-channel SAR for ground moving target indication”, in Academic Press Library in Signal Processing: Communications and Radar Signal Processing, vol. 2, 1st ed., R. Chellappa and S. Theodoridis, Eds. New York, NY, USA: Academic, 2014, ch. 18, pp. 911-986.

[4].

S. V. Baumgartner, G. Krieger, “SAR traffic monitoring using time-frequency analysis for detection and parameter estimation”, in Proc. IEEE Int. Geosci. Remote Sense. Symp., Boston, MA, USA, Jul. 2008, pp. 25-28.

[5].

J. Wu, Y. Jiang, G. Kuang, J. Lu, and Z. Li, “Parameter estimation for SAR moving target detection using fractional Fourier transform”, in Proc. IEEE Int. Geosci. Remote Sense. Symp., Quebec City, QC, Canada, Jul. 2014, pp. 596-599.

[6].

S. V. Baumgartner, G. Krieger, “Fast GMTI algorithm for traffic monitoring based on a priori knowledge”, IEEE Trans. Geosci. Remote Sens., vol. 50, no. 11, pp. 4626-4641, Nov. 2012.

[7].

P. Huang, G. Liao, Z. Yang, X. Xia, J. Ma, and J. Zheng, “Ground maneuvering target imaging and high-order motion parameter estimation based on second-order keystone and generalized Hough-HAF transform”, IEEE Trans. Geosci. Remote Sens., vol. 55, no. 1, pp. 320-335, Jan. 2017.

[8].

J. Wen, M. Bao, and S. Lin, “Estimating radial velocity for ground moving targets from multi-channel SAR by spectral estimation”, in Proc. IET Int. Radar Conf., Apr. 14-16, 2013, pp. 1-4.

[9].

J. Tian, W. Cui, and X. Xia, “Parameter estimation of ground moving targets based on SKT-DLVT processing”, IEEE Trans. Comput. Imag., vol. 2, no. 1, pp. 13-26, Mar. 2016.

Author Information

배 창 식

jkiees-28-2-139-i1

  • 2015년 2월: 광운대학교 전자융합공학과 (공학사)

  • 2015년 3월~현재: 광운대학교 전파공학과 석사과정

  • [주 관심분야] 무선통신 시스템, 레이다 신호처리

전 현 무

jkiees-28-2-139-i2

  • 2012년 2월: 광운대학교 전자공학과 (공학사)

  • 2012년 3월~현재: 광운대학교 전파공학과 석박통합과정

  • [주 관심분야] 무선통신 시스템, 레이다 신호처리

양 동 혁

jkiees-28-2-139-i3

  • 1997년 2월: 공군사관학교 전자공학과 (공학사)

  • 2006년 6월: 미국 Wichita State University 전기전자공학과 (공학석사)

  • 2015년 3월~현재: 광운대학교 방위사업학과 박사과정

  • [주 관심분야] SAR ATR, SAR GMTI, Maximin algorithm, GPS 좌표 보정

양 훈 기

jkiees-28-2-139-i4

  • 1985년 2월: 연세대학교 전자공학과 (공학사)

  • 1987년 5월: 미국 SUNY at Buffalo University 전기컴퓨터공학과 (공학석사)

  • 1992년 4월: 미국 SUNY at Buffalo University 전기컴퓨터공학과 (공학박사)

  • 1993년 3월~현재: 광운대학교 전파공학과 교수

  • [주 관심분야] 무선통신 시스템, RFID, 레이다 신호처리


Notice

I don't want to open this window for a day.