contour_fourier_complex — SHAPE2D `descriptor` op

• 資料種類:pairs → matrix

• 呼叫: import fullseye as fs; fs.ledger.contour_fourier_complex(points, n_harmonics=None, parametrisation='index')(要直接呼叫實作,import fourierdesc; fourierdesc.contour_fourier_complex(points, n_harmonics=None, parametrisation='index');從台帳取用則 opsshape2d.get("contour_fourier_complex"))

用法

> 該運算子的說明尚無譯文,以下照原文給出。

閉輪郭を複素フーリエ級数の係数にする —— `z(t) = Σ c_k exp(2πi k t)`。

引数:

points: (N,2) の閉輪郭頂点。(row, col) 順(`fourierdesc` と同じ)。

n_harmonics: 残す次数 `K(|k| <= K)。None` なら表せる上限まで。

parametrisation: `"index"` = 点の添字を等間隔の媒介変数にする /

`"arclength"` = 弧長で等間隔に打ち直してから変換する。

返り値: `(2K+1, 3) の実配列(matrix)。行は [k, Re c_k, Im c_k]` で

`k は 0, 1, -1, 2, -2, …`(低い周波数から)。

★媒介変数の取り方を既定で黙って決めないのがこの op の肝。不均一に標本化された

輪郭を `"index"` で読むと、点が密なところに余分な時間が割り当てられ、

絵は似ているのに係数が別物になる(実測で `k=3` の係数が 47 倍ちがった)。

既定は `"index"`(FFT の素の意味)だが、形の記述として使うなら

`"arclength"` を選ぶこと。

★★`"arclength"` を選ぶと op の内部で輪郭を打ち直すので、返る係数は

渡した配列ではなく打ち直した輪郭を記述する。だから

:func:contour_fourier_truncation_energy の予言は「打ち直した輪郭に対して」

厳密で、呼び手が手元の配列で誤差を測ると合わない(実測: 512 点で 1.15 ずれた)。

手元の配列に対して厳密な予言が欲しいときは、先に弧長で打ち直してから

`"index"` で呼ぶこと。

★弧長の打ち直しは冪等ではない: 標本を弦で結ぶので、すでに等弧長の輪郭に

もう 1 度かけると点が最大 3.73 px 動き(平均 0.84 px)、長さがさらに 3.25 %

縮む。二度かけないこと。

詳細使用指南

• shape_description_2d 族使用指南

參考(範例資料・文獻)

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

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

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

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

• poc_one_stroke_epicycles — py -3.11 examples/poc_one_stroke_epicycles.py

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

contour_epicycle_chain · contour_fourier_truncation_energy

同類別(`descriptor`)

elliptic_fourier · reconstruct · invariants · normalize · descriptor_distance · fourier_smooth · from_xld · contour_epicycle_chain

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

contour_fourier_complex — SHAPE2D descriptor op

用法