observe_surface — OPTICS scene op

Data kinds: nonetable (an op determined by its arguments alone — it takes no image or data input)

Call: import optscene; optscene.observe_surface(material: 'str' = 'al', finish: 'str' = 'hairline', pitch_um: 'float' = 90.0, depth_um: 'float' = 1.0, roughness_um: 'float' = 0.08, focal_mm: 'float' = 25.0, working_distance_mm: 'float' = 200.0, f_number: 'float' = None, na: 'float' = None, pixel_um: 'float' = 3.45, resolution=(256, 256), wavelength_nm: 'float' = 550.0, illumination: 'str' = 'coaxial', tilt_deg: 'float' = 0.0, source_size_mm: 'float' = 40.0, supersample: 'int' = 2, exposure: 'str' = 'auto', defects: 'dict' = None, seed: 'int' = 0) -> 'dict' (or opsoptics.get("observe_surface"))

Usage

Assemble an observation optic and shoot the finished surface of a named material (a one-line entry point).

> The detailed description below is the original text — the summary and the headings are translated.

ユーザー要望(2026-09-05)「基本的な観察用の光学系をレイアウトした場合を想定して、

それらのパラメータを与えて、アルミなどの素材を指定した上で、ヘアライン画像を

生成できるといいな。最初は同軸照明でも良い」。

光学系は NA でも F 値でも指定できる。作動距離・焦点距離・画素ピッチ・波長から

:func:optical_budget が倍率・エアリー半径・被写界深度・分解限界を出し、その

まま撮像に効く(回折ぼけは :func:diffraction_blur、ぼけの深さ方向は

:func:depth_of_field)。素材は `glassmirror.METALS`、仕上げは

`surface_finish` の種類。

`illumination` は coaxial / ring / dome / bar / backlight。同軸は面が鏡なら

器具が映って明るく(明視野)、低角のリングは映らないので暗い(暗視野)。

返り値 dict: `image(線形 RGB)/ camera / scene / light` /

`budget(上の光学バジェット)/ defect_mask(defects` を渡したとき)。

既定の引数は 117 秒かかる。 「一行で使える入口」と書いてあるが、

既定 `resolution=(256, 256) / supersample=2` は 2 分の計算である

(2026-09-06 実測、この機械):

============== ============= =============

resolution supersample 1 supersample 2

============== ============= =============

32 x 32 0.25 s —

64 x 64 0.69 s —

128 x 128 2.40 s 19.10 s

256 x 256 19.12 s 116.93 s

============== ============= =============

画素あたりの費用が一定ではない(32→256 で画素は 64 倍なのに時間は 76 倍)。

内訳は `_light_background` が 16.2 秒中の 11.0 秒で、面光源 196 点 x 鏡面

候補 25 回のループ。発光点をまとめて処理する書き換えは試したが 1.3 倍

にしかならず(律速は Python のループではなくメモリ帯域)、まだ解いていない

—— `docs/KNOWN_ISSUES.md` に残してある。

試すときは小さく呼ぶこと: ``observe_surface(resolution=(64, 64),

supersample=1)`` なら 0.69 秒で、仕上げの見え方の違いは十分わかる。

連鎖ファザーもこの理由で 32 x 32 / supersample 1 に固定してある

(`typed_catalog.OP_PARAM_HINTS`)。

Family-wide input contract (fail-closed)

Every optics op validates its input before computing (nothing slips through silently):

Units are baked into the argument name_mm / _um / _deg / _mrad. Confusing mm with µm does not crash; it yields a plausible-looking wrong answer, so the name prevents it. Nothing here guesses the unit from the magnitude.

• **Strings raise ValueError** — float('50') succeeds, so an unparsed configuration value would slip through as a length (measured: thin_lens('50', '200') returned a plausible 66.667 mm). bool is refused too, as the implicit promotion True == 1.

• **complex / masked arrays raise ValueError (real-valued slots only; silently dropping the imaginary part or peeling off the mask is refused). NaN/Inf raises ValueError on every input.**

Division by zero and its relatives are refused by name: focal length 0, radius of curvature 0, refractive index <= 0, a fully opaque aperture (all zeros, so the normalisation is 0/0), a PSF whose sum is <= 0, a Stokes vector with S0 = 0, and an object sitting at the front focal point (the image is at infinity).

Only two ops return a non-finite value, and both state it as a contract: depth_of_field returns far_mm = inf beyond the hyperfocal distance (that is what the hyperfocal distance means), and gaussian_beam returns wavefront_radius_mm = inf at the waist (the radius of curvature of a plane wavefront). Both also return a finite companion (far_is_infinite / curvature_per_mm). **Any other silent NaN/Inf is detected internally and raises ValueError** — "float64 overflowed" and "the answer is infinite" are different claims, so the first is never returned wearing the face of the second.

Size caps: generated grids are capped by optics.MAX_GRID (4096); supplied fields/PSFs/apertures by optics.MAX_FIELD_ELEMENTS (2^24); ABCD element chains by optics.MAX_SYSTEM_ELEMENTS (1024); Zernike by MAX_ZERNIKE_TERMS (512) / MAX_ZERNIKE_ORDER (40) / MAX_ZERNIKE_BASIS (2^25). This closes, fail-closed, the paths where a small argument triggers a huge internal allocation (measured: n_max=40 × 4096² needs 108 GB).

Physically impossible states are refused too: a Stokes vector with degree of polarisation > 1, negative transmittance, negative intensity, and invalid Zernike indices such as n-|m| odd.

Detailed usage guide

optics_imaging family guide

Background guides (the physics and conventions behind this op)

dataset_conventions — 学習データセット規約の知識 — COCO / YOLO / VOC と外観検査での落とし穴

mv_cameras — 産業用カメラメーカー(センサとの紐付け・ラインスキャン / TDI)

mv_illumination_practice — 照明の実務知識 — 波長・偏光・点灯方式・外光・安全

mv_image_sensors — 産業用イメージセンサ(現行品中心)

virtual_machine_vision — 仮想マシンビジョン — パラメータの洗い出しとオブジェクト模型

References (sample data, literature)

• Sample-data catalog (download URLs / licences) — 2-D uses skimage.data (BSD/public domain) plus synthetic images; 3-D lists download URLs for real data sources (Stanford, PDS, …).

• Operator provenance and references — the sources of the research/methods this op family came from.

• The canonical algorithm (author, year) and its uses are named in the family usage guide above.

Runnable examples (verified samples that actually call this op)

vision_layout_from_catalogpy -3.11 examples/vision_layout_from_catalog.py

Ops the type connects to (they accept table as input)

abcd_matrix · wavefront_stats · paraxial_trace · seidel_coefficients · spot_stats · tolerance_analysis · wavefront_from_opd · spot_diagram

Same category (scene)

scene_material · scene_plane · scene_sphere · scene_box · scene_cylinder · surface_defect · surface_finish · random_defects


*Provenance: optscene.py — OPTICS operator registry. This per-op note is generated by tools/opdocs.py md (do not hand-edit).*

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