fit_poly_surface — 3D surface_fit op

데이터 종류: image2d × image2d × image2dpoly_surface

호출: import fullseye as fs; fs.ledger.fit_poly_surface(x, y, z, degree=2)(구현을 직접 호출하려면 import match3d; match3d.fit_poly_surface(x, y, z, degree=2), 원장에서 가져오려면 ops3d.get("fit_poly_surface"))

사용법

산포 (x,y,z) → z=f(x,y) 다항 최소제곱. 반환값 model(coef/powers/degree/rms/pv).

> 아래 상세 설명은 원문입니다 —— 요약과 제목은 번역되어 있습니다.

基底 `{x^i·y^j : i + j <= degree}(項数 (degree+1)(degree+2)/2`。degree=1 で 3 項の平面、

2 で 6 項の 2 次曲面)を `lstsq で当てる。x, y, z` は同じ要素数なら形は問わない(内部で

ravel。格子なら `np.mgrid` の出力をそのまま)。

返り値 dict: `coef (T,) 係数、powers は各係数の (i, j)(項 x**i * y**j`)、

`degreerms は残差 RMS、pv` は残差の peak-to-valley(max − min)。単位は z。

• 点数が項数より少ないと最小ノルム解が黙って返る。x, y の桁が大きいと高次で条件が悪くなる

(座標を中心化・正規化してから)。NaN の検証は無い。

後段: `eval_poly_surface(model, x, y) で任意点を評価。格子の高さ場なら surface_form_error`

/ `background_flatten` がこれを内部で呼ぶ。

배경 지식 가이드(이 연산자 앞에 있는 물리·규약)

blas_threads_and_memory — 行列分解が遅い理由の知識 — BLAS スレッド・キャッシュ・メモリ配置

참고(샘플 데이터·문헌)

• 샘플 데이터 카탈로그(DL URL / 라이선스) —— 2-D 는 skimage.data(BSD/public)+ 합성, 3-D 는 실데이터 소스(Stanford/PDS 등)의 DL URL.

• 연산자의 내력·참고문헌 —— 이 연산자 족의 바탕이 된 연구/기법의 출처.

실행 가능한 예제(이 연산자를 실제로 호출하는 검증된 샘플)

contours_to_terrainpy -3.11 examples_3d/contours_to_terrain.py

타입이 이어지는 다음 연산자(poly_surface 를 입력으로 받는 것)

fuse_to_voxel · eval_poly_surface

같은 카테고리(surface_fit)

eval_poly_surface · surface_form_error · background_flatten


*Provenance: match3d.py — 3D 연산자 레지스트리. 이 op 노트는 tools/opdocs.py md 가 자동 생성합니다(직접 편집하지 마세요).*

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