基于由塞巴斯蒂安Endrikat启发//www.tatmou.com/matlabcentral/fileexchange/52343-find-vortices-in-velocity-fields

底层数学提出用L Graftieaux“结合PIV，POD和涡流识别算法不稳定湍流涡旋的研究流”https://iopscience.iop.org/article/10.1088/0957-0233/12/9/307/pdf
非常感谢！

塞巴斯蒂安的代码是很慢的，但是，因为它使用嵌套的for循环4。生产迦玛地图可以通过认识到在Graftieaux的论文中描述的操作可以重新安排成两个卷积更有效地在Matlab进行。此代码这样做，基本上是这样。它还通过拟合样条曲线找到最大伽马1，因此没有由矢量场的网格分辨率的限制。

其缺点是，它会发现只有一个涡核，这可能不是在某些应用中的多个涡流形成不理想。一个也可以，但是，使用第二个输出（G1）执行自己的峰值寻找/阈值处理来查找多个漩涡核心。

细节写了上https://zigunov.com/2019/07/08/finding-vortex-cores-with-gamma-1/

希望这是有用的！
-FZ