In [5]:
%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 [14]:
# It appears that rayopt requires distance rather than
# thickness quoted in patent literature
# So we have to move the thickness value one row down
#
# The first surface represents the Object (O)
# The last surface represents the Image (I)
description = "Doublegauss-Nikkor-Haruo-58mmf1.4-JP2013-019993A"
columns = "type distance roc diameter material"
# Following uses Hikari glass types where possible
# The aspherical elements do not appear to have a good match
# probably because they have molded plastic?
# There is also a discrepancy in the element with
# index 1.883 - the catalog lists vd as 40.69
#
text = """
O 20 0 100 AIR
S 5.0 52.8577 47.2 1.74443/49.53
S 6.0 229.3475 47.2 AIR
S 0.1 40.3738 39.6 J-LASKH2
S 6.0 354.9744 39.6 J-FK5
S 1.5 42.4134 35.1 AIR
S 4.1038 290.8467 34.0 J-SF8
S 1.5 31.6359 33.0 AIR
A 6.0 0 28.4 AIR
S 6.0 -30.7873 30.0 J-SF10
S 1.7 35.1427 33.5 J-LASF08A
S 7.0 -131.1407 33.5 AIR
S 0.1 118.7661 33.8 J-LASF08A
S 6.0 -44.2318 33.8 J-LLF6
S 1.5 44.2683 33.8 1.74443/49.53
S 6.0 -77.2943 33.8 AIR
I 38.7 0 42.42 AIR
"""
text = """
O 20 0 100 AIR
S 5.0 52.8577 47.2 1.74443/49.53
S 6.0 229.3475 47.2 AIR
S 0.1 40.3738 39.6 1.755/52.34
S 6.0 354.9744 39.6 1.48749/70.31
S 1.5 42.4134 35.1 AIR
S 4.1038 290.8467 34.0 1.68893/31.16
S 1.5 31.6359 33.0 AIR
A 6.0 0 28.4 AIR
S 6.0 -30.7873 30.0 1.72825/28.46
S 1.7 35.1427 33.5 1.883/40.77
S 7.0 -131.1407 33.5 AIR
S 0.1 118.7661 33.8 1.883/40.66
S 6.0 -44.2318 33.8 1.53172/48.78
S 1.5 44.2683 33.8 1.74443/49.53
S 6.0 -77.2943 33.8 AIR
I 38.7 0 42.42 AIR
"""
In [15]:
s = ro.system_from_text(text, columns.split(),
description=description)
s.object.angle = np.deg2rad(20.0)
# κ A4 A6 A8 A10
# 1st surface 0.5721 1.10084E-07 6.21998E-10 -4.25694E-13 0.00000E+00
# κ A4 A6 A8 A10
# 15th 14.1597 8.65514E-06 4.15194E-09 1.25812E-11 1.22728E-14
## Aspheric data
## Note 1.0 has been subtracted from k because rayopt adds it
s[1].conic = -0.4279
s[1].aspherics = [0, 1.10084E-07, 6.21998E-10, -4.25694E-13]
s[15].conic = 13.1597
s[15].aspherics = [0, 8.65514E-06, 4.15194E-09, 1.25812E-11, 1.22728E-14]
s.fields = 0, .7, 1.
#s.wavelengths = [587.5618e-9]
In [16]:
s.update()
#s.paraxial.resize()
print(s)
In [17]:
ro.Analysis(s)
Out[17]:
In [ ]: