typed op• Data kinds: counts → counts
• Call: fullseye.apply(img, "tb_tcspc_coates_correct", a=0.5, b=0.5) (2-D model = एक image + दो scalar knobs a,b∈[0,1])
*Figure synthetic 128×128 input पर सच में चलाकर मिला output है। बाईं ओर input, दाईं ओर output। Point clouds ऊपर से देखे गए scatter के रूप में (brightness = z), 1-D series line के रूप में, volumes z दिशा में maximum-intensity projection के रूप में, videos बीच के frame के रूप में, complex images amplitude के रूप में, और जो return values चित्र नहीं बन सकते वे स्वयं values के रूप में दिखाए जाते हैं।*
*Knob a output नहीं बदलता (मापा गया: 0.1 / 0.5 / 0.9 पर एक जैसा)।*
*Knob b output नहीं बदलता (मापा गया: 0.1 / 0.5 / 0.9 पर एक जैसा)।*
Stages (पहले आने वाले ops → यह op, बाएँ से क्रम में):
▸ tb_tcspc_coates_correct: stages का figure (docs site)
दूसरी images पर भी (synthetic scene / photo / coins। ऊपर की row में inputs, नीचे की row में उनके outputs। Knobs default पर):
▸ tb_tcspc_coates_correct: दूसरी images पर output (docs site)
TCSPC pile-up को exact रूप से उलटता है (Coates's estimator) — early-photon bias।
> नीचे का विस्तृत विवरण मूल text में है — सारांश और headings का अनुवाद हो चुका है।
Classical TCSPC records at most one photon per excitation cycle: the first one. Late bins are therefore starved, because the cycles in which an early photon arrived never reach them, and the measured histogram is biased toward short arrival times — a dToF depth read straight off a piled-up histogram is *too close*, and a fluorescence lifetime is *too short*.
Coates's estimator inverts that exactly. With N_k the measured counts in bin k and C the number of excitation cycles, the number of cycles that survived to reach bin k is D_k = C - sum_{j<k} N_j, the per-cycle detection probability in that bin is p_k = N_k / D_k and the pile-up-free per-cycle intensity is lambda_k = -ln(1 - p_k). This op returns C * lambda_k — the histogram the same scene would have produced if the detector could record every photon — so it is directly comparable to the measured one.
This is an exact inverse, not a linearisation: build a histogram from a known lambda through the forward model N_k = C * exp(-sum_{j<k} lambda_j) * (1 - exp(-lambda_k)) and Coates returns lambda to machine precision (measured max relative error 1.6e-15 in the tests, on a pile-up so severe that the last bin was suppressed to 14.8% of its true counts).
*hist* is the 1-D measured histogram (counts per bin); *cycles* the number of excitation cycles (laser pulses) that produced it.
Raises ValueError: negative, non-finite or non-1-D *hist*, a non-positive or non-integer *cycles*, a histogram whose total exceeds *cycles* (impossible: at most one photon per cycle — a sure sign that *cycles* is wrong or the data are not first-photon TCSPC), and any bin that consumed every remaining cycle (p_k = 1, where -ln(0) is inf).
Typed bridge of the photon op tcspc_coates_correct into the 2-D evolution registry: the same implementation, called under the op(v, a, b) convention. This op has no tunable parameter; a and b are unused.
• 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 के स्रोत।
नीचे का program सच में चलता है, यह जाँचा गया है (figure वाला ही input)। Studio की help में यह block buttons बन जाता है, जिनसे इसे वहीं load करके चलाया जा सकता है।
img_to_counts 0.50 0.50 tb_tcspc_coates_correct 0.50 0.50
▸ यह pipeline load करें · Load करके चलाएँ
नीचे के उदाहरण मूल ledger op tcspc_coates_correct को call करते हैं। यह bridge op वही implementation है, बस fn(v, a, b) convention के अनुसार ढाला गया है, इसलिए व्यवहार वैसा ही लागू होता है (केवल call का रूप अलग है)।
• photon_timeresolved — py -3.11 examples/photon_timeresolved.py
• poc_dtof_ranging — py -3.11 examples/poc_dtof_ranging.py
counts को input के रूप में लेते हैं)identity · tb_spad_deadtime_apply · tb_spad_deadtime_correct · tb_tcspc_irf_convolve · tb_tcspc_background_subtract · tb_dtof_depth · tb_countrate_to_counts · tb_counts_to_countrate
typed)tb_points_to_voxel · tb_estimate_point_normals · tb_iss_keypoints · tb_project_points · tb_render_point_depth · tb_statistical_outlier_removal · tb_radius_outlier_removal · tb_voxel_grid_downsample
*Provenance: ops.py — 2D operator registry. यह per-op note tools/opdocs.py md से अपने-आप generate होता है (हाथ से edit न करें)।*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.