monogenic_orientation — QUAT riesz op

資料種類:qimageimage2d

呼叫:import quatimage; quatimage.monogenic_orientation(qimage, display: 'bool' = False) -> 'np.ndarray'(或 opsquat.get("monogenic_orientation"))

用法

單演訊號的局部方向 `atan2(R2, R1)`。→ (H, W)。

> 以下的詳細說明為原文 —— 摘要與標題已翻譯。

Radians in `[0, pi): an orientation is defined modulo pi` (a grating at

10 degrees and one at 190 degrees are the same grating), and the value is

folded into that range rather than left in `(-pi, pi]` where the same

structure would read as two different numbers on either side of a contrast

reversal. `display=True maps it to [0, 1]`.

Continuous, not quantised — the angle is read directly from two filters,

for any angle, where a steerable bank with `K` orientations interpolates

between its `K`. Measured against eight grid-exact grating orientations the

error is at most 3.6e-15 rad, including the obliques. (Whether that

buys anything downstream is a separate question, and the measured answer is

mostly *no* — see :func:riesz_displacement.)

Where it is undefined, and the mask is not the one you expect. The

orientation dies where the *Riesz vector* dies, which is at every

even-symmetric point — local phase 0 or pi, the crest of a bright or dark

line — and the amplitude is at full strength there. Measured on a 45-degree

grating, the worst orientation error over the whole frame is 0.2764 rad, at a

pixel where `|R| = 6.8e-16 and :func:monogenic_amplitude` reads

`1.0000`. So masking on the amplitude does not protect you; mask on

`hypot(q[..., 1], q[..., 2])`, the Riesz magnitude. With that mask the

error over the same eight orientations is at most 3.6e-15 rad.

Where the Riesz vector is exactly zero, `atan2(0, 0) = 0` is returned —

a *value*, not a measurement.

Raises `ValueError`: the input is not a valid quaternion field, or its

`k` component is non-zero; *display* is not a bool.

詳細使用指南

quaternion_monogenic 族使用指南

參考(範例資料・文獻)

• 範例資料目錄(下載 URL / 授權) —— 2-D 用 skimage.data(BSD/公有領域)加合成圖,3-D 給出真實資料源(Stanford/PDS 等)的下載 URL。

• 運算子來歷與參考文獻 —— 該運算子族所依據的研究/方法出處。

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

可執行的範例(實際呼叫該運算子並已驗證的樣例)

quaternion_monogenicpy -3.11 examples/quaternion_monogenic.py

型別可銜接的下一個運算子(可接受 image2d 作為輸入)

riesz_transform · monogenic_signal

同類別(riesz)

riesz_transform · monogenic_signal · monogenic_amplitude · monogenic_phase


*Provenance: quatimage.py — QUAT 運算子登記表。本條目由 tools/opdocs.py md 自動產生(請勿手動編輯)。*

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