def initDemoSystem(self):
     s = OpticalSystem()
     
     lc0 = s.addLocalCoordinateSystem(LocalCoordinates(name="object", decz=0.0), refname=s.rootcoordinatesystem.name)                
     lc1 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf1", decz=2.0), refname=lc0.name) # objectDist
     lc2 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf2", decz=3.0), refname=lc1.name)
     lc3 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf3", decz=5.0, tiltx=0.0*math.pi/180.0), refname=lc2.name)
     lc4 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf4", decz=3.0), refname=lc3.name)
     lc5 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf5", decz=3.0), refname=lc4.name)
     lc6 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf6", decz=2.0), refname=lc5.name)
     lc7 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf7", decz=3.0), refname=lc6.name)
     lc8 = s.addLocalCoordinateSystem(LocalCoordinates(name="image", decz=19.0), refname=lc7.name)
     
     
     objectsurf = Surface(lc0)
     surf1 = Surface(lc1, shape=surfShape.Conic(lc1, curv=1/-5.922))
     surf2 = Surface(lc2, shape=surfShape.Conic(lc2, curv=1/-3.160))
     surf3 = Surface(lc3, shape=surfShape.Conic(lc3, curv=1/15.884))
     surf4 = Surface(lc4, shape=surfShape.Conic(lc4, curv=1/-12.756))
     stopsurf = Surface(lc5)
     surf6 = Surface(lc6, shape=surfShape.Conic(lc6, curv=1/3.125))
     surf7 = Surface(lc7, shape=surfShape.Conic(lc7, curv=1/1.479))
     image = Surface(lc8)
     
     
     elem = OpticalElement(lc0, label="lenssystem")
     
     glass = material_isotropic.ConstantIndexGlass(lc0, n=1.7)
     glass2 = material_isotropic.ConstantIndexGlass(lc0, n=1.5)
     
     elem.addMaterial("glass", glass)
     elem.addMaterial("glass2", glass2)
     
     elem.addSurface("object", objectsurf, (None, None))
     elem.addSurface("surf1", surf1, (None, "glass"))
     elem.addSurface("surf2", surf2, ("glass", None))
     elem.addSurface("surf3", surf3, (None, "glass"))
     elem.addSurface("surf4", surf4, ("glass", None))
     elem.addSurface("stop", stopsurf, (None, None))
     elem.addSurface("surf6", surf6, (None, "glass2"))
     elem.addSurface("surf7", surf7, ("glass2", None))
     elem.addSurface("image", image, (None, None))
     
     s.addElement("lenssys", elem)
     
     return s
Beispiel #2
0
lc0 = s.addLocalCoordinateSystem(LocalCoordinates(name="object", decz=0.0), refname=s.rootcoordinatesystem.name)

#air = AnisotropicMaterial(lc0, myeps)  # tests for anisotropic mirror
air = ConstantIndexGlass(lc0, 1.0)
s.material_background = air

lc1 = s.addLocalCoordinateSystem(LocalCoordinates(name="m1", decz=50.0, tiltx=-math.pi/8), refname=lc0.name) # objectDist
lc2 = s.addLocalCoordinateSystem(LocalCoordinates(name="m2_stop", decz=-50.0, decy=-20, tiltx=math.pi/16), refname=lc1.name)
lc3 = s.addLocalCoordinateSystem(LocalCoordinates(name="m3", decz=50.0, decy=-30, tiltx=3*math.pi/32), refname=lc2.name)
lc4 = s.addLocalCoordinateSystem(LocalCoordinates(name="image1", decz=-50, decy=-15, tiltx=-math.pi/16), refname=lc3.name)
lc5 = s.addLocalCoordinateSystem(LocalCoordinates(name="oapara", decz=-100, decy=-35), refname=lc4.name)
lc5ap = s.addLocalCoordinateSystem(LocalCoordinates(name="oaparaap", decz=0, decy=35), refname=lc5.name)
lc6 = s.addLocalCoordinateSystem(LocalCoordinates(name="image2", decz=55, tiltx=1*math.pi/32), refname=lc5.name)
lc7 = s.addLocalCoordinateSystem(LocalCoordinates(name="image3", decz=5), refname=lc6.name)

