spectral op• Data kinds: pairs → pairs
• Call: import fullseye as fs; fs.ledger.kubelka_munk_reflectance(absorption, scattering_per_mm=5.0, thickness_mm=1.0, background_reflectance=0.0, diffuse_factor=2.0) (implementation को सीधे call करना हो तो import spectralrad; spectralrad.kubelka_munk_reflectance(absorption, scattering_per_mm=5.0, thickness_mm=1.0, background_reflectance=0.0, diffuse_factor=2.0), registry से लेना हो तो opsoptics.get("kubelka_munk_reflectance"))
Kubelka–Munk के closed form में scattering layer की diffuse reflectance (absorption coefficient की curve से)।
> नीचे का विस्तृत विवरण मूल text में है — सारांश और headings का अनुवाद हो चुका है।
K = *diffuse_factor* · α(拡散光は斜めに通るので平行光の吸収係数の約 2 倍、Kubelka 1948)、S = *scattering_per_mm*。*thickness_mm*(既定 1 mm)が None なら無限厚 R∞ = 1 + K/S − √((K/S)² + 2K/S)、有限の厚さ X なら背景 R_g の上の層:
R = [1 − R_g (a − b coth(bSX))] / [a − R_g + b coth(bSX)], a = 1 + K/S, b = √(a² − 1)
(Kubelka, J. Opt. Soc. Am. 38, 448 (1948))。極限: X → ∞ で R∞、K → 0 で(SX(1 − R_g) + R_g)/(1 + SX(1 − R_g))、R_g = R∞ なら厚さに依らず R∞。逆は K/S = (1 − R∞)²/(2R∞)。
*absorption* は pairs(λ nm, α 1/mm)(materials.material_nir_absorption の返り値)。*scattering_per_mm* は数か、*absorption* と同じ波長の pairs。
返り値: pairs(λ nm, R)。界面の鏡面反射は含まない(spectral_reflectance_from_material が Saunderson の補正で足す)。
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_hyperspectral_foreign_matter — py -3.11 examples/poc_hyperspectral_foreign_matter.py
pairs को input के रूप में लेते हैं)mtf50 · energy_to_photons · spectral_product · spectral_integral · spectral_stats · hyperspectral_capture · film_thickness_fft · film_thickness_fit
spectral)blackbody_spectrum · cie_illuminant_a · led_spectrum · bandpass_spectrum · silicon_photodiode_qe · energy_to_photons · spectral_product · spectral_integral
*Provenance: spectralrad.py — OPTICS operator registry. यह per-op note tools/opdocs.py md से अपने-आप generate होता है (हाथ से edit न करें)।*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.