fullseye

ring_artifact_apply — TOMOGRAPHY artifact op

使い方

Give the detector a per-bin gain error — the source of ring artefacts.

A detector bin whose gain is g reports I = g I_true, so after the logarithm the line integral picks up a constant offset -ln g at that bin, the same at every angle. Back-projecting a constant column smears it around the rotation axis, and the reconstruction grows a ring at the radius that bin’s rays are tangent to. One bad pixel, one perfect circle.

The offsets are drawn once from N(0, gain_sigma) with a fixed seed and applied to every row, because the whole point is that the error does not vary with angle — that is what distinguishes a ring from noise, and what makes :func:ring_artifact_remove possible.

:param sinogram: (n_angles, n_detectors). :param gain_sigma: standard deviation of the per-bin offset, >= 0. :param seed: RNG seed; there is no None (determinism is a contract here). :param offsets: explicit (n_detectors,) offsets; overrides the random draw. :returns: (n_angles, n_detectors) float64 sinogram. :raises ValueError: on a negative sigma, a non-int seed, or an offsets whose length is not the detector count.

参考(サンプルデータ・文献)

実行できる例(この op を実際に呼ぶ検証済みサンプル)

型が繋がる次の op(sinogram を入力に取れる)

backproject_sinogram · filtered_backprojection · sart_reconstruct · beam_hardening_apply · beam_hardening_correct · ring_artifact_remove · metal_trace_interpolate · sinogram_center_of_rotation

同カテゴリ(artifact)

beam_hardening_apply · beam_hardening_correct · ring_artifact_remove · metal_trace_interpolate


Provenance: tomography.py — TOMOGRAPHY operator registry. この per-op ノートは tools/opdocs.py md が自動生成(手編集しない)。

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