objectsurf = Surface(lc0)
m1surf = Surface(lc1, shape=surfShape.Conic(lc1, curv=-0.01), apert=CircularAperture(lc1, 20.))
m2surf = Surface(lc2, shape=surfShape.Conic(lc2, curv=0.01), apert=CircularAperture(lc2, 12.7))
m3surf = Surface(lc3, shape=surfShape.Conic(lc3, curv=-0.006), apert=CircularAperture(lc3, 12.7))
image1 = Surface(lc4)
oapara = Surface(lc3, shape=surfShape.Conic(lc5, curv=0.01, cc=-1.), apert=CircularAperture(lc5ap, 30.0))
image2 = Surface(lc6, apert=CircularAperture(lc6, 20.0))
image3 = Surface(lc7, apert=CircularAperture(lc7, 20.0))


elem = OpticalElement(lc0, label="TMA")

elem.addMaterial("air", air)

elem.addSurface("object", objectsurf, (None, None))
elem.addSurface("m1", m1surf, (None, None))
Beispiel #3
0
lc1 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf1",
                                                  decz=-dropletradius),
                                 refname=lccomprism.name)  # objectDist
lc2 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf2",
                                                  decz=dropletradius + 0.0001),
                                 refname=lccomprism.name)
lc3 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf3",
                                                  decz=-dropletradius,
                                                  tiltx=math.pi),
                                 refname=lccomprism.name)
lc4 = s.addLocalCoordinateSystem(LocalCoordinates(name="image",
                                                  decz=-2. * dropletradius),
                                 refname=lccomprism.name)

stopsurf = Surface(lc0, apert=CircularAperture(lc0, dropletradius))
frontsurf = Surface(lc1,
                    shape=surfShape.Conic(lc1,
                                          curv=1. / (10 * dropletradius),
                                          cc=0.1),
                    apert=CircularAperture(lc1, dropletradius))
rearsurf = Surface(lc2,
                   shape=surfShape.Conic(lc2, curv=-1. / (20 * dropletradius)),
                   apert=CircularAperture(lc2, dropletradius))
frontsurf2 = Surface(lc3,
                     shape=surfShape.Conic(lc3,
                                           curv=-1. / (30 * dropletradius),
                                           cc=0.01),
                     apert=CircularAperture(lc3, dropletradius))

image = Surface(lc4, apert=CircularAperture(lc4, 7. * dropletradius))
import logging
logging.basicConfig(level=logging.DEBUG)

wavelength = 0.5876e-3

# definition of optical system
s = OpticalSystem(name='os') 

lc0 = s.addLocalCoordinateSystem(LocalCoordinates(name="stop", decz=0.0), refname=s.rootcoordinatesystem.name)
lc1 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf1", decz=-1.048), refname=lc0.name) # objectDist
lc2 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf2", decz=4.0), refname=lc1.name)
lc3 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf3", decz=2.5), refname=lc2.name)
lc4 = s.addLocalCoordinateSystem(LocalCoordinates(name="image", decz=97.2), refname=lc3.name)


stopsurf = Surface(lc0, name="stopsurf")
frontsurf = Surface(lc1, name="frontsurf", shape=surfShape.Conic(lc1, curv=1./62.8, name='conic1'), apert=CircularAperture(lc1, 12.7))
cementsurf = Surface(lc2, name="cementsurf", shape=surfShape.Conic(lc2, curv=-1./45.7, name='conic2'), apert=CircularAperture(lc2, 12.7))
rearsurf = Surface(lc3, name="rearsurf", shape=surfShape.Conic(lc3, curv=-1./128.2, name='conic3'), apert=CircularAperture(lc3, 12.7))
image = Surface(lc4, name="imagesurf")


elem = OpticalElement(lc0, name="thorlabs_AC_254-100-A")

rnd_data1 = np.random.random((3, 3)) #np.eye(3)
rnd_data2 = np.random.random((3, 3))#np.zeros((3, 3))#
rnd_data3 = np.random.random((3, 3)) #np.eye(3)
rnd_data4 = np.random.random((3, 3))#np.zeros((3, 3))#

# isotropic tests
Beispiel #5
0
lccomprism = s.addLocalCoordinateSystem(LocalCoordinates(name="prismcenter",
                                                         decz=50.0),
                                        refname=lc0.name)

lc1 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf1",
                                                  decz=-10.0,
                                                  tiltx=30. * deg),
                                 refname=lccomprism.name)  # objectDist
lc2 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf2",
                                                  decz=10.0,
                                                  tiltx=-30. * deg),
                                 refname=lccomprism.name)
lc3 = s.addLocalCoordinateSystem(LocalCoordinates(name="image", decz=50.0),
                                 refname=lccomprism.name)

