energy_radii — OPTICS imaging op

• Data kinds: image2d → table

• Call: import fullseye as fs; fs.ledger.energy_radii(psf, pixel_pitch_um=1.0, fractions=(0.5, 0.8, 0.9), kind='circle', centre='centroid') (implementation को सीधे call करना हो तो import optics; optics.energy_radii(psf, pixel_pitch_um=1.0, fractions=(0.5, 0.8, 0.9), kind='circle', centre='centroid'), registry से लेना हो तो opsoptics.get("energy_radii"))

उपयोग

वह त्रिज्या जिसमें point-spread function (PSF) की ऊर्जा का दिया गया अनुपात आता है (table)।

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

Solves the encircled (kind='circle') or ensquared (kind='square') energy curve for each fraction in *fractions* (strictly between 0 and 1, at most 16, ascending or not) by bisection on the exact edge-weighted curve of encircled_energy / ensquared_energy: the usual "radius at 50 / 80 / 90 %" figures of merit, in micrometres.

Returns {"kind", "fractions", "radius_um", "diameter_um", "centre_px"}: for a circle the radius, for a square the half-width (diameter_um is twice that, the pixel side); centre_px is (row, col). A fraction that the array never reaches (the window cuts the wings off) raises rather than returning the array edge.

Ground truth: the Airy pattern's radii at 50 / 80 / 90 % are the roots of 1 - J0(v)^2 - J1(v)^2 = fraction with v = 2 pi r / (lambda N) (tested against an independent root find).

Raises ValueError: bad *psf* / *pixel_pitch_um* / *centre* (see encircled_energy); *kind* not 'circle' or 'square'; a fraction outside (0, 1) or non-finite; a fraction the array's energy does not reach inside the largest circle / square that fits.

पूरी 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

पृष्ठभूमि guides (इस op के पीछे की physics और conventions)

• measurement_uncertainty — Measurement uncertainty और calibration — 'हम माप रहे हैं' कहने के लिए क्या चाहिए

• mv_cameras — Industrial camera निर्माता (sensor mapping, line scan / TDI)

• virtual_machine_vision — Virtual machine vision — parameters की सूची और एक object model

संदर्भ (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 करते हैं)

• optical_design_metrics — py -3.11 examples/optical_design_metrics.py

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

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

इसी category के (imaging)

psf_to_mtf · mtf_diffraction · encircled_energy · ensquared_energy · mtf_tangential_sagittal · wavefront_stats · edge_spread · sfr_from_edge


*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.