polarization op• Data kinds: कोई नहीं → matrix (केवल arguments से तय होने वाला op — यह कोई image या data input नहीं लेता)
• Call: import fullseye as fs; fs.ledger.mueller_partial_polarizer(angle_deg=0.0, extinction_ratio=1000.0, transmission=1.0) (implementation को सीधे call करना हो तो import optics; optics.mueller_partial_polarizer(angle_deg=0.0, extinction_ratio=1000.0, transmission=1.0), registry से लेना हो तो opsoptics.get("mueller_partial_polarizer"))
• GPU: नहीं (केवल CPU reference implementation)
सीमित extinction ratio वाले linear polariser (linear diattenuator) का Mueller matrix। → (4, 4), angle field हो तो (*shape, 4, 4)।
एक वास्तविक polariser अपनी axis के साथ polarised प्रकाश के लिए intensity transmittance T1 = transmission और उसके लंबवत polarised प्रकाश के लिए T2 = T1 / extinction_ratio देता है। P = T1 + T2, D = T1 - T2, Q = 2 sqrt(T1 T2) और c = cos 2a, s = sin 2a (axis *angle_deg* पर) लिखने पर matrix है
0.5 * [[P, D c, D s, 0],
[D c, P c^2 + Q s^2, (P - Q) c s, 0],
[D s, (P - Q) c s, P s^2 + Q c^2, 0],
[0, 0, 0, Q]]
— textbook diattenuator का block-diagonal रूप, *angle_deg* तक घुमाया हुआ। यह Jones matrix R(a) diag(sqrt T1, sqrt T2) R(-a) के Mueller matrix के भी ठीक बराबर है (शुद्ध diattenuator में न retardance होता है न depolarisation), और tests इस समानता को सामान्य Jones → Mueller रूपांतरण (एक अलग code path) से जाँचते हैं।
सत्य मान (closed form, tests में जाँचे गए): unpolarised प्रकाश intensity (T1 + T2)/2 और degree of polarisation (ER - 1)/(ER + 1) के साथ निकलता है; ऐसे दो polarisers को parallel / crossed रखने पर transmission (T1^2 + T2^2)/2 / T1 T2 है — जोड़ी से मापा गया extinction ratio ER नहीं बल्कि (ER^2 + 1)/(2 ER) है; axis से theta कोण पर पूर्णतः polarised input के लिए Malus का नियम T2 + (T1 - T2) cos^2(theta) बन जाता है; ER -> अनंत पर यह T1 * mueller_element ("polarizer") की ओर जाता है; ER = 1 एक neutral-density filter T1 * I है।
*angle_deg* एक array भी हो सकता है (हर entry के लिए एक matrix, mueller_element की तरह); *extinction_ratio* और *transmission* scalars हैं। कोण module की परिपाटी मानते हैं (+x से +y = rows में नीचे की ओर, दिखाई गई image में clockwise)।
Raises ValueError: non-finite / non-real / bool / string arguments; *extinction_ratio* < 1 (1 से कम "extinction ratio" का अर्थ है कि axis दूसरी है — reciprocal देने के बजाय 90 degree घुमाएँ) या अनंत (ideal polariser के लिए mueller_element); *transmission* (0, 1] से बाहर (passive element amplify नहीं कर सकता, और zero transmission वाले element की कोई axis नहीं); MAX_FIELD_ELEMENTS / 16 से अधिक entries वाला angle field।
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 में दिए गए हैं।
• polarization_camera_pipeline — py -3.11 examples/polarization_camera_pipeline.py
matrix को input के रूप में लेते हैं)abcd_trace · emva_filter_distortion · mueller_apply · mueller_checks · equirect_rotate · pixel_shift_capture · pixel_shift_merge · pixel_shift_motion_mask
polarization)jones_element · jones_apply · stokes_from_jones · mueller_element · mueller_apply · stokes_analyze · polarization_demosaic · polarization_demosaic_color
*Provenance: optics.py — OPTICS operator registry. यह per-op note tools/opdocs.py md से अपने-आप generate होता है (हाथ से edit न करें)।*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.