Nice example, thank you Alex!
It seems that the multi-grid window deformation algorithm `windef` is a more preferred method for my images, which I did not use before. I compare below my old method (as described in the first tutorial based on extended_area_search_piv), simple_piv and windef.
I tried to find more information about windef in the documentation site (readthedocs), but only found a short tutorial (
https://openpiv.readthedocs.io/en/latest/src/windef.html), which does not say much about how it works and what advantages to expect. windef also generates many intermediate results (images) which are supposed to be informative of what is going on in the backend, but I don't find a guide about how to use this information. `windef` might be a good place to improve the documentation.
I also notice that when I do the quiver plot using `matplotlib.pyplot.quiver(x, y, u, v)` (left panel), the vertical velocity components all look inverted. As a result, I have to do `quiver(x, y, u, -v)` to make the vector field consistent with the motion (shown below). I might misunderstand the convention, so it would be helpful if someone can help me understand this.
I also notice something I don't understand in the last two output images of windef: the y-axis direction is flipped, from [350, 0] to [0, 350], but the velocity field remains the same.