separate_direct_global — OPTICS light_transport op

• Data kinds: images → table

• Call: import fullseye as fs; fs.ledger.separate_direct_global(images, fill=0.5, black_level=0.0, scheme='shift') (implementation को सीधे call करना हो तो import lighttransport; lighttransport.separate_direct_global(images, fill=0.5, black_level=0.0, scheme='shift'), registry से लेना हो तो opsoptics.get("separate_direct_global"))

उपयोग

खिसकाए गए उच्च-आवृत्ति पैटर्नों की छवियों से, हर pixel के L_max और L_min द्वारा direct और global (interreflection) घटकों को अलग करना (Nayar आदि 2006, table)।

> नीचे का विस्तृत विवरण मूल text में है — सारांश और headings का अनुवाद हो चुका है।

With a fraction α of the projector lit in a pattern of high frequency, a pixel whose scene point is lit sees L_d + α L_g (the global part only depends on the average power when it is smooth over the pattern period) and an unlit one sees only the global part.

• scheme="shift" — every image is the same pattern (fill α) shifted: L_max = L_d + (α + b(1−α)) L_g, L_min = b L_d + (α + b(1−α)) L_g, so L_d = (L_max − L_min)/(1 − b).

• scheme="complement" — a pattern and its complement (Nayar et al. eq. 8): L_max = L_d + (α + b(1−α)) L_g, L_min = b L_d + (1 − α + αb) L_g (assumes L_max is the lit image).

Both reduce to L_d = L_max − L_min, L_g = 2 L_min for α = 0.5, b = 0. *black_level* b is the relative brightness of an unlit projector pixel.

Returns a dict (a table) with direct, global, l_max, l_min (each (H, W)).

Raises ValueError: fewer than 2 images, non-2-D or non-finite; *fill* outside (0, 1); *black_level* outside [0, 1); unknown *scheme*.

पूरी family का input contract (fail-closed)

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।

उपयोग की विस्तृत guide

• optics_imaging family guide

संदर्भ (sample data, साहित्य)

• 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 में दिए गए हैं।

चलाए जा सकने वाले उदाहरण (verified samples जो सच में यह op call करते हैं)

• poc_interreflection_fringe_error — py -3.11 examples/poc_interreflection_fringe_error.py

वे अगले ops जिनसे type जुड़ता है (table को input के रूप में लेते हैं)

abcd_matrix · wavefront_stats · sfr_from_edge · paraxial_trace · seidel_coefficients · spot_stats · tolerance_analysis · wavefront_from_opd

इसी category के (light_transport)

groove_profile · form_factors_2d · transport_matrix_radiosity · high_frequency_patterns · render_transport · fringe_multipath_phase_error · predict_fringe_phase_error · micro_phase_decode


*Provenance: lighttransport.py — OPTICS operator registry. यह per-op note tools/opdocs.py md से अपने-आप generate होता है (हाथ से edit न करें)।*

© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.