#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))
#Add include directories to default path list
from sys import path
path.append('../')
#Include SARIT toolset
from ritsar import phsRead
from ritsar import phsTools
from ritsar import imgTools
#Define directory containing *.au2 and *.phs files
directory = './data/DIRSIG/'
#Import phase history and create platform dictionary
[phs, platform] = phsRead.DIRSIG(directory)
#Correct for reisdual video phase
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)
#Import image plane dictionary from './parameters/img_plane'
img_plane = imgTools.img_plane_dict(platform, res_factor = 1.0, aspect = 1.0)
#Apply polar format algorithm to phase history data
img_wk = imgTools.omega_k(phs_fixed, platform, taylor = 13, upsample = 2)
#img_bp = imgTools.backprojection(phs_corr, platform, img_plane, taylor = 13, upsample = 6)
#img_pf = imgTools.polar_format(phs_corr, platform, img_plane, taylor = 13)
#Output image
imgTools.imshow(img_wk, [-20,0])
