oct op• Data kinds: कोई नहीं → table (केवल arguments से तय होने वाला op — यह कोई image या data input नहीं लेता)
• Call: import fullseye as fs; fs.ledger.oct_interferogram(sample=None, center_nm=840.0, bandwidth_nm=50.0, mode='sd', n_pixels=2048, wl_min_nm=None, wl_max_nm=None, spectrometer_fwhm_px=1.0, linewidth_nm=None, sweep_nonlinearity=0.0, source_spectrum=None, photons_ref_per_pixel=20000.0, photons_mirror_per_pixel=None, eta=0.8, sigma_d=10.0, K=1.0, full_well=None, quantize=False, reference_phase_rad=0.0, mismatch_material=None, mismatch_length_mm=0.0, focus_um=None, beam_waist_um=None, oversample=8, noise=True, seed=0, grid='lambda') (implementation को सीधे call करना हो तो import octscan; octscan.oct_interferogram(sample=None, center_nm=840.0, bandwidth_nm=50.0, mode='sd', n_pixels=2048, wl_min_nm=None, wl_max_nm=None, spectrometer_fwhm_px=1.0, linewidth_nm=None, sweep_nonlinearity=0.0, source_spectrum=None, photons_ref_per_pixel=20000.0, photons_mirror_per_pixel=None, eta=0.8, sigma_d=10.0, K=1.0, full_well=None, quantize=False, reference_phase_rad=0.0, mismatch_material=None, mismatch_length_mm=0.0, focus_um=None, beam_waist_um=None, oversample=8, noise=True, seed=0, grid='lambda'), registry से लेना हो तो opsoptics.get("oct_interferogram"))
एक layer stack (sample) के एक OCT measurement का spectral interferogram (table)।
> नीचे का विस्तृत विवरण मूल text में है — सारांश और headings का अनुवाद हो चुका है।
sample = {"ambient": "air", "surface_um": 表面の光路(参照腕との差、空気中 µm、負なら鏡像の側),
"layers": [{"material", "thickness_um", "attenuation_per_mm"(強度の減衰係数 µ_t、往復で exp(−2µd)), "scatterers_per_mm"(層の中の点散乱体の密度)、"scatter_amplitude"(1 個の振幅反射率)}, …],
"substrate": 材料}。省略時は空気 / PMMA 25 µm / ポリカーボネート 40 µm / ガラス相当の基板。
場の和: E_s(k) = Σ_j r_j(k) T_j A_j exp(2ik Δ_j(k))、Δ_j(k) = surface + Σ_(上の層) n_i(k) d_i (位相屈折率で位相、群遅延は自然に出る)、r_j = (n_上 − n_下)/(n_上 + n_下)(各 k で)、T_j = Π(1 − r²)、A_j = exp(−Σ µ_i d_i)。参照光(振幅 1、位相 reference_phase_rad)と足して |E_r + s·E_s|² ×光源 S(k) (s² = photons_mirror/photons_ref)—— 相互相関・自己相関・DC が全部入る。
検出: mode="sd" は分光器(grid="lambda" = λ 等間隔の画素、"k" = k 等間隔の画素。画素の箱で積分 + 分解能 spectrometer_fwhm_px のガウス)、mode="ss" は時間等間隔の波長掃引(k(t) = k_lo + Δk (t + ε t(1−t))、ε = sweep_nonlinearity、瞬時線幅 linewidth_nm のガウス、各標本は瞬時の値 —— 画素の箱が無いので roll-off は線幅だけで決まる)。細かい等間隔 k(1 画素 oversample 点)で作ってから、分光器は画素の区間で積分、掃引は標本の k で読む。
雑音: ショット(電子の平方根)+ 読み出し σ_d e⁻、×K で DN、full_well で飽和、quantize で整数化。
分散の不一致: 試料腕だけに mismatch_material が mismatch_length_mm(群遅延は参照腕で合わせる)。
焦点: focus_um(物理深さ)と beam_waist_um を渡すと、共焦点の結合 (1 + ((z − z_f)/z_R)²)^(−1/2)。
返り値: dict(kind="oct_interferogram", spectrum (N,) DN, wavelength_nm, k_rad_per_um(画素の中心の真の k)、reference_spectrum(参照腕だけ・雑音なし)、sample_only_spectrum(試料腕だけ・雑音なし = 自己相関の素)、interfaces(真値: 物理深さ・群遅延の光路・反射率)、settings)。
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_bscan_simulate · oct_from_spectra · oct_k_linearize · oct_dispersion_estimate · oct_dispersion_compensate · oct_ascan
*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.