#Create platform dictionary
platform = plat_dict()
#Create image plane dictionary
img_plane = imgTools.img_plane_dict(platform)
#Simulate phase history, if needed
##############################################################################
nsamples = platform['nsamples']
npulses = platform['npulses']
x = img_plane['u']; y = img_plane['v']
points = [[0,0,0],
[0,-100,0],
[200,0,0]]
amplitudes = [1,1,1]
phs = phsTools.simulate_phs(platform, points, amplitudes)
##############################################################################
#Apply RVP correction
phs_corr = phsTools.RVP_correct(phs, platform)
#Demodulate phase history with constant reference, if needed
phs_fixed = phsTools.phs_to_const_ref(phs_corr, platform, upchirp = 1)
#Apply algorithm of choice to phase history data
img_pf = imgTools.polar_format(phs_corr, platform, img_plane, taylor = 43)
#img_wk = imgTools.omega_k(phs_fixed, platform, taylor = 43, upsample = 2)
#img_bp = imgTools.backprojection(phs_corr, platform, img_plane, taylor = 0, upsample = 2)
#Output image
plt.imshow(np.abs(img_pf))
#Create platform dictionary platform = plat_dict() #Create image plane dictionary img_plane = imgTools.img_plane_dict(platform, aspect=1) #Simulate phase history, if needed ############################################################################## nsamples = platform['nsamples'] npulses = platform['npulses'] x = img_plane['u'] y = img_plane['v'] points = [[0, 0, 0], [0, -100, 0], [200, 0, 0]] amplitudes = [1, 1, 1] phs = phsTools.simulate_phs(platform, points, amplitudes) ############################################################################## #Apply RVP correction phs_corr = phsTools.RVP_correct(phs, platform) #Demodulate phase history with constant reference, if needed phs_fixed = phsTools.phs_to_const_ref(phs_corr, platform, upchirp=1) #Apply algorithm of choice to phase history data img_pf = imgTools.polar_format(phs_corr, platform, img_plane, taylor=17) #img_wk = imgTools.omega_k(phs_fixed, platform, taylor = 17, upsample = 2) #img_bp = imgTools.backprojection(phs_corr, platform, img_plane, taylor = 17, upsample = 2) #Output image imgTools.imshow(img_pf, dB_scale=[-25, 0])