コード例 #1
0
ファイル: run_SubaruFrontend.py プロジェクト: jessmos/MEDIS
            subplt_cols=sp.spectra_cols,
            vlim=(1e-4, 1e-1),
            dx=fp_sampling)

    # Plotting Timeseries in White Light
    if sp.show_tseries:
        img_tseries = np.sum(focal_plane, axis=1)  # sum over wavelength
        view_timeseries(img_tseries,
                        title=f"White Light Timeseries\n"
                        f"AO={tp.use_ao}",
                        subplt_cols=sp.tseries_cols,
                        logZ=True,
                        vlim=(1e-6, 1e-3))
        # dx=fp_sampling[0])

    # Plotting Selected Plane
    if sp.show_planes:
        vlim = [(None, None), (None, None), (None, None), (None, None),
                (None, None)]
        # vlim = [(None,None), (None,None), (None,None), (1e-2,1e-1)]  # (7e-4, 6e-4)
        logZ = [True, False, False, True]
        if sp.save_list:
            plot_planes(cpx_sequence,
                        title=f"White Light through Optical System",
                        subplt_cols=2,
                        vlim=vlim,
                        logZ=logZ,
                        dx=sampling)

    test = 1
コード例 #2
0
ファイル: run_SCExAO-DM.py プロジェクト: jessmos/MEDIS
    # Plotting Timeseries in White Light
    if sp.show_tseries:
        view_timeseries(img_tseries, cdi, title=f"White Light Timeseries\n"  #f"AO={tp.use_ao}. CDI={cdi.use_cdi}"
                                                f'n_probes={cdi.n_probes}, phase intervals = {cdi.phs_intervals/np.pi:.2f}'\
                                                + r'$\pi$',
                        subplt_cols=sp.tseries_cols,
                        logZ=True,
                        vlim=(1e-7, 1e-4),
                        dx=fp_sampling[0])

    # Plotting Selected Plane
    if sp.show_planes:
        # vlim = [(None, None), (None, None), (None, None), (1e-7,1e-3), (1e-7,1e-3), (1e-7,1e-3)]
        vlim = [(None, None), (1e-7, 1e-4), (1e-8, 1e-4),
                (None, None)]  # (1e-2,1e-1) (7e-4, 6e-4)
        logZ = [False, True, True, True, True, True]
        if sp.save_list:
            plot_planes(cpx_sequence,
                        title=f"White Light through Optical System"
                        f"\nFirst Timestep",
                        subplt_cols=2,
                        vlim=vlim,
                        logZ=logZ,
                        first=True,
                        dx=sampling)

    plt.show()

##