In [1]:
%pylab inline
import warnings
import numpy as np
import matplotlib.pyplot as plt
import rayopt as ro
# ignore matplotlib and numpy warning each other
warnings.simplefilter("ignore", FutureWarning)
# ignore floating point exceptions
np.seterr(divide="ignore", invalid="ignore")
# by default only print 4 significant digits
np.set_printoptions(precision=4)
In [69]:
description = "Nikkor-Zoom-80-400mm-f4.5-5.6-US20200049962A1"
columns = "type roc distance diameter material"
# 80mm
text80 = """
O 0 0 100 AIR
S 182.816 5.0 70 1.90265/35.7
S 92.566 2.5 70 1.49782/82.6
S -707.416 10 70 AIR
S 83.365 0.1 70 1.437/95
S 1420.361 9.2 70 AIR
S 117.082 8.225 45 1.801/34.9
S -117.044 6.4 45 1.49782/82.6
S 61.183 2.2 45 AIR
S -265.081 5.81 37 1.816/46.6
S 30.785 2 37 1.80518/25.5
S 92.264 4.6 37 AIR
S -56.342 6.2 35 1.83481/42.7
S 158.965 2 35 AIR
S 112.252 27.059 35 1.59319/67.9
S -78.685 4.6 35 AIR
S 67.612 0.1 35 1.90366/31.3
S 35.499 1.8 35 1.59319/67.9
S -238.177 6.4 35 AIR
S -58.467 5.388 30 1.72916/54.6
S 38.999 1.6 30 1.90265/35.7
S 146.9 3.6 30 AIR
A 0 26.684 26.25 AIR
S 124.142 2 30 1.744/44.8
S -124.142 3.4 30 AIR
S 26.615 0.1 30 1.48749/70.3
S 0 6.8 30 1.95/29.4
S 26.437 2 30 1.51742/52.2
S 0 4.8 30 AIR
S 176.178 17.6 27 1.64769/33.7
S -19.703 6 27 1.603/65.4
S 0 1.6 27 AIR
S -22.131 11.27 30 1.83481/42.7
S -33.748 1.6 30 AIR
I 0 52.8 43.2 AIR
"""
# 400mm
text400 = """
O 0 0 100 AIR
S 182.816 5.0 70 1.90265/35.7
S 92.566 2.5 70 1.49782/82.6
S -707.416 10 70 AIR
S 83.365 0.1 70 1.437/95
S 1420.361 9.2 70 AIR
S 117.082 64.292 45 1.801/34.9
S -117.044 6.4 45 1.49782/82.6
S 61.183 2.2 45 AIR
S -265.081 5.81 37 1.816/46.6
S 30.785 2 37 1.80518/25.5
S 92.264 4.6 37 AIR
S -56.342 6.2 35 1.83481/42.7
S 158.965 2 35 AIR
S 112.252 3.056 35 1.59319/67.9
S -78.685 4.6 35 AIR
S 67.612 0.1 35 1.90366/31.3
S 35.499 1.8 35 1.59319/67.9
S -238.177 6.4 35 AIR
S -58.467 29.391 30 1.72916/54.6
S 38.999 1.6 30 1.90265/35.7
S 146.9 3.6 30 AIR
A 0 2.382 26.25 AIR
S 124.142 2 30 1.744/44.8
S -124.142 3.4 30 AIR
S 26.615 0.1 30 1.48749/70.3
S 0 6.8 30 1.95/29.4
S 26.437 2 30 1.51742/52.2
S 0 4.8 30 AIR
S 176.178 17.6 27 1.64769/33.7
S -19.703 6 27 1.603/65.4
S 0 1.6 27 AIR
S -22.131 11.27 30 1.83481/42.7
S -33.748 1.6 30 AIR
I 0 77.1 43.2 AIR
"""
text=text400
In [70]:
# wavelengths=[588e-9],
s = ro.system_from_text(text, columns.split(),
description=description)
s.object.angle = np.deg2rad(3.1)
s.fields = 0, .7, 1.
In [71]:
s.update()
#s.paraxial.resize()
print(s)
In [72]:
ro.Analysis(s)
Out[72]:
In [73]:
text=text80
In [74]:
s = ro.system_from_text(text, columns.split(),
description=description)
s.object.angle = np.deg2rad(14.8)
s.fields = 0, .7, 1.
In [75]:
s.update()
#s.paraxial.resize()
print(s)
In [76]:
ro.Analysis(s)
Out[76]:
In [ ]: