def test_auto_count_n(self):
out = core.auto_count(10, n = 5)
self.assertTrue(np.allclose(out,auto_count_10[0:5]))
out = core.auto_count(t1, n = 5)
self.assertTrue(np.allclose(out,auto_count_t1[:5]))
out = core.auto_count(t1, aout = out)
self.assertTrue(np.allclose(out,2*auto_count_t1[:5]))
def test_auto_count(self):
out = core.auto_count(10)
self.assertTrue(np.allclose(out,auto_count_10))
out = core.auto_count(t1)
self.assertTrue(np.allclose(out,auto_count_t1))
out = core.auto_count(t1, aout = out)
self.assertTrue(np.allclose(out,auto_count_t1*2))
"""Random triggered single camera video
"""
from cddm.core import auto_count
from cddm.sim import simple_brownian_video, create_random_time, adc
from cddm.viewer import VideoViewer
from cddm.video import crop, multiply
from examples.paper.simple_video.conf import NFRAMES_RANDOM, SIMSHAPE, BACKGROUND, DELTA, INTENSITY, PERIOD_RANDOM,\
NUM_PARTICLES, SIGMA, SHAPE, DUST1_PATH, BIT_DEPTH, VMAX, SATURATION, READOUT_NOISE, NOISE_MODEL, APPLY_DUST, DT_RANDOM
import matplotlib.pyplot as plt
import numpy as np
while True:
#make valid random time (all passible times present)
t = create_random_time(NFRAMES_RANDOM, n=PERIOD_RANDOM)
count = auto_count(t, NFRAMES_RANDOM)
#make sure all times are present, if so break, else try again
if np.all(count):
break
#t = create_random_time(NFRAMES_RANDOM, n=PERIOD, dt_min = PERIOD//2)
#: this creates a brownian motion frame iterator.
#: each element of the iterator is a tuple holding a single numpy array (frame)
video = simple_brownian_video(t,
shape=SIMSHAPE,
background=BACKGROUND,
num_particles=NUM_PARTICLES,
dt=DT_RANDOM,
sigma=SIGMA,
delta=DELTA,