python雷达成像（SAR）仿真：（二）距离压缩

回波信号承接上文：python雷达成像（SAR）仿真：（一）生成单目标数据

本节主要是对回波进行距离压缩

 

 

########### 利用回波数据成像 距离压缩###########
S_range = np.fft.fft(s_echo, NFFT_r)

# plt.pcolor(S_range.real, cmap='jet')
# plt.colorbar()
# pylab.title('S_range.real')
# plt.show()
#
# plt.pcolor(abs(S_range), cmap='jet')
# plt.colorbar()
# pylab.title('S_range.abs')
# plt.show()

t_ref = np.mat(np.arange(-Nr / 2, Nr / 2, 1) / Fr)  # 用来生成距离MF的距离时间轴
t_ref_mtx = np.dot(np.ones((Naz, 1)), t_ref)  # 矩阵形式
w_ref = np.mat(np.kaiser(Nr, 2.5))  # 距离向，构建Kaiser窗，此为列向量。
w_ref = np.dot(np.mat(np.ones((Naz, 1))), w_ref)  # 构成矩阵形式，每一行都相同的加窗。
s_ref = np.exp((1j * np.pi * Kr) * (np.power(t_ref_mtx, 2)))
# s_ref = np.pad(s_ref, (0, Nrg-int(Nr)), 'constant')  # 对复制脉冲，后端补零。
S_ref = np.fft.fft(s_ref, NFFT_r)
H_range = np.conj(S_ref)                  # 距离向匹配滤波器，零频在两端。
S_range_c = S_range*H_range           # 乘以匹配滤波器，零频在两端。
s_rc = np.fft.ifft(S_range_c)            # 完成距离压缩，回到二维时域。

plt.pcolor(S_range_c.real, cmap='jet')
plt.colorbar()
pylab.title('S_range_c.real')
plt.show()

plt.pcolor(abs(S_range_c), cmap='jet')
plt.colorbar()
pylab.title('S_range_c.abs')
plt.show()


plt.pcolor(s_rc.real, cmap='jet')
plt.colorbar()
pylab.title('s_rc.real')
plt.show()

plt.pcolor(abs(s_rc), cmap='jet')
plt.colorbar()
pylab.title('s_rc.abs')
plt.show()