oct op• Data kinds: image2d → table
• Call: import fullseye as fs; fs.ledger.oct_bscan_boundaries(image, n_boundaries=2, um_per_px=1.0, n_group=None, min_gap_px=4, max_jump_px=2, jump_penalty=0.0, smooth_sigma=1.0, mode='peak', top_px=None, bottom_px=None, reference_rows=None, offset_px=None) (implementation को सीधे call करना हो तो import octscan; octscan.oct_bscan_boundaries(image, n_boundaries=2, um_per_px=1.0, n_group=None, min_gap_px=4, max_jump_px=2, jump_penalty=0.0, smooth_sigma=1.0, mode='peak', top_px=None, bottom_px=None, reference_rows=None, offset_px=None), registry से लेना हो तो opsoptics.get("oct_bscan_boundaries"))
B-scan image पर interfaces को पार्श्व दिशा में trace करना (image2d → table)। instrument के output चित्र पर भी काम करता है।
> नीचे का विस्तृत विवरण मूल text में है — सारांश और headings का अनुवाद हो चुका है।
各界面 = 列ごとに 1 行を選ぶ道で、隣の列との行の差 ≤ max_jump_px、費用 = −(明るさの z 値)の和 + jump_penalty × |行の差| が最小のもの(動的計画法。網膜の層分割の Chiu ら 2010、塗膜の Dong ら 2016 と同じ「最短路」の考え方の 1 本ずつ版)。mode="peak" は明るさ(界面は明るい線)、"edge" は下向きの明るさの増え方(暗 → 明、散乱の多い層の上の面)。見つけた道の上下 min_gap_px は次の道で使わない。
探す範囲: top_px / bottom_px(数か列ごとの (W,))、または reference_rows(既に見つけた面、例えば表面)からの相対 offset_px = (下限, 上限)。★相対で探すと、面と平行な界面に罰が掛からず、表面の傾きに従わない平らな線 —— 自己相関のアーチファクト(層の光学厚の深さに水平に出る)—— は jump_penalty で不利になる。行は明るさの放物線で行以下に詰める。
厚さ = 隣の界面の行の差 × um_per_px(縦の 1 画素の µm、光路)÷ n_group(層ごと、None なら光学厚)。reference_rows を渡したときは、それを 1 本目の界面として行の表に含める。
返り値: dict(kind="oct_boundaries", rows (n, W)、thickness_um (n − 1, W)、mean_um、std_um、contrast)。画像が平ら・範囲が狭すぎて道が置けないときは拒否。
optics के सभी ops गणना से पहले input की जाँच करते हैं (कुछ भी चुपचाप पास नहीं होता):
• Unit argument के नाम में ही लिखी होती है — _mm / _um / _deg / _mrad। mm और µm की गड़बड़ी crash नहीं करती, बल्कि "सही दिखने वाला पर गलत उत्तर" देती है, इसलिए नाम से इसे रोका जाता है। मान के आकार से unit का अनुमान बिल्कुल नहीं लगाया जाता।
• String से ValueError — float('50') सफल हो जाता है, इसलिए बिना parse किया config value length के रूप में निकल जाता (मापा गया: thin_lens('50', '200') सही दिखने वाला 66.667 mm लौटा रहा था)। bool को भी True == 1 वाले implicit promotion के रूप में अस्वीकार किया जाता है।
• complex / masked array से ValueError (केवल real-valued slots; imaginary part को चुपचाप छोड़ना या mask हटाना अस्वीकार)। हर input पर NaN/Inf से ValueError।
• शून्य से भाग और उससे जुड़े मामलों को नाम लेकर अस्वीकार किया जाता है: focal length 0, radius of curvature 0, refractive index <= 0, पूरी तरह opaque aperture (सब 0, इसलिए normalisation 0/0), ऐसा PSF जिसका योग <= 0 हो, S0 = 0 वाला Stokes vector, और front focal point पर रखा object (image अनंत पर)।
• केवल 2 ops non-finite value लौटाते हैं, और वह भी contract के रूप में लिखा हुआ: depth_of_field hyperfocal distance से आगे far_mm = inf (hyperfocal distance की परिभाषा यही है) और gaussian_beam waist पर wavefront_radius_mm = inf (plane wavefront की radius of curvature)। दोनों साथ में एक finite साथी (far_is_infinite / curvature_per_mm) भी लौटाते हैं। इसके अलावा कोई भी चुपचाप आया NaN/Inf अंदर ही पकड़ा जाता है और ValueError देता है — "float64 overflow हुआ" और "उत्तर अनंत है" अलग-अलग दावे हैं, इसलिए पहले को दूसरे का चेहरा पहनाकर नहीं लौटाया जाता।
• Size की सीमा: generate होने वाले grid के लिए optics.MAX_GRID (4096), दिए गए field/PSF/aperture के लिए optics.MAX_FIELD_ELEMENTS (2^24), ABCD element chain के लिए optics.MAX_SYSTEM_ELEMENTS (1024), Zernike के लिए MAX_ZERNIKE_TERMS (512) / MAX_ZERNIKE_ORDER (40) / MAX_ZERNIKE_BASIS (2^25)। छोटे argument से बहुत बड़ा internal allocation होने वाले रास्तों (मापा गया: n_max=40 × 4096² पर 108 GB) को fail-closed ढंग से बंद किया गया है।
• भौतिक रूप से असंभव states भी अस्वीकार: degree of polarisation > 1 वाला Stokes vector, ऋणात्मक transmittance, ऋणात्मक intensity, और n-|m| odd जैसे अमान्य Zernike indices।
• Sample data catalog (DL URL / license) — 2-D के लिए skimage.data (BSD/public) + synthetic, 3-D के लिए real data sources (Stanford/PDS आदि) के DL URL।
• Operators का स्रोत और संदर्भ — इस op family के मूल research/methods के स्रोत।
• Algorithm का मानक स्रोत (लेखक, वर्ष) और उपयोग ऊपर की family usage guide में दिए गए हैं।
• poc_industrial_oct — py -3.11 examples/poc_industrial_oct.py
table को input के रूप में लेते हैं)abcd_matrix · wavefront_stats · sfr_from_edge · paraxial_trace · seidel_coefficients · spot_stats · tolerance_analysis · wavefront_from_opd
oct)oct_design · oct_dispersion_broadening · oct_interferogram · oct_bscan_simulate · oct_from_spectra · oct_k_linearize · oct_dispersion_estimate · oct_dispersion_compensate
*Provenance: octscan.py — OPTICS operator registry. यह per-op note tools/opdocs.py md से अपने-आप generate होता है (हाथ से edit न करें)।*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.