piv_divergence — PIV field op

数据种类:flow2dimage2d

调用: import fullseye as fs; fs.ledger.piv_divergence(flow, spacing=1.0)(要直接调用实现,import pivops; pivops.piv_divergence(flow, spacing=1.0);从台账取用则 opspiv.get("piv_divergence"))

用法

> 该算子的说明尚无译文,以下照原文给出。

発散 `d(dy)/dy + d(dx)/dx`。非圧縮なら 0 —— 独立な検算に使える。

式: `gradient(flow[0], spacing, axis=0) + gradient(flow[1], spacing, axis=1)`

(`numpy.gradient`: 内部は中央差分、格子の端は片側差分)。

`flow[0] が行方向の変位 dyflow[1] が列方向の変位 dx`、

`x = 列、y` = 行。行が下向きでも発散の符号は座標の向きに依らない

(渦度 `piv_vorticity` とは違い、反転しない)。正 = 湧き出し、負 = 吸い込み。

• `flow: (2, h, w)、成分 (dy, dx) [px/frame] (piv_cross_correlate` /

`piv_multipass の返り値)。h, w >= 2` でないと微分できないので

`ValueError`(窓が大きすぎて 1 本しかベクトルが無い場合)。

有限性は検査しない ―― `piv_replace_outliers(method="nan")` の NaN は

隣接 2 セルへ広がる。

• `spacing: ベクトル間隔 [px]、正の有限値。info["step"]` を渡す

(既定 1.0 のままだと単位が「1/ベクトル間隔」になる)。

• 返り値: `(h, w) float64、単位 1/frame(spacing` が px のとき)。

`piv_to_velocity で m/s にした場と spacing` を m にすれば 1/s。

2-D PIV で平面外の速度成分があると非零になるので、「非圧縮なのに発散が

大きい」は面外流か外れ値(`piv_outlier_mask`)の指標。

`piv_velocity_gradient` は 4 成分すべてを返す。

详细使用指南

piv_displacement 族使用指南

参考(示例数据・文献)

• 示例数据目录(下载 URL / 许可证) —— 2-D 用 skimage.data(BSD/公有领域)加合成图,3-D 给出真实数据源(Stanford/PDS 等)的下载 URL。

• 算子来历与参考文献 —— 该算子族所依据的研究/方法出处。

• 算法的正典(作者・年份)与用途见上面的族使用指南

可运行的示例(实际调用该算子并已验证的样例)

piv_flow_from_particlespy -3.11 examples/piv_flow_from_particles.py

类型可衔接的下一个算子(可接受 image2d 作为输入)

piv_cross_correlate · piv_multipass · piv_deform_pass · strain_from_displacement · correlation_quality · speckle_quality

同类别(field)

piv_vorticity · piv_flow_magnitude · piv_to_velocity · piv_velocity_gradient · piv_q_criterion · piv_swirling_strength · piv_strain_rate


*Provenance: pivops.py — PIV 算子登记表。本条目由 tools/opdocs.py md 自动生成(请勿手工编辑)。*

© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.