6 from scipy.ndimage.measurements
import center_of_mass
11 os.getenv(
"COMPASS_ROOT") +
"/shesha/data/par/par4bench/scao_sh_16x16_8pix.py")
12 sup.config.p_controllers[0].delay = 0
18 os.getenv(
"COMPASS_ROOT") +
"/shesha/data/par/par4bench/scao_sh_16x16_8pix.py")
19 sup2.config.p_controllers[0].delay = 0
25 rtc.add_centroider(sup.context, sup.config.p_wfs0._nvalid,
26 sup.config.p_wfs0.npix / 2 - 0.5, sup.config.p_wfs0.pixsize,
False, 0,
"cog")
27 rtc.add_controller(sup.context, sup.config.p_wfs0._nvalid, sup.config.p_wfs0._nvalid * 2,
28 sup.config.p_controller0.nactu, sup.config.p_controller0.delay, 0,
29 "generic", idx_centro=np.zeros(1), ncentro=1)
30 centro = rtc.d_centro[0]
31 control = rtc.d_control[0]
32 rtc.d_centro[0].set_npix(sup.config.p_wfs0.npix)
33 xvalid = np.array(sup.rtc.rtc.d_centro[0].d_validx)
34 yvalid = np.array(sup.rtc.rtc.d_centro[0].d_validy)
35 rtc.d_centro[0].load_validpos(xvalid, yvalid, xvalid.size)
36 cmat = sup.rtc.get_command_matrix(0)
37 rtc.d_control[0].set_cmat(cmat)
38 rtc.d_control[0].set_gain(sup.config.p_controller0.gain)
52 sup.dms.set_command(np.array(sup._sim.rtc.d_control[0].d_voltage))
53 sup._sim.raytraceTar(0,
"all")
54 sup._sim.comp_tar_image(0,compLE=
False)
55 sup._sim.compStrehl(0)
56 SR32.append(sup.get_strehl(0)[0])
57 frame = sup2.get_wfs_image(0)
59 rtc.d_centro[0].load_img(frame, frame.shape[0])
60 rtc.d_centro[0].calibrate_img()
64 sup2.dms.set_command(np.array(rtc.d_control[0].d_voltage))
65 sup2._sim.raytraceTar(0,
"all")
66 sup2._sim.comp_tar_image(0,compLE=
False)
67 sup2._sim.compStrehl(0)
68 SR16.append(sup2.get_strehl(0)[0])
71 for _
in tqdm(range(niter)):