stopsurf = Surface(lc0)
frontsurf = Surface(lc1,
                    shape=surfShape.Conic(lc1, curv=0),
                    apert=CircularAperture(lc1, 20.0))
rearsurf = Surface(lc2,
                   shape=surfShape.Conic(lc2, curv=0),
                   apert=CircularAperture(lc2, 20.0))
image = Surface(lc3)

elem = OpticalElement(lc0, label="prism")

glass = ModelGlass(lc1)

elem.addMaterial("glass", glass)

elem.addSurface("stop", stopsurf, (None, None))
Beispiel #6
0
wavelength = 0.5876e-3

# definition of optical system
s = OpticalSystem()

lc0 = s.addLocalCoordinateSystem(LocalCoordinates(name="stop", decz=0.0),
                                 refname=s.rootcoordinatesystem.name)
lc1 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf1", decz=5.0),
                                 refname=lc0.name)  # objectDist
lc2 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf2", decz=20.0),
                                 refname=lc1.name)
lc3 = s.addLocalCoordinateSystem(LocalCoordinates(name="image", decz=100.0),
                                 refname=lc2.name)

stopsurf = Surface(lc0)
frontsurf = Surface(lc1,
                    shape=surfShape.Conic(lc1),
                    apert=CircularAperture(lc1, 12.7))
backsurf = Surface(lc2,
                   shape=surfShape.Asphere(lc2,
                                           curv=-1. / 50.0,
                                           cc=-1.,
                                           acoeffs=[0.0, 0.0, 0.0]),
                   apert=CircularAperture(lc2, 12.7))
image = Surface(lc3)

elem = OpticalElement(lc0, label="asphereelement")

bk7 = ConstantIndexGlass(lc1, n=1.5168)
Beispiel #7
0
lcD4prime = s.addLocalCoordinateSystem(LocalCoordinates(
    name="D4prime",
    decy=-23.067,
    decz=36.667,
    tiltx=-si * 50.668 * degree,
    tiltThenDecenter=False),
                                       refname=lc0.name)
lcimage = s.addLocalCoordinateSystem(LocalCoordinates(name="image",
                                                      decy=-30.892,
                                                      decz=43.083,
                                                      tiltx=-si * 50.668 *
                                                      degree,
                                                      tiltThenDecenter=False),
                                     refname=lc0.name)

objsurf = Surface(lc0)
D1surf = Surface(lcD1)
#S1surf = Surface(lcS1, shape=Conic(lcS1, curv=si*1./108.187, cc=0), apert=CircularAperture(lcS1, 100.0))
S1surf = Surface(lcS1,
                 shape=Biconic(lcS1,
                               curvy=si * 1. / 108.187,
                               curvx=si * 1. / 73.105,
                               coefficients=[(0., 0.), (-si * 5.542e-7, -0.08),
                                             (-si * 8.176e-11, -1.379)]),
                 apert=CircularAperture(lcS1, 40.0))
D1Psurf = Surface(lcD1prime)
D2surf = Surface(lcD2)
#S2surf = Surface(lcS2, shape=Conic(lcS2, curv=si*1./69.871, cc=-0.1368), apert=CircularAperture(lcS2, 60.0))
S2surf = Surface(lcS2,
                 shape=Biconic(lcS2,
                               curvy=si * 1. / 69.871,
Beispiel #8
0
lc3 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf3",
                                                  decz=5.0,
                                                  tiltx=2.5 * math.pi / 180.0),
                                 refname=lc2.name)
lc4 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf4", decz=3.0),
                                 refname=lc3.name)
lc5 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf5", decz=3.0),
                                 refname=lc4.name)
lc6 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf6", decz=2.0),
                                 refname=lc5.name)
lc7 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf7", decz=3.0),
                                 refname=lc6.name)
lc8 = s.addLocalCoordinateSystem(LocalCoordinates(name="image", decz=19.0),
                                 refname=lc7.name)

objectsurf = Surface(lc0)
surf1 = Surface(lc1, shape=surfShape.Conic(lc1, curv=1 / -5.922))
surf2 = Surface(lc2, shape=surfShape.Conic(lc2, curv=1 / -3.160))
surf3 = Surface(lc3, shape=surfShape.Conic(lc3, curv=1 / 15.884))
surf4 = Surface(lc4, shape=surfShape.Conic(lc4, curv=1 / -12.756))
stopsurf = Surface(lc5)
surf6 = Surface(lc6, shape=surfShape.Conic(lc6, curv=1 / 3.125))
surf7 = Surface(lc7, shape=surfShape.Conic(lc7, curv=0.1 * 1 / 1.479))
image = Surface(lc8)

