###############################################################################
# Define top level directory containing *.mat file
# and choose polarization and starting azimuth
pol = 'HH'
directory = './data/AFRL/pass1'
start_az = 1
# Import phase history and create platform dictionary
[phs, platform] = phsRead.AFRL(directory, pol, start_az, n_az=4)
# Create image plane dictionary
img_plane = imgTools.img_plane_dict(platform, res_factor=1.0, upsample=True, aspect=1.0)
# full backprojection
start = time()
img_bp = imgTools.backprojection(phs, platform, img_plane, taylor=17, upsample=2)
bp_time = time() - start
# Fast-factorized backprojection without multi-processing
start = time()
img_FFBP = imgTools.FFBP(phs, platform, img_plane, taylor=17, factor_max=2)
fbp_time = time() - start
# Fast-factorized backprojection with multi-processing
start = time()
img_FFBP = imgTools.FFBPmp(phs, platform, img_plane, taylor=17, factor_max=2)
fbpmp_time = time() - start
# Output image
u = img_plane['u'];
v = img_plane['v']
from ritsar import imgTools
#Define top level directory containing *.mat file
#and choose polarization and starting azimuth
pol = 'HH'
directory = './data/AFRL/pass1'
start_az = 1
#Import phase history and create platform dictionary
[phs, platform] = phsRead.AFRL(directory, pol, start_az, n_az = 3)
#Create image plane dictionary
img_plane = imgTools.img_plane_dict(platform, res_factor = 1.5, upsample = True, aspect = 1.0)
#Apply algorithm of choice to phase history data
img_bp = imgTools.backprojection(phs, platform, img_plane, taylor = 43, upsample = 6)
#img_pf = imgTools.polar_format(phs, platform, img_plane, taylor = 43)
#Output image
plt.imshow(np.abs(img_bp)**0.5, cmap = cm.Greys_r)
plt.title('Backprojection')
#Autofocus image
print('autofocusing')
#img_af, af_ph = imgTools.autoFocus(img_bp, win = 0, win_params = [300,0.8])
img_af, af_ph = imgTools.autoFocus(img_bp, win = 0, win_params = [300,0.8])
#Output autofocused image
plt.figure()
plt.imshow(np.abs(img_af)**0.5, cmap = cm.Greys_r)
plt.title('Autofocused Polar Format')
directory = './data/AFRL/pass1'
start_az = 1
#Import phase history and create platform dictionary
[phs, platform] = phsRead.AFRL(directory, start_az, pol=pol, n_az=3)
#Create image plane dictionary
img_plane = imgTools.img_plane_dict(platform,
res_factor=1.4,
upsample=True,
aspect=1.0)
#Apply algorithm of choice to phase history data
img_bp = imgTools.backprojection(phs,
platform,
img_plane,
taylor=20,
upsample=6,
prnt=50)
#img_pf = imgTools.polar_format(phs, platform, img_plane, taylor = 20)
#Output image
imgTools.imshow(img_bp, dB_scale=[-30, 0])
plt.title('Backprojection')
#%%
#Define top level directory containing *.mat file
#and choose polarization and starting azimuth
directory = './data/AFRL/Wide_Angle_SAR'
start_az = 214
#Import phase history and create platform dictionary