vol_richardson_lucy — 3D `restoration` op

• Datenarten: voxel × voxel → voxel

• Aufruf: import fullseye as fs; fs.ledger.vol_richardson_lucy(vol, psf, iterations=10, clip_tiny=1e-12) (die Implementierung direkt: import volrestore; volrestore.vol_richardson_lucy(vol, psf, iterations=10, clip_tiny=1e-12); aus dem Register: ops3d.get("vol_richardson_lucy"))

Verwendung

Richardson–Lucy-Entfaltung eines nichtnegativen Volumens mit bekannter PSF.

> Die ausführliche Beschreibung unten ist der Originaltext — Zusammenfassung und Überschriften sind übersetzt.

*psf* is a 3-D non-negative kernel (any odd/even size smaller than the

volume; it is normalised to sum 1 internally so overall intensity is

preserved). *iterations* trades sharpness against noise amplification —

5-30 is the practical range (see the module notes on semi-convergence).

Negative voxels are refused (RL is a Poisson model) — except *rounding

dust*: values no lower than `-NEGATIVE_DUST_TOL * max|vol|` (1e-9

relative; an FFT-blurred observation typically carries -1e-16) are clipped

to 0 instead of rejected, so the module's own forward model feeds back in.

Returns the deblurred `(D, H, W)` float64 volume (non-negative).

Measured on the test scene (binary sphere pair blurred by a sigma-2

Gaussian): the RMSE to ground truth falls to 0.81x the blurred

observation's at 10 iterations and 0.68x at 50 — genuine but *gradual*,

because the residual is dominated by the spheres' hard edges, which RL

recovers slowly. What converges fast is the *forward consistency*:

re-blurring the estimate reproduces the observation almost exactly (that

is the quantity the RL update actually optimises).

Referenzen (Beispieldaten, Literatur)

• Katalog der Beispieldaten (Download-URLs / Lizenzen) — 2-D nutzt skimage.data (BSD/Public Domain) plus synthetische Bilder, 3-D nennt Download-URLs echter Datenquellen (Stanford, PDS, …).

• Herkunft und Literatur der Operatoren — die Quellen der Forschung/Verfahren, auf denen diese Operatorfamilie beruht.

Ausführbare Beispiele (verifizierte Samples, die diesen Operator wirklich aufrufen)

• deconv_fft_restore — py -3.11 examples_3d/deconv_fft_restore.py

Typkompatible Folge-Operatoren (nehmen `voxel` als Eingabe)

voxel_to_mips · voxel_to_mesh · signed_distance_field · to_points · sobel3d · hessian3d · curvature_maps · edt_jfa

Gleiche Kategorie (`restoration`)

vol_gaussian_psf

*Provenance: volrestore.py — 3D Operator-Registry. Diese Notiz wird von tools/opdocs.py md erzeugt (nicht von Hand bearbeiten).*

vol_richardson_lucy — 3D restoration op