grating_interf op• Data kinds: image2d → images
• Call: import fullseye as fs; fs.ledger.grating_stepping_simulate(transmission=None, dpc_rad=None, darkfield=None, n_steps=8, ref_visibility=0.25, ref_phase_rad=0.0, steps_per_period=None, flux=None, seed=0, shape=(64, 64)) (implementation को सीधे call करना हो तो import gratinginterf; gratinginterf.grating_stepping_simulate(transmission=None, dpc_rad=None, darkfield=None, n_steps=8, ref_visibility=0.25, ref_phase_rad=0.0, steps_per_period=None, flux=None, seed=0, shape=(64, 64)), registry से लेना हो तो opsoptics.get("grating_stepping_simulate"))
• Ledger के ज़रिए return value: fullseye.ledger.grating_stepping_simulate(...) केवल घोषित out type images का value लौटाता है (मूल function सहायक जानकारी भी लौटाता है)। जो हिस्सा छोड़ दिया गया है उसकी ज़रूरत हो तो fullseye.ledger.grating_stepping_simulate.raw(...) का उपयोग करें, या gratinginterf.grating_stepping_simulate को सीधे call करें।
तीन maps (transmission, differential phase, dark field) से phase-stepping image stack (N, H, W) बनाता है (forward model)।
> नीचे का विस्तृत विवरण मूल text में है — सारांश और headings का अनुवाद हो चुका है।
各画素の k 枚目(k = 0..N−1)は I_k = flux·T·(1 + V0·D·cos(φ0 + φ − 2πk/P))。T = transmission(吸収)、φ = dpc_rad(縞の位相のずれ = 微分位相。屈折角 α とは φ = 2π d α/p2)、D = darkfield(可視度の比、1 = 散乱なし)、V0 = ref_visibility(試料なしの縞の可視度)、φ0 = ref_phase_rad(格子のむらによる基準の位相。スカラーか図)、P = steps_per_period (1 周期のステップ数。None は N = 1 周期を N 等分)。3 枚とも None なら試料なし(基準)の画像列になる: T = 1、φ = 0、D = 1。
flux: None なら雑音なしの期待値(flux = 1)。数を渡すと、T = 1 の画素の 1 枚あたりの平均光子数として Poisson 雑音を足した整数の計数を float で返す(seed で再現)。
返り値: (N, H, W) float64。grating_three_image_retrieve に試料ありと基準の 2 組を渡すと 3 枚が戻る。Raises ValueError: 図が 2-D でない / 形がそろわない / NaN・Inf、T < 0、D が [0, 1] の外、V0 が (0, 1] の外、n_steps が 3..256 の外、P ≤ 0、flux ≤ 0。
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 में दिए गए हैं।
• grating_interferometry — py -3.11 examples/grating_interferometry.py
• poc_grating_interferometry — py -3.11 examples/poc_grating_interferometry.py
images को input के रूप में लेते हैं)film_thickness_map · grating_three_image_retrieve · grating_step_period_estimate · render_transport · separate_direct_global · micro_phase_decode
grating_interf)talbot_grating_design · grating_three_image_retrieve · darkfield_scatter_relation · talbot_carpet · talbot_source_visibility · grating_step_period_estimate · dpc_integrate
*Provenance: gratinginterf.py — OPTICS operator registry. यह per-op note tools/opdocs.py md से अपने-आप generate होता है (हाथ से edit न करें)।*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.