sensorchar op• Data kinds: signal × signal → table
• Call: import fullseye as fs; fs.ledger.emva_bspline_regression(x, y, n_intervals=8, x_max=None, x_eval=None) (implementation को सीधे call करना हो तो import sensorchar; sensorchar.emva_bspline_regression(x, y, n_intervals=8, x_max=None, x_eval=None), registry से लेना हो तो opsoptics.get("emva_bspline_regression"))
least-squares cubic B-spline fit और उसकी ढलान (General खंड 6.6, समीकरण 44-49)।
> नीचे का विस्तृत विवरण मूल text में है — सारांश और headings का अनुवाद हो चुका है।
区間 [0, x_max] を P 等分(Δ = x_max/P)、基底 f_p(x) = β₃(x/Δ − (p − 1))(p = 0..P+2 の P + 3 個、式 44・46)。係数 a_p は Σ_q (Σ_p a_p f_p(x_q) − y_q)² の最小(式 47、Q > P + 3 点が要る)。傾きは Σ_p a_p β₃'(x/Δ − (p − 1)) / Δ(式 48 に連鎖律の 1/Δ を足した —— 本文は正規化した変数での傾き)。General は「露光の段数が推奨どおりなら 6〜8 区間」と書く。
*x_max*: 省略すると max(x)。*x_eval*: 当てはめを評価する追加の点(区間 [0, x_max] の中だけ —— 外挿しない)。返り値 {"coef", "n_intervals", "x_max", "delta", "fit", "slope", "residual_rms", "x_eval", "fit_eval", "slope_eval"}。
Raises ValueError: 1-D でない・長さが違う / NaN・inf / P が 1 以上の整数でない / 点が P + 3 以下 / x が [0, x_max] の外 / 区間に点が無く係数が決まらない(計画行列の階数が P + 3 未満)/ x_eval が区間の外。
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_emva1288_in_camera_processing — py -3.11 examples/poc_emva1288_in_camera_processing.py
table को input के रूप में लेते हैं)abcd_matrix · wavefront_stats · sfr_from_edge · paraxial_trace · seidel_coefficients · spot_stats · tolerance_analysis · wavefront_from_opd
sensorchar)emva_pair_statistics · emva_photon_transfer · emva_quantum_efficiency · emva_linearity_error · emva_snr_curve · emva_sensitivity_threshold · emva_dynamic_range · emva_dark_current
*Provenance: sensorchar.py — OPTICS operator registry. यह per-op note tools/opdocs.py md से अपने-आप generate होता है (हाथ से edit न करें)।*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.