fly_lamina_filter — FLYVISION lamina op

• Data kinds: matrix → matrix

• Call: import fullseye as fs; fs.ledger.fly_lamina_filter(movie, dt_s, tau_adapt_s=0.2, tau_lp_s=0.02, mode='divisive', floor=0.001) (implementation को सीधे call करना हो तो import flyvision; flyvision.fly_lamina_filter(movie, dt_s, tau_adapt_s=0.2, tau_lp_s=0.02, mode='divisive', floor=0.001), registry से लेना हो तो opsflyvision.get("fly_lamina_filter"))

उपयोग

photoreceptor adaptation + lamina का band-pass: intensities इनपुट, contrast आउटपुट।

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

The first thing the optic lobe does to a picture is throw away its brightness. A photoreceptor adapts to the running mean light level and the large monopolar cells (L1/L2) report the *deviation* from it, so the same scene at dawn and at noon arrives at the motion detectors as the same signal. Two closed-form stages, in that order:

1. adaptation — a first-order low-pass of time constant *tau_adapt_s* per ommatidium is the adaptation state a(t). mode="divisive" returns the Weber contrast (x - a)/(a + eps) (eps = floor * mean(x), so it scales with the picture and a dark ommatidium cannot divide by zero); mode="subtractive" returns x - a, which is the same high-pass without the gain control.

2. membrane — a first-order low-pass of time constant *tau_lp_s*, the cell's own bandwidth.

movie: (T, n) intensities, rows = time. mode="divisive" refuses a negative entry (a negative light level is not a measurement) and an all-zero movie (its contrast is 0/0, which would be fabricated rather than measured).

dt_s: sample interval, seconds. tau_adapt_s / tau_lp_s: the two time constants, seconds. floor: the divisive guard, relative to the mean intensity.

Returns (T, n) float64 contrast.

Ground truth, both exact rather than approximate:

• Weber invariance. In "divisive" mode, scaling the whole movie by any positive constant returns *the same array* — both a and eps scale with it. That is the point of the stage and the tests pin it to machine precision.

• The transfer is the product of the two first-order filters. With A = 1 - exp(-dt/tau_adapt) and B = 1 - exp(-dt/tau_lp), the steady-state gain at angular frequency w is |1 - H_A(w)| * |H_B(w)| where H(w) = C/(1 - (1-C) exp(-i w dt)) — a band-pass that blocks DC exactly and is measured at four frequencies in the tests.

Raises ValueError: a non-2-D / too-short / non-finite *movie*, a movie

over MAX_MOVIE_ELEMENTS, a non-positive *dt_s* / *tau_adapt_s* / *tau_lp_s*, a non-positive *floor*, an unknown *mode*, and (divisive only) a negative or all-zero movie.

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

• fly_vision 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_fly_optomotor_steering — py -3.11 examples/poc_fly_optomotor_steering.py

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

fly_onoff_split · fly_t4t5_field · fly_flow_from_directions · fly_egomotion_from_flow · fly_hs_readout

इसी category के (lamina)

fly_onoff_split


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

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