Esempio n. 1
0
import os

from odysseus import imageio, imageprocess
from odysseus.fitfermions import fit_img, find_ellipticity

dirname = 'datafiles'
#dirname = 'c:\\Data\\2008-12-12' # double backslashes on windows

# import a single image (for finding each TIF file in a directory,
# see example_filetools.py)
fname = 'raw9.11.2008 7;35;23 PM.TIF'
img_name = os.path.join(dirname, fname)
rawframes = imageio.import_rawframes(img_name)

#calculate the transmission and OD images
transimg, odimg = imageprocess.calc_absimage(rawframes)
# set the ROI
transimg = transimg[120:350, 50:275]

# find the ellipticity if the expansion is not spherically symmetric
ellip = find_ellipticity(transimg)
# do the fit (this pops up an image, unless you set showfig=False)
ToverTF, N, ans = fit_img(transimg, elliptic=(ellip, 0), showfig=True)
Esempio n. 2
0
 def test_find_ellipticity(self):
     tol = 2e-3
     assert find_ellipticity(self.img, tol=tol) - self.ellip < tol
Esempio n. 3
0
 def test_find_ellipticity(self):
     tol = 2e-3
     assert find_ellipticity(self.img, tol=tol) - self.ellip < tol
Esempio n. 4
0
import os

from odysseus import imageio, imageprocess
from odysseus.fitfermions import fit_img, find_ellipticity


dirname = 'datafiles'
#dirname = 'c:\\Data\\2008-12-12' # double backslashes on windows

# import a single image (for finding each TIF file in a directory,
# see example_filetools.py)
fname = 'raw9.11.2008 7;35;23 PM.TIF'
img_name = os.path.join(dirname, fname)
rawframes = imageio.import_rawframes(img_name)

#calculate the transmission and OD images
transimg, odimg = imageprocess.calc_absimage(rawframes)
# set the ROI
transimg = transimg[120:350, 50:275]

# find the ellipticity if the expansion is not spherically symmetric
ellip = find_ellipticity(transimg)
# do the fit (this pops up an image, unless you set showfig=False)
ToverTF, N, ans = fit_img(transimg, elliptic=(ellip, 0), showfig=True)