elem = OpticalElement(lc0, name="lenssystem")

glass = ConstantIndexGlass(lc0, n=1.7)
glass2 = ConstantIndexGlass(lc0, n=1.5)

elem.addMaterial("glass", glass)
Beispiel #9
0
# definition of optical system
s = OpticalSystem() 

dropletradius = 0.1

lc0 = s.addLocalCoordinateSystem(LocalCoordinates(name="stop", decz=0.0), refname=s.rootcoordinatesystem.name)
lccomprism = s.addLocalCoordinateSystem(LocalCoordinates(name="dropletcenter", decz=2.*dropletradius), refname=lc0.name)

lc1 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf1", decz=-dropletradius), refname=lccomprism.name) # objectDist
lc2 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf2", decz=dropletradius), refname=lccomprism.name)
lc3 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf3", decz=0), refname=lccomprism.name)
lc4 = s.addLocalCoordinateSystem(LocalCoordinates(name="image", decz=-2.*dropletradius), refname=lccomprism.name)


stopsurf = Surface(lc0, apert=CircularAperture(lc0, 7*dropletradius))
frontsurf = Surface(lc1, shape=Asphere(lc1, curv=1./dropletradius), apert=CircularAperture(lc1, dropletradius))
rearsurf = Surface(lc2, shape=Asphere(lc2, curv=-1./dropletradius), apert=CircularAperture(lc2, dropletradius))
midsurf = Surface(lc3, shape=Asphere(lc3, curv=0), apert=CircularAperture(lc3, dropletradius))

image = Surface(lc4, apert=CircularAperture(lc4, 7.*dropletradius))


elem = OpticalElement(lc0, name="droplet")



try:
    database_basepath = "refractiveindex.info-database/database"     
    shelf = "3d"
    book  = "liquids"
Beispiel #10
0
s = OpticalSystem()

lc0 = s.addLocalCoordinateSystem(LocalCoordinates(name="obj", decz=0.0),
                                 refname=s.rootcoordinatesystem.name)
lc1 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf1",
                                                  decz=10.0,
                                                  tiltx=5. * math.pi / 180.0),
                                 refname=lc0.name)  # objectDist
lc2 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf2",
                                                  decz=20.0,
                                                  tiltx=10. * math.pi / 180.0),
                                 refname=lc1.name)
lc3 = s.addLocalCoordinateSystem(LocalCoordinates(name="image", decz=10.0),
                                 refname=lc2.name)

stopsurf = Surface(lc0)
surf1 = Surface(lc1,
                shape=surfShape.Conic(lc1, curv=1. / 24.0),
                apert=CircularAperture(lc1, 5.0))
surf2 = Surface(lc2,
                shape=surfShape.Conic(lc2, curv=-1. / 24.0),
                apert=CircularAperture(lc2, 5.0))
image = Surface(lc3)

elem = OpticalElement(lc0, label="grinelement")

grin_strength = 0.5


def nfunc(x):
    return grin_strength * np.exp(
# definition of optical system
s = OpticalSystem()

lc0 = s.addLocalCoordinateSystem(LocalCoordinates(name="stop", decz=0.0),
                                 refname=s.rootcoordinatesystem.name)
lc1 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf1", decz=-1.048),
                                 refname=lc0.name)  # objectDist
lc2 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf2", decz=4.0),
                                 refname=lc1.name)
lc3 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf3", decz=2.5),
                                 refname=lc2.name)
lc4 = s.addLocalCoordinateSystem(LocalCoordinates(name="image", decz=97.2),
                                 refname=lc3.name)

stopsurf = Surface(lc0)
frontsurf = Surface(lc1,
                    shape=surfShape.Conic(lc1, curv=1. / 62.8),
                    apert=CircularAperture(lc1, 12.7))
cementsurf = Surface(lc2,
                     shape=surfShape.Conic(lc2, curv=-1. / 45.7),
                     apert=CircularAperture(lc2, 12.7))
rearsurf = Surface(lc3,
                   shape=surfShape.Conic(lc3, curv=-1. / 128.2),
                   apert=CircularAperture(lc3, 12.7))
image = Surface(lc4)

elem = OpticalElement(lc0, label="thorlabs_AC_254-100-A")

rnd_data1 = np.random.random((3, 3))  #np.eye(3)
rnd_data2 = np.random.random((3, 3))  #np.zeros((3, 3))#