How to use the aotools.opticalpropagation.angularSpectrum function in aotools
To help you get started, we’ve selected a few aotools examples, based on popular ways it is used in public projects.
Secure your code as it's written. Use Snyk Code to scan source code in minutes - no build needed - and fix issues immediately.
AOtools / aotools / test / testOpticalprop.py View on Github 
def test_propagation_split_propagation():
screen = phasescreen.ft_phase_screen(0.16, 512, 4.2 / 512, 100, 0.01)
# Input E Field
E = numpy.exp(1j * screen)
prop1 = opticalpropagation.angularSpectrum(E, 500e-9, 4.2 / 512, 4.2 / 512, 5000.)
prop2 = opticalpropagation.angularSpectrum(prop1, 500e-9, 4.2 / 512, 4.2 / 512, 5000.)
prop3 = opticalpropagation.angularSpectrum(E, 500e-9, 4.2 / 512, 4.2 / 512, 10000.)
assert numpy.allclose(prop2, prop3)def test_propagation_round_trip():
screen = phasescreen.ft_phase_screen(0.16, 512, 4.2 / 512, 100, 0.01)
# Input E Field
E = numpy.exp(1j * screen)
prop1 = opticalpropagation.angularSpectrum(E, 500e-9, 4.2 / 512, 4.2 / 512, 10000.)
prop2 = opticalpropagation.angularSpectrum(prop1, 500e-9, 4.2 / 512, 4.2 / 512, -10000.)
assert numpy.allclose(E, prop2)def test_propagation_round_trip():
screen = phasescreen.ft_phase_screen(0.16, 512, 4.2 / 512, 100, 0.01)
# Input E Field
E = numpy.exp(1j * screen)
prop1 = opticalpropagation.angularSpectrum(E, 500e-9, 4.2 / 512, 4.2 / 512, 10000.)
prop2 = opticalpropagation.angularSpectrum(prop1, 500e-9, 4.2 / 512, 4.2 / 512, -10000.)
assert numpy.allclose(E, prop2)def test_propagation_split_propagation():
screen = phasescreen.ft_phase_screen(0.16, 512, 4.2 / 512, 100, 0.01)
# Input E Field
E = numpy.exp(1j * screen)
prop1 = opticalpropagation.angularSpectrum(E, 500e-9, 4.2 / 512, 4.2 / 512, 5000.)
prop2 = opticalpropagation.angularSpectrum(prop1, 500e-9, 4.2 / 512, 4.2 / 512, 5000.)
prop3 = opticalpropagation.angularSpectrum(E, 500e-9, 4.2 / 512, 4.2 / 512, 10000.)
assert numpy.allclose(prop2, prop3)def test_propagation_conserves_intensity():
screen = phasescreen.ft_phase_screen(0.16, 512, 4.2 / 512, 100, 0.01)
# Input E Field
E = numpy.exp(1j * screen)
Em = E * circle(150, 512)
sum1 = (abs(Em) ** 2).sum()
prop1 = opticalpropagation.angularSpectrum(Em, 500e-9, 4.2 / 512, 4.2 / 512, 10000.)
sum2 = (abs(prop1) ** 2).sum()
assert numpy.allclose(sum1, sum2)def test_propagation_split_propagation():
screen = phasescreen.ft_phase_screen(0.16, 512, 4.2 / 512, 100, 0.01)
# Input E Field
E = numpy.exp(1j * screen)
prop1 = opticalpropagation.angularSpectrum(E, 500e-9, 4.2 / 512, 4.2 / 512, 5000.)
prop2 = opticalpropagation.angularSpectrum(prop1, 500e-9, 4.2 / 512, 4.2 / 512, 5000.)
prop3 = opticalpropagation.angularSpectrum(E, 500e-9, 4.2 / 512, 4.2 / 512, 10000.)
assert numpy.allclose(prop2, prop3)aotools
A set of useful functions for Adaptive Optics in Python
Package Health Score
46 / 100Popular aotools functions
- aotools._version.NotThisMethod
- aotools.circle
- aotools.CovarianceMatrix
- aotools.fouriertransform
- aotools.fouriertransform.ft2
- aotools.fouriertransform.ift2
- aotools.functions
- aotools.functions.circle
- aotools.functions.pupil
- aotools.functions.pupil.circle
- aotools.functions.zernikeArray
- aotools.functions.zernikeRadialFunc
- aotools.image_processing
- aotools.image_processing.centreOfGravity
- aotools.image_processing.correlation_centroid
- aotools.image_processing.encircled_energy
- aotools.opticalpropagation.angularSpectrum
- aotools.turbulence
- aotools.turbulence.get_tps_time_axis
- aotools.turbulence.phasescreen.ft_phase_screen