def test_minimal_piv(self):
params = WorkPIV.create_default_params()
# for a very short computation
params.piv0.shape_crop_im0 = 32
params.piv0.grid.overlap = -3
# still buggy
# params.piv0.nb_peaks_to_search = 2
params.multipass.number = 2
params.fix.displacement_max = 2
params.fix.threshold_diff_neighbour = 2
piv = WorkPIV(params=params)
with stdout_redirected():
result = piv.calcul(self.serie)
result.piv0.save(self.path_tmp)
result.save(self.path_tmp)
path_file = next(self.path_tmp.iterdir())
MultipassPIVResults(path_file)
result.save(self.path_tmp, "uvmat")
light = result.make_light_result()
light.save(self.path_tmp)
path_file = next(self.path_tmp.glob("*light*"))
LightPIVResults(str_path=str(path_file))
def test_piv_list(self):
params = WorkPIV.create_default_params()
# for a very short computation
params.piv0.shape_crop_im0 = [33, 44]
params.piv0.grid.overlap = -3
# params.piv0.nb_peaks_to_search = 2
params.multipass.use_tps = False
piv = WorkPIV(params=params)
with stdout_redirected():
result = piv.calcul(self.serie)
piv0 = result.piv0
im0, im1 = piv0.get_images()
with stdout_redirected():
DisplayPIV(im0, im1, piv0)
d = DisplayPIV(im0, im1, piv0, show_interp=True, hist=True)
d.switch()
d.select_arrow([0], artist=d.q)
d.select_arrow([0], artist=None)
event = MyObj()
event.key = "alt+h"
event.inaxes = None
d.onclick(event)
event.inaxes = d.ax1
event.key = "alt+s"
d.onclick(event)
event.key = "alt+left"
d.onclick(event)
event.key = "alt+right"
d.onclick(event)
event.artist = None
d.onpick(event)
event.artist = d.q
event.ind = [0]
d.onpick(event)
from fluidimage import SeriesOfArrays from fluidimage.works.piv import WorkPIV params = WorkPIV.create_default_params() params.multipass.number = 2 params.multipass.use_tps = True params.piv0.shape_crop_im0 = 32 params.piv0.displacement_max = 5 params.fix.correl_min = 0.2 params.fix.threshold_diff_neighbour = 8 params.mask.strcrop = '30:250, 100:' work = WorkPIV(params=params) path = '../../image_samples/Oseen/Images' # path = '../../image_samples/Karman/Images' series = SeriesOfArrays(path, 'i+1:i+3') serie = series.get_serie_from_index(0) piv = work.calcul(serie) # piv.display(show_interp=True, scale=0.3, show_error=True) piv.display(show_interp=False, scale=1, show_error=True) # result.save()