def _synth_field(fl_i):
    with synth.Synth(n_channels=n_channels, n_cycles=n_cycles, dim=dim) as s:
        peaks = synth.PeaksModelGaussianCircular(n_peaks=n_peaks).locs_randomize().widths_uniform(psf_width).amps_constant(gain)
        synth.CameraModel(bg_mean=bg_mean, bg_std=bg_std)
        synth.HaloModel()
        synth.IlluminationQuadraticFalloffModel()
    chcy_ims = s.render_chcy(0)
    for ch_i in range(chcy_ims.shape[0]):
        for cy_i in range(chcy_ims.shape[1]):
            np.save(str((source_folder / f'area_{fl_i:03d}_cell_000_{ch_i:03d}nm_{cy_i:03d}.npy')), chcy_ims[(ch_i, cy_i)])
def _synth_field(fl_i):
    with synth.Synth(n_channels=n_channels, n_cycles=n_cycles, dim=dim) as s:
        peaks = synth.PeaksModelGaussianCircular(n_peaks=n_peaks).locs_randomize().widths_uniform(psf_width).amps_constant(gain)
        synth.CameraModel(bg_mean=bg_mean, bg_std=bg_std)
        synth.HaloModel()
        synth.IlluminationQuadraticFalloffModel()
    chcy_ims = s.render_chcy(0)
    for ch_i in range(chcy_ims.shape[0]):
        for cy_i in range(chcy_ims.shape[1]):
            np.save(str((source_folder / f'area_{fl_i:03d}_cell_000_{ch_i:03d}nm_{cy_i:03d}.npy')), chcy_ims[(ch_i, cy_i)])