sensorchar op• Data kinds: image2d → table
• Call: import fullseye as fs; fs.ledger.emva_spectrogram(image, axis='horizontal', band=None) (implementation को सीधे call करना हो तो import sensorchar; sensorchar.emva_spectrogram(image, axis='horizontal', band=None), registry से लेना हो तो opsoptics.get("emva_spectrogram"))
क्षैतिज और ऊर्ध्वाधर spectrogram (खंड 8.6, समीकरण 70-72)।
> नीचे का विस्तृत विवरण मूल text में है — सारांश और headings का अनुवाद हो चुका है।
各行(axis="vertical" なら各列)から全体の平均を引き、ŷ[m][v] = (1/√N) Σ_n y[m][n] exp(−2πinv/N)(式 70)、P[v] = (1/M) Σ_m ŷ ŷ*(式 71)。白い雑音なら平らで、その高さが空間分散 s²(Parseval: P の全 v の平均 = 画像の分散)。*band* = (f_lo, f_hi) [cycles/pixel] を渡すと、その帯の周期模様の振幅 a = 2 ((1/N) Σ_{v=vmin}^{vmax} P[v])^{1/2} (式 72、片側の和)を返す。式 72 は帯の中の白い雑音の力も足す —— 照合用(規格の数ではない)に、P の中央値を白い床とみて帯の点の数 × 床を引いた amplitude_above_floor も返す。★8.6 節は高域フィルタ(8.1 節)をかけた DSNU 画像 ⟨y_dark⟩ と PRNU 画像 ⟨y_50⟩ − ⟨y_dark⟩ から求めると書く。
返り値 {"freq": v/N (v = 0..N/2), "power": P[v], "sqrt_power", "mean_power", "variance", "white_floor", "amplitude", "amplitude_above_floor", "band", "axis"} (mean_power は全 v の P の平均、variance は平均を引いた画像の分散 —— 両者は等しい)。
Raises ValueError: 2-D の実数でない / NaN・inf / axis が違う / 変換する軸が 4 未満 / band が [0, 0.5] の中の lo ≤ hi でない / 帯に周波数の点が無い。
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_sensor — py -3.11 examples/poc_emva1288_sensor.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.