Mirror and lens in an imperfect system #165
-
|
I'm attempting to build a complex system to include some Monte Carlo tolerancing analysis. The bits that are getting me stuck are mostly the effect of a mirror and what happens to the elements after the mirror. I've simplified my optical chain to those parts. I would very much like to solve this issue. I also need to put in my own array of starting rays so I've included some of the code for that as well. `opm = OpticalModel() osp['pupil'] = PupilSpec(osp, key=['object', 'pupil'], value=40) BFL = 10 opm.system_spec.dimensions = 'mm' sm.gaps[0].thi=1e-3 sm.add_surface([1e6, BFL, 'air']) opm.add_dummy_plane(idx=1) sm.add_surface([1e6, BFL/2, 'air']) sm.add_surface([1e4, 10, 'N-BK7','Schott', 10]) opm.update_model() fig = plt.figure(FigureClass=InteractiveLayout, opt_model=opm, is_dark=False, do_draw_rays=False, do_draw_beams=False, do_draw_edge_rays=False) def plot_ray(ray): for m, pos in enumerate([-1, 0, 1]): This is the result. Clearly reflection off the mirror is not correct and the lens just disappears.
|
Beta Was this translation helpful? Give feedback.
Replies: 7 comments 31 replies
-
|
Ok - almost got something working. It appears that the ray tracing may be accurately represented but the lens after the mirror is not drawn correctly. This isn't a problem for getting results, but makes it harder to convince others that you've done the right job! opm = OpticalModel() osp['pupil'] = PupilSpec(osp, key=['object', 'pupil'], value=40) BFL = 20 opm.system_spec.dimensions = 'mm' sm.gaps[0].thi = BFL/2 opm.add_dummy_plane(idx=1) sm.gaps[-1].thi = BFL/2 # need to move the mirror by ~mirror_radius/sqrt(2) sign = 1 # change to -1 when mirror in play sm.add_surface([sign1e4, sign10, 'N-BK7','Schott', 10]) opm.update_model() fig = plt.figure(FigureClass=InteractiveLayout, opt_model=opm, is_dark=False, do_draw_rays=False, do_draw_beams=False, do_draw_edge_rays=False) def plot_ray(ray): for m, pos in enumerate([-1, 0, 1]): If the mirror is not inserted this is the lens picture: |
Beta Was this translation helpful? Give feedback.
-
|
I'm looking thru the code and there are a lot of places where the YZ plane is baked into the code. The profile classes always sample in the local YZ plane, for instance. If you have spheres or rotationally symmetric surfaces, this should be OK. However, the APIs used for getting this info would need to be modified to include which slice to use for the general surface profile. There is little useful graphics data generated for a XY plane view, i.e. looking along the z axis. The circular edges you'd expect looking at a lens face-on are not currently generated. I'm not sure what I'm looking at, so I'm not sure what's missing. If you could post the model you're working with, that would help. Regards, |
Beta Was this translation helpful? Give feedback.
-
|
Ok - here's the whole code: `%matplotlib widget initializationfrom rayoptics.environment import * Input lens:n1 = 1.444 # Fused Silica, Corning 7980 OF C to C gap between L1 and L2d1 = 23.0 Center lens:n2 = 1.444 # Fused Silica, Corning 7980 OF C to C gap between L2 and L3d2 = 27.0 Output lens:n3 = 1.444 # Fused Silica, Corning 7980 OF BFL = 40 opm = OpticalModel() FFL = 10 L1_tilt = 5 mirror_sign = 1 opm.system_spec.dimensions = 'mm' sm.gaps[0].thi = BFL/2 sm.gaps[-1].thi = BFL/2 # need to move the mirror by ~mirror_radius/sqrt(2) sm.add_surface([mirror_sign1e6, mirror_signFFL/2, 'air']) sm.add_surface([mirror_signR11, mirror_signt1, g1, m1, D1/1]) sm.add_surface([mirror_signR21, mirror_signt2, g2, m2, D2/1]) sm.add_surface([mirror_signR31, mirror_signt3, g3, m3, D3/1]) sm.gaps[-1].thi = mirror_sign2FFL osp['pupil'] = PupilSpec(osp, key=['object', 'pupil'], value=40) opm.update_model() fig = plt.figure(FigureClass=InteractiveLayout, opt_model=opm, is_dark=False, do_draw_rays=False, do_draw_beams=False, do_draw_edge_rays=False) #, clip_rays=True) def plot_ray(ray): for m, pos in enumerate([-1, 0, 1]): I also attempted to replace my pupil spec: with the code below, but had errors so I left it as is for now. ' osp['pupil'] = PupilSpec(osp, key=['object', 'epd'], value=2.) |
Beta Was this translation helpful? Give feedback.
-
|
The next optic is a mirror that will reflect into the x-z plane. I might forge ahead with the model, but it's tricky to see what's happening. Another solution is to build three models each rotated by 90 degrees. Updating those could be error prone though. |
Beta Was this translation helpful? Give feedback.
-
|
Do you have a way to initiate rays traveling along x or y? |
Beta Was this translation helpful? Give feedback.
-
|
Hello @justintorgerson, initializationfrom rayoptics.environment import *
from opticalglass.glassfactory import create_glass
from opticalglass.rindexinfo import summary_plots
SiO2_url = 'https://refractiveindex.info/database/data-nk/main/SiO2/Malitson.yml'
SiO2 = create_glass(SiO2_url, "rindexinfo")# Input lens:
n1 = 1.444 # Fused Silica, Corning 7980 OF
g1 = SiO2
m1 = ''
R11 = 190.0 # (+) for CX, (-) for CC
t1 = 6.0
R12 = 50.0 # (+) for CX, (-) for CC
D1 = 12.0
#C to C gap between L1 and L2
d1 = 23.0
#Center lens:
n2 = 1.444 # Fused Silica, Corning 7980 OF
g2 = SiO2
m2 = ''
R21 = -25.0 # (+) for CX, (-) for CC
t2 = 10.0
R22 = -45.0 # (+) for CX, (-) for CC
D2 = 12.0
#C to C gap between L2 and L3
d2 = 27.0
#Output lens:
n3 = 1.444 # Fused Silica, Corning 7980 OF
g3 = SiO2
m3 = ''
R31 = -200.0 # (+) for CX, (-) for CC
t3 = 6.0
R32 = 55.0 # (+) for CX, (-) for CC
D3 = 12.0
BFL = 40set up comments
opm = OpticalModel()
sm = opm['seq_model']
osp = opm['optical_spec']
pm = opm['parax_model']
em = opm['ele_model']
pt = opm['part_tree']
ar = opm['analysis_results']
FFL = 10
L1_tilt = 5
L1_shift = 1
L2_tilt = -5
L2_shift = -1
L3_tilt = 5
L3_shift = 1
mirror_sign = 1
opm.system_spec.dimensions = 'mm'
opm.radius_mode = True
sm.do_apertures = False
sm.gaps[0].thi = BFL/2
opm.add_dummy_plane()
sm.set_stop() # target for the rays so ray origin and direction is fixed
sm.gaps[-1].thi = BFL/2 # need to move the mirror by ~mirror_radius/sqrt(2)
sm.add_surface([.0, -BFL/2, 'refl', 10]) # need to move away from the mirror by ~mirror_radius/sqrt(2)
mirror_sign = -1*mirror_sign # change to -1 for odd bounces off mirror
mirror_surface = sm.cur_surface
sm.ifcs[mirror_surface].decenter = srf.DecenterData('bend', alpha=-45.) # rotate the mirror about x-axis (into page)
sm.add_surface([0, mirror_sign*FFL/2, 'air'])
sm.add_surface([mirror_sign*R11, mirror_sign*t1, g1, m1, D1/1])
sm.ifcs[sm.cur_surface].decenter = srf.DecenterData('decenter')
sm.ifcs[sm.cur_surface].decenter.euler[0] = L1_tilt
sm.ifcs[sm.cur_surface].decenter.dec[1] = L1_shift
sm.add_surface([-mirror_sign*R12, 0, 'air', D1/1])
sm.ifcs[sm.cur_surface].decenter = srf.DecenterData('reverse')
sm.ifcs[sm.cur_surface].decenter.euler[0] = L1_tilt
sm.ifcs[sm.cur_surface].decenter.dec[1] = L1_shift + mirror_sign*t1*math.sin(np.deg2rad(L1_tilt))
#opm.add_dummy_plane()
sm.gaps[-1].thi = mirror_sign*d1
sm.add_surface([mirror_sign*R21, mirror_sign*t2, g2, m2, D2/1])
sm.ifcs[sm.cur_surface].decenter = srf.DecenterData('decenter')
sm.ifcs[sm.cur_surface].decenter.euler[0] = L2_tilt
sm.ifcs[sm.cur_surface].decenter.dec[1] = L2_shift
sm.add_surface([-mirror_sign*R22, 0, 'air', D2/1])
sm.ifcs[sm.cur_surface].decenter = srf.DecenterData('reverse')
sm.ifcs[sm.cur_surface].decenter.euler[0] = L2_tilt
sm.ifcs[sm.cur_surface].decenter.dec[1] = L2_shift + mirror_sign*t2*math.sin(np.deg2rad(L2_tilt))
sm.gaps[-1].thi = mirror_sign*d2
sm.add_surface([mirror_sign*R31, mirror_sign*t3, g3, m3, D3/1])
sm.ifcs[sm.cur_surface].decenter = srf.DecenterData('decenter')
sm.ifcs[sm.cur_surface].decenter.euler[0] = L3_tilt
sm.ifcs[sm.cur_surface].decenter.dec[1] = L3_shift
sm.add_surface([-mirror_sign*R32, 0, 'air', D3/1])
sm.ifcs[sm.cur_surface].decenter = srf.DecenterData('reverse')
sm.ifcs[sm.cur_surface].decenter.euler[0] = L3_tilt
sm.ifcs[sm.cur_surface].decenter.dec[1] = L3_shift + mirror_sign*t3*math.sin(np.deg2rad(L3_tilt))
sm.gaps[-1].thi = mirror_sign*2*FFL
osp['pupil'] = PupilSpec(osp, key=['object', 'epd'], value=2.)
osp['fov'] = FieldSpec(osp, key=('object', 'angle'), value=1.,
flds=[0., 1., -1.],
is_relative=True)
osp['wvls'] = WvlSpec([(940.0, 1.0)], ref_wl=0)
sm.gaps[0].thi = 1.0e+10
opm.update_model()sm.list_model()Use list_sg to get the decenter data with the rest of the sequential model. sm.list_sg()Use list_gbl_tfrms to get the rotation matrices and translation ventors for each interface with respect to the 1st interface. sm.list_gbl_tfrms()You can list the global data with respect to any interface by using the sm.list_tfrms(sm.compute_global_coords(glo=4));fig = plt.figure(FigureClass=InteractiveLayout, opt_model=opm, is_dark=False) #, clip_rays=True)
fig.plot()
Ray tracingIt's not clear to me what type of system you are trying to model. What I set up was a system with collimated input over a 2mm stop diameter. I trace collimated rays at 0, 1 and -1 degrees thru the stop surface. Another approach is to have a finite object and a f/# or numerical aperture (NA) specifiation in object space. This might be typical of a microscope system, for example. Here is the optical specification data listobj(osp)You can get the ray trace output listed in both local and global coordinates. optical_axis_ray = trace_ray(opm, [0, 0], osp['fov'].fields[0], osp['wvls'].central_wvl)
list_ray(optical_axis_ray)The list_ray function takes a list_ray(optical_axis_ray, tfrms=sm.gbl_tfrms) |
Beta Was this translation helpful? Give feedback.
-
|
Thank you. I will include your suggestions in my code now. I have another question too. I would like the 45 degree fold mirror to be able to adjust over a small range so that I can assess some tolerancing issues. However, it moves the optic axis when I change it's angle. How do I change the mirror angle but not change the position and orientation of the lenses above it? |
Beta Was this translation helpful? Give feedback.
-
|
I have another question now too. I'm wanting to see the progression of 'sm.gbl_tfrms' as I add elements to the system but if I use 'opm.update_model()' to update the 'sm.gbl_tfrms' after each element, I will either get an error saying that "AttributeError: 'NoneType' object has no attribute 'name'" or some of the elements won't show up in the figure when I plot it. When they don't show up in the figure, it appears that they are still in the system because the rays are still refracting or reflecting. |
Beta Was this translation helpful? Give feedback.
-
|
Hello @justintorgerson, initializationopm = OpticalModel()
sm = opm['seq_model']
osp = opm['optical_spec']
pm = opm['parax_model']
em = opm['ele_model']
pt = opm['part_tree']
ar = opm['analysis_results']
FFL = 10
bend_angle = -45
tilt_error = 0
L1_tilt = 5
L1_shift = 1
L2_tilt = -5
L2_shift = -1
L3_tilt = 5
L3_shift = 1
mirror_sign = 1
opm.system_spec.dimensions = 'mm'
opm.radius_mode = True
sm.do_apertures = False
sm.gaps[0].thi = 1.0e+10
opm.add_dummy_plane(t=BFL/2) # need to move the mirror by ~mirror_radius/sqrt(2)
sm.set_stop() # target for the rays so ray origin and direction is fixedSet up to perturb fold mirrorThere are several ways that you could use decenters to achieve this. IMO, the best way to model this is to separate out the fold mirror perturbation from the nominal fold mirror bend angle. For this we don't use the I used the element form of entry, e.g. opm.add_dummy_plane()
sm.ifcs[sm.cur_surface].decenter = srf.DecenterData('decenter')
sm.ifcs[sm.cur_surface].decenter.euler[0] = bend_angle
opm.add_mirror(sd=10)
mirror_sign = -1*mirror_sign # change to -1 for odd bounces off mirror
mirror_surface = sm.cur_surface
sm.ifcs[sm.cur_surface].decenter = srf.DecenterData('dec and return')
sm.ifcs[sm.cur_surface].decenter.euler[0] = tilt_error
opm.add_dummy_plane(t=-BFL/2)
sm.ifcs[sm.cur_surface].decenter = srf.DecenterData('decenter')
sm.ifcs[sm.cur_surface].decenter.euler[0] = bend_angleadd following lensessm.add_surface([0, mirror_sign*FFL/2, 'air'])
sm.add_surface([mirror_sign*R11, mirror_sign*t1, g1, m1, D1/1])
sm.ifcs[sm.cur_surface].decenter = srf.DecenterData('decenter')
sm.ifcs[sm.cur_surface].decenter.euler[0] = L1_tilt
sm.ifcs[sm.cur_surface].decenter.dec[1] = L1_shift
sm.add_surface([-mirror_sign*R12, 0, 'air', D1/1])
sm.ifcs[sm.cur_surface].decenter = srf.DecenterData('reverse')
sm.ifcs[sm.cur_surface].decenter.euler[0] = L1_tilt
sm.ifcs[sm.cur_surface].decenter.dec[1] = L1_shift + mirror_sign*t1*math.sin(np.deg2rad(L1_tilt))
sm.gaps[sm.cur_surface].thi = mirror_sign*d1
sm.add_surface([mirror_sign*R21, mirror_sign*t2, g2, m2, D2/1])
sm.ifcs[sm.cur_surface].decenter = srf.DecenterData('decenter')
sm.ifcs[sm.cur_surface].decenter.euler[0] = L2_tilt
sm.ifcs[sm.cur_surface].decenter.dec[1] = L2_shift
sm.add_surface([-mirror_sign*R22, 0, 'air', D2/1])
sm.ifcs[sm.cur_surface].decenter = srf.DecenterData('reverse')
sm.ifcs[sm.cur_surface].decenter.euler[0] = L2_tilt
sm.ifcs[sm.cur_surface].decenter.dec[1] = L2_shift + mirror_sign*t2*math.sin(np.deg2rad(L2_tilt))
sm.gaps[sm.cur_surface].thi = mirror_sign*d2
sm.add_surface([mirror_sign*R31, mirror_sign*t3, g3, m3, D3/1])
sm.ifcs[sm.cur_surface].decenter = srf.DecenterData('decenter')
sm.ifcs[sm.cur_surface].decenter.euler[0] = L3_tilt
sm.ifcs[sm.cur_surface].decenter.dec[1] = L3_shift
sm.add_surface([-mirror_sign*R32, 0, 'air', D3/1])
sm.ifcs[sm.cur_surface].decenter = srf.DecenterData('reverse')
sm.ifcs[sm.cur_surface].decenter.euler[0] = L3_tilt
sm.ifcs[sm.cur_surface].decenter.dec[1] = L3_shift + mirror_sign*t3*math.sin(np.deg2rad(L3_tilt))
sm.gaps[sm.cur_surface].thi = mirror_sign*2*FFL
osp['pupil'] = PupilSpec(osp, key=['object', 'epd'], value=2.)
osp['fov'] = FieldSpec(osp, key=('object', 'angle'), value=1.,
flds=[0., 1., -1.],
is_relative=True)
osp['wvls'] = WvlSpec([(940.0, 1.0)], ref_wl=0)
sm.gaps[0].thi = 1.0e+10
opm.update_model()sm.list_sg()sm.list_gbl_tfrms()fig1 = plt.figure(FigureClass=InteractiveLayout, opt_model=opm, is_dark=False).plot()
Add a tilt error to the fold mirrorNote that the positions of the lens elements is unchanged but the rays follow a different path. tilt_error = 2
sm.ifcs[mirror_surface].decenter.euler[0] = tilt_error
opm.update_model()fig2 = plt.figure(FigureClass=InteractiveLayout, opt_model=opm, is_dark=False).plot()
sm.list_sg()sm.list_gbl_tfrms() |
Beta Was this translation helpful? Give feedback.
-
|
I've been charging ahead with some other issues and decided to get back to your previous comments and suggestions to be on the same page. However, I'm executing this code coped from above from 4 days ago: And am getting this error: UnboundLocalError: local variable 'slp0' referenced before assignment |
Beta Was this translation helpful? Give feedback.
-
|
The whole error chain looks like this: |
Beta Was this translation helpful? Give feedback.
-
|
Turn out the error is from the pupil type 'epd'. If I use 'pupil' I get the rays to plot, but no optics are shown. Here is my code: And this is the plot: Also getting some artifacts. I haven't dug into them yet. Do they add to the optical path? |
Beta Was this translation helpful? Give feedback.
-
|
I'm having trouble duplicating your issues. I've added a few print stmts to list versions. %matplotlib inline
# initialization
from rayoptics.environment import *import sys
print(f"python {sys.version}")
print(f"rayoptics v{rayoptics.__version__}")from opticalglass.glassfactory import create_glass
from opticalglass.rindexinfo import summary_plots
SiO2_url = 'https://refractiveindex.info/database/data-nk/main/SiO2/Malitson.yml'
SiO2 = create_glass(SiO2_url, "rindexinfo")
# summary_plots(SiO2)
# Input lens:
n1 = 1.444 # Fused Silica, Corning 7980 OF
g1 = SiO2
m1 = ''
R11 = 190.0 # (+) for CX, (-) for CC
t1 = 6.0
R12 = 50.0 # (+) for CX, (-) for CC
D1 = 12.0
# C to C gap between L1 and L2
d1 = 23.0
# Center lens:
n2 = 1.444 # Fused Silica, Corning 7980 OF
g2 = SiO2
m2 = ''
R21 = -25.0 # (+) for CX, (-) for CC
t2 = 10.0
R22 = -45.0 # (+) for CX, (-) for CC
D2 = 12.0
# C to C gap between L2 and L3
d2 = 27.0
# Output lens:
n3 = 1.444 # Fused Silica, Corning 7980 OF
g3 = SiO2
m3 = ''
R31 = -200.0 # (+) for CX, (-) for CC
t3 = 6.0
R32 = 55.0 # (+) for CX, (-) for CC
D3 = 12.0
BFL = 40
opm = OpticalModel()
sm = opm['seq_model']
osp = opm['optical_spec']
pm = opm['parax_model']
em = opm['ele_model']
pt = opm['part_tree']
ar = opm['analysis_results']
FFL = 10
bend_angle = -45
tilt_error = 0
L1_tilt = 5
L1_shift = 1
L2_tilt = -5
L2_shift = -1
L3_tilt = 5
L3_shift = 1
mirror_sign = 1
opm.system_spec.dimensions = 'mm'
opm.radius_mode = True
sm.do_apertures = False
osp['wvls'] = WvlSpec([(940.0, 1.0)], ref_wl=0)
#osp['pupil'] = PupilSpec(osp, key=['object', 'pupil'], value=2)
osp['pupil'] = PupilSpec(osp, key=['object', 'epd'], value=2.)
osp['fov'] = FieldSpec(osp, key=('object', 'angle'), value=1.,
flds=[0., 1., -1.],
is_relative=True)
sm.gaps[0].thi = 1.0e+10
opm.add_dummy_plane(t=BFL/2) # need to move the mirror by ~mirror_radius/sqrt(2)
sm.set_stop() # target for the rays so ray origin and direction is fixed
opm.add_dummy_plane()
sm.ifcs[sm.cur_surface].decenter = srf.DecenterData('decenter')
sm.ifcs[sm.cur_surface].decenter.euler[0] = bend_angle
opm.add_mirror(sd=10)
mirror_sign = -1*mirror_sign # change to -1 for odd bounces off mirror
mirror_surface = sm.cur_surface
sm.ifcs[sm.cur_surface].decenter = srf.DecenterData('dec and return')
sm.ifcs[sm.cur_surface].decenter.euler[0] = tilt_error
opm.add_dummy_plane(t=-BFL/2)
sm.ifcs[sm.cur_surface].decenter = srf.DecenterData('decenter')
sm.ifcs[sm.cur_surface].decenter.euler[0] = bend_angle
sm.add_surface([0, mirror_sign*FFL/2, 'air'])
sm.add_surface([mirror_sign*R11, mirror_sign*t1, g1, m1, D1/1])
sm.ifcs[sm.cur_surface].decenter = srf.DecenterData('decenter')
sm.ifcs[sm.cur_surface].decenter.euler[0] = L1_tilt
sm.ifcs[sm.cur_surface].decenter.dec[1] = L1_shift
sm.add_surface([-mirror_sign*R12, 0, 'air', D1/1])
sm.ifcs[sm.cur_surface].decenter = srf.DecenterData('reverse')
sm.ifcs[sm.cur_surface].decenter.euler[0] = L1_tilt
sm.ifcs[sm.cur_surface].decenter.dec[1] = L1_shift + mirror_sign*t1*math.sin(np.deg2rad(L1_tilt))
sm.gaps[sm.cur_surface].thi = mirror_sign*d1
sm.add_surface([mirror_sign*R21, mirror_sign*t2, g2, m2, D2/1])
sm.ifcs[sm.cur_surface].decenter = srf.DecenterData('decenter')
sm.ifcs[sm.cur_surface].decenter.euler[0] = L2_tilt
sm.ifcs[sm.cur_surface].decenter.dec[1] = L2_shift
sm.add_surface([-mirror_sign*R22, 0, 'air', D2/1])
sm.ifcs[sm.cur_surface].decenter = srf.DecenterData('reverse')
sm.ifcs[sm.cur_surface].decenter.euler[0] = L2_tilt
sm.ifcs[sm.cur_surface].decenter.dec[1] = L2_shift + mirror_sign*t2*math.sin(np.deg2rad(L2_tilt))
sm.gaps[sm.cur_surface].thi = mirror_sign*d2
sm.add_surface([mirror_sign*R31, mirror_sign*t3, g3, m3, D3/1])
sm.ifcs[sm.cur_surface].decenter = srf.DecenterData('decenter')
sm.ifcs[sm.cur_surface].decenter.euler[0] = L3_tilt
sm.ifcs[sm.cur_surface].decenter.dec[1] = L3_shift
sm.add_surface([-mirror_sign*R32, 0, 'air', D3/1])
sm.ifcs[sm.cur_surface].decenter = srf.DecenterData('reverse')
sm.ifcs[sm.cur_surface].decenter.euler[0] = L3_tilt
sm.ifcs[sm.cur_surface].decenter.dec[1] = L3_shift + mirror_sign*t3*math.sin(np.deg2rad(L3_tilt))
sm.gaps[sm.cur_surface].thi = mirror_sign*2*FFL
opm.update_model()
fig = plt.figure(FigureClass=InteractiveLayout, opt_model=opm, is_dark=False).plot() #, clip_rays=True)
|
Beta Was this translation helpful? Give feedback.
-
Beta Was this translation helpful? Give feedback.
-
|
Euler angles are tricky things. Each new primitive rotation is done on the previously modified coordinate system. In your case, the alpha tilt is applied and then the gamma rotation is applied around the newly rotated z axis. This has no effect on the tilt of the fold mirror. Euler angles are a brilliant solution for some problems and a nightmare for others. At ORA, we could be seen at our desks, our right hands formed into a set of coordinate axes, twisting our wrists to see how the coordinates changed for each rotation angle. This is one of the difficulties that LightTools was developed to solve. If what you're trying to achieve is to "clock" the fold mirrors around the z axis, this can be done by adding additional decenter surfaces before and after the fold mirrors and applying the gamma rotations to these decenters. If that isn't quite it, try to describe what you're trying to model and I can suggest approaches. Mike %matplotlib inline
from rayoptics.environment import *
from opticalglass.glassfactory import create_glass
from rayoptics.util import misc_math
import numpy as np
from scipy.spatial.transform import Rotationopm = OpticalModel()
sm = opm['seq_model']
osp = opm['optical_spec']
pm = opm['parax_model']
em = opm['ele_model']
pt = opm['part_tree']
ar = opm['analysis_results']
opm.system_spec.dimensions = 'mm'
opm.radius_mode = True
sm.do_apertures = False
osp['wvls'] = WvlSpec([(940.0, 1.0)], ref_wl=0)
# osp['pupil'] = PupilSpec(osp, key=['object', 'pupil'], value=2)
osp['pupil'] = PupilSpec(osp, key=['object', 'epd'], value=2.)
osp['fov'] = FieldSpec(osp, key=('object', 'angle'), value=1.,
flds=[0.,], # just 1 field to see what's going on
#flds=[0., 1., -1.],
is_relative=True)
def shiftr(euler, vec, noz=False):
rom = misc_math.euler2rot3d(euler)
# rom = Rotation.from_euler('xyz', euler, degrees=True).as_matrix()
xlat = np.dot(rom, vec)
if noz: xlat[2] = 0
return xlat
def addMirror(sd, m, euler_design=[-45,0,0], euler_clocking=[0,0,0], euler_off=[0,0,0], shift=[0,0,0], clock_after=False):
if not clock_after:
opm.add_dummy_plane()
sm.ifcs[sm.cur_surface].decenter = srf.DecenterData('decenter')
sm.ifcs[sm.cur_surface].decenter.euler = np.array(euler_clocking)
opm.add_dummy_plane()
sm.ifcs[sm.cur_surface].decenter = srf.DecenterData('decenter')
sm.ifcs[sm.cur_surface].decenter.euler = np.array(euler_design)
opm.add_mirror(sd=sd)
mir1 = sm.cur_surface
sm.ifcs[mir1].decenter = srf.DecenterData('dec and return')
sm.ifcs[mir1].decenter.euler = np.array(euler_off)
sm.ifcs[mir1].decenter.dec = shiftr(np.array(euler_design), np.array(shift))
opm.add_dummy_plane()
sm.ifcs[sm.cur_surface].decenter = srf.DecenterData('decenter')
sm.ifcs[sm.cur_surface].decenter.euler = np.array(euler_design)
if clock_after:
opm.add_dummy_plane()
sm.ifcs[sm.cur_surface].decenter = srf.DecenterData('decenter')
sm.ifcs[sm.cur_surface].decenter.euler = np.array(euler_clocking)
return -1*m, mir1 # mirror_sign = -1*mirror_sign # change to -1 for odd bounces off mirror
m_sign = 1
sm.gaps[0].thi = 1e10
opm.add_dummy_plane(t=10)
sm.set_stop() # target for the rays so ray origin and direction is fixedm_sign,mir1 = addMirror(5, m_sign, [-45,0,0], [0,0,-30], [0,0,0], [0,0,0])
sm.add_surface([0, m_sign*10, 'air'])
m_sign,mir2 = addMirror(5, m_sign, [45,0,0], [0,0,30], [0,0,0], [0,0,0], clock_after=True)
sm.add_surface([0, m_sign*20, 'air'])opm.update_model()sm.list_model()
sm.list_decenters(full=True)sm.list_gbl_tfrms()fig = plt.figure(FigureClass=InteractiveLayout, opt_model=opm, is_dark=False).plot() #, clip_rays=True)
|
Beta Was this translation helpful? Give feedback.
-
|
In short, I'm trying to create a system that I can use to make a complex geometry where not all the fold mirrors make 90 degree folds. There are some lenses here and there and some targets that get illuminated. I would like to put the system together as designed and then start adding some shifts and offsets and quantify the effects. I would like to have a mirror function that can be put in an arbitrary design position. I can use either intrinsic or extrinsic Euler angles as inputs. It appeared to me that Rayoptics uses extrinsic rotations because [-45,0,-90] gave me an optic with a normal that was the same as a [0,45,0] rotation. That is a normal that points in the original XZ plane with the vector [.7, 0, -.7]. This is itself confusing because the first rotation of [-45,0,0] shouldn't point along the +y axis but rather along the -y axis. I've been operating under the assumption that this is the default coordinate system with k initially pointing along the +z axis.
|
Beta Was this translation helpful? Give feedback.
-
|
Hello @justintorgerson, %matplotlib inlineimport numpy as np
import transforms3d as t3d
from scipy.spatial.transform import Rotation
from rayoptics.util import misc_mathdef print_rot(rot_mat):
print(np.array_str(rot_mat, precision=4, suppress_small=True))Use built-in ray-optics function euler2rot3drot0 = misc_math.euler2rot3d(np.array([-45, 0, 0]))
print_rot(rot0)beta_rot = misc_math.euler2rot3d(np.array([0, 45, 0]))
print_rot(beta_rot)gamma_rot = misc_math.euler2rot3d(np.array([0, 0, -90]))
print_rot(gamma_rot)rot1 = misc_math.euler2rot3d(np.array([-45, 0, -90]))
print_rot(rot1)print_rot(np.matmul(rot0.T, gamma_rot))The 2 matrices above should match, but they don't.Call euler2mat with the intrinsic
|
Beta Was this translation helpful? Give feedback.
-
|
Thank you, this helps. I wasn't aware of the optical convention either. I am still a bit confused about the nature of the sm.gbl_tfrms. Are they relative to the global coordinate system? Also, is propagation always along the +z intrinsic axis? If so, then a rotation about the intrinsic z-axis should not affect the propagation direction on-axis or show up in sm.gbl_tfrms. |
Beta Was this translation helpful? Give feedback.
-
|
The gbl_tfrms attribute in seq_model gives the position and orientation of every interface wrt the 1st interface. It is populated by the function compute_global_coords. This function has an argument, The z axis in the sequential model is used to closely model the propagation of light, but its really just coordinates systems positioning geometry in space. When a system has mirrors, you must flip the sign of the thickness to follow the propagation of light. After a mirror reflection, the light is now propagating along the -Z axis. The thickness for the gap following the mirror is negative so that the interface following the mirror is positioned to receive the reflected light. The |
Beta Was this translation helpful? Give feedback.
-
|
Is there a way to restart the model without restarting my python kernel? |
Beta Was this translation helpful? Give feedback.
-
|
Most of my issues have been resolved by moving to python 3.12 and getting the most recent version of rayoptics 0.95. For compatibility with some other packages, I was trying to stick with python 3.8. |
Beta Was this translation helpful? Give feedback.
-
|
How can I access the ray data for the rays that are traced automatically in the plots? |
Beta Was this translation helpful? Give feedback.
-
|
assume you have the following: The boundary defining rays for each field point are saved in a list of RayBundles called Each RayBundle has a dict called If you use ray fans, they're saved for each field in a list of RayFanBundles called Each RayFanBundle has a fan of rays in the fan called |
Beta Was this translation helpful? Give feedback.
assume you have the following:
The boundary defining rays for each field point are saved in a list of RayBundles called
ray_shapes.Each RayBundle has a dict called
rayset. It contains 5 rays, labeled '00', '+Y', '-Y', '+X', '-X'.If you use ray fans, they're saved for each field in a list of RayFanBundles called
rayfan_shapes.Each RayFanBundle has a fan of rays in the fan called
ray_fan. This is a RayFan and the list of rays is the fi…