Following Hanser et. al.,
Noise will quickly corrupt the phase retrieval. The background has to be suppressed well.
TODO:
In [1]:
%load_ext autoreload
%autoreload 2
import os
from skimage.external import tifffile as tif
%pylab inline
from pyotf import HanserPSF
from pyotf.phaseretrieval import *
from pyotf.utils import *
from dphutils import scale, fft_pad
from dphplotting import mip, slice_plot, display_grid
from matplotlib.colors import LogNorm
In [2]:
test_data = tif.imread("../fixtures/psf_wl520nm_z300nm_x130nm_na0.85_n1.0.tif")
params = dict(
size = test_data.shape[-1],
zsize = test_data.shape[0],
na = 0.85,
res = 130,
zres = 300,
wl = 520,
ni = 1.0,
vec_corr="none",
condition="none"
)
test_data = tif.imread("../fixtures/BIL113_zstack_z300nm_current.tif") params = dict( size = test_data.shape[-1], zsize = test_data.shape[0], na = 1.1, res = 90.5, zres = 300, wl = 605, ni = 1.33, vec_corr="none", condition="none" )
In [3]:
junk = prep_data_for_PR(test_data, None, 1.5)
params["size"] = 256
result = retrieve_phase(junk, params, 25, pupil_tol=nan, mse_tol=nan, phase_only=False)
result.plot()
result.plot_convergence()
zd_result = result.fit_to_zernikes(120)
zd_result.plot_coefs()
matshow(zd_result.phase(slice(4, None, None)), cmap="seismic")
colorbar()
Out[3]:
In [4]:
junk = prep_data_for_PR(test_data, 512, 1.5)
params["size"] = 512
result = retrieve_phase(junk, params, 25, pupil_tol=nan, mse_tol=nan, phase_only=False)
result.plot()
result.plot_convergence()
zd_result2 = result.fit_to_zernikes(120)
zd_result2.plot_coefs()
matshow(zd_result2.phase(slice(4, None, None)), cmap="seismic")
Out[4]:
In [4]:
junk = prep_data_for_PR(test_data, 128, 1.5)
params["size"] = 128
result = retrieve_phase(junk, params, 25, pupil_tol=nan, mse_tol=nan, phase_only=False)
result.plot()
result.plot_convergence()
result.fit_to_zernikes(120)
result.zd_result.plot_coefs()
matshow(result.zd_result.phase(slice(4, None, None)), cmap="seismic")
colorbar()
result.zd_result.plot_named_coefs()
Out[4]:
In [28]:
for s in (128, 256, 512, 1024):
mip(result.generate_psf(size=s)[:, s // 2 - 32:s // 2 + 32, s // 2 - 32:s // 2 + 32])
In [11]:
plot(zd_result2.pcoefs - zd_result.pcoefs, label="2 - 1")
plot(zd_result3.pcoefs - zd_result.pcoefs, label="3 - 1")
plot(zd_result2.pcoefs - zd_result3.pcoefs, label="2 - 3")
legend()
Out[11]:
In [17]:
print(zd_result.zerns.shape)
print(zd_result2.zerns.shape)
print(zd_result3.zerns.shape)
Out[17]:
In [31]:
import tqdm
In [34]:
for i in tqdm.tqdm((0.5, 0.8, 1.0, 1.2, 1.5, 2.0)):
junk = prep_data_for_PR(test_data, None, i)
params["size"] = 256
result = retrieve_phase(junk, params, 25, pupil_tol=nan, mse_tol=nan, phase_only=False)
result.plot()
result.plot_convergence()
result.fit_to_zernikes(120)
result.zd_result.plot_coefs()
matshow(result.zd_result.phase(slice(4, None, None)), cmap="seismic")
colorbar()