def _test(self): p = '/Users/ross/Pychrondata_demo/data/snapshots/scan6/007.jpg' g = Graph() g.new_plot() for scan_i, z, idxs in [ # 1, # 2, # 3, 4, 5, # (6, [20, 30, 40, 50, 60, 70, 80, 90, 100, ], # [2, 3, 4, 5, 6, 7, 8, 9, 10] # ), (6, [10], [1] ), # (6, [100, 90, 80, 70, 60, 50, 40, 30, 20], # [11, 12, 13, 14, 15, 16, 17, 18, 19] # ) ]: dxs = [] zs = [] root = '/Users/ross/Pychrondata_demo/data/snapshots/scan{}'.format(scan_i) for zi, idx in zip(z, idxs): pn = os.path.join(root, '{:03n}.jpg'.format(idx)) d = load_image(pn) dx = self._calculate_spacing(d) dxs.append(dx) zs.append(zi) g.new_series(zs, dxs, type='scatter') coeffs = polyfit(zs, dxs, 2) print 'parabolic intercept {}'.format(coeffs[-1]) xs = linspace(0, max(zs)) ys = polyval(coeffs, xs) g.new_series(xs, ys) fitfunc = lambda p, x: p[0] * exp(p[1] * x) + p[2] lr = LeastSquaresRegressor(fitfunc=fitfunc, initial_guess=[1, 0.1, 0], xs=zs, ys=dxs ) xs = linspace(0, max(zs)) ys = lr.predict(xs) print 'exponential intercept {}'.format(lr.predict(0)) g.new_series(xs, ys) invoke_in_main_thread(g.edit_traits)
def _test(self): p = '/Users/ross/Pychrondata_demo/data/snapshots/scan6/007.jpg' g = Graph() g.new_plot() for scan_i, z, idxs in [ # 1, # 2, # 3, 4, 5, # (6, [20, 30, 40, 50, 60, 70, 80, 90, 100, ], # [2, 3, 4, 5, 6, 7, 8, 9, 10] # ), (6, [10], [1]), # (6, [100, 90, 80, 70, 60, 50, 40, 30, 20], # [11, 12, 13, 14, 15, 16, 17, 18, 19] # ) ]: dxs = [] zs = [] root = '/Users/ross/Pychrondata_demo/data/snapshots/scan{}'.format( scan_i) for zi, idx in zip(z, idxs): pn = os.path.join(root, '{:03n}.jpg'.format(idx)) d = load_image(pn) dx = self._calculate_spacing(d) dxs.append(dx) zs.append(zi) g.new_series(zs, dxs, type='scatter') coeffs = polyfit(zs, dxs, 2) print 'parabolic intercept {}'.format(coeffs[-1]) xs = linspace(0, max(zs)) ys = polyval(coeffs, xs) g.new_series(xs, ys) fitfunc = lambda p, x: p[0] * exp(p[1] * x) + p[2] lr = LeastSquaresRegressor(fitfunc=fitfunc, initial_guess=[1, 0.1, 0], xs=zs, ys=dxs) xs = linspace(0, max(zs)) ys = lr.predict(xs) print 'exponential intercept {}'.format(lr.predict(0)) g.new_series(xs, ys) invoke_in_main_thread(g.edit_traits)
def _test2(self): self.test_image.setup_images(3, # (475, 613) (640, 480) ) root = '/Users/ross/Pychrondata_demo/data/snapshots/scan6' p = os.path.join(root, '008.jpg') im = load_image(p) im = grayspace(im) self.test_image.set_image(im) nim = zeros_like(im) nim[((im > 100) & (im < 200))] = 255 self.test_image.set_image(nim, 1)
def _test2(self): self.test_image.setup_images( 3, # (475, 613) (640, 480)) root = '/Users/ross/Pychrondata_demo/data/snapshots/scan6' p = os.path.join(root, '008.jpg') im = load_image(p) im = grayspace(im) self.test_image.set_image(im) nim = zeros_like(im) nim[((im > 100) & (im < 200))] = 255 self.test_image.set_image(nim, 1)