restoration op• 데이터 종류: measurement → voxel
• 호출: import fullseye as fs; fs.ledger.vol_gaussian_psf(sigma, truncate=4.0)(구현을 직접 호출하려면 import volrestore; volrestore.vol_gaussian_psf(sigma, truncate=4.0), 원장에서 가져오려면 ops3d.get("vol_gaussian_psf"))
총합 1 로 정규화한 3-D 가우스 PSF 커널. *sigma* 는 스칼라 또는 축별.
> 아래 상세 설명은 원문입니다 —— 요약과 제목은 번역되어 있습니다.
`(sz, sy, sx) in voxels; the kernel spans +-truncate*sigma` per axis
(odd size, centre at the middle voxel). The convenient companion to
:func:vol_richardson_lucy when the instrument PSF is well approximated
as Gaussian. Kernels above `PSF_MAX_ELEMENTS` (~256^3) raise
`ValueError` — a PSF that large is an input mistake, not an instrument.
形: 各軸の半径 `r = ceil(truncate * sigma)、辺長 2r + 1`(奇数)の
`(2rz+1, 2ry+1, 2rx+1) float64 配列。軸順は (z, y, x) で、sigma` を
`(sz, sy, sx)` で渡すとその順に対応する(異方 PSF)。3 本の 1-D Gaussian
`exp(-x^2 / (2 sigma^2))` の外積を全体和で割るので、合計はちょうど 1、中心
voxel が最大値。単位は voxel(mm の PSF 幅は spacing で割ってから渡す)。
引数と検証(`ValueError`):
• `sigma`: 正の有限スカラ、または長さ 3 の正の有限値。2 乗が 0 に
アンダーフローするほど小さい値(~1.5e-154 未満)は 0/0 になるので拒否。
• `truncate`: 正の有限値(既定 4.0 = 片側 4σ で打ち切り)。
• 要素数 `prod(2r+1) が PSF_MAX_ELEMENTS(2**24`)を超えると拒否
(σ≈300 で 64 GB を確保しかけた実測が動機)。
使いどころ: `vol_richardson_lucy(vol, psf) の psf`。RL 側でも合計 1 に
再正規化するので、ここで作ったカーネルをそのまま渡してよい。
• 샘플 데이터 카탈로그(DL URL / 라이선스) —— 2-D 는 skimage.data(BSD/public)+ 합성, 3-D 는 실데이터 소스(Stanford/PDS 등)의 DL URL.
• 연산자의 내력·참고문헌 —— 이 연산자 족의 바탕이 된 연구/기법의 출처.
• deconv_fft_restore — py -3.11 examples_3d/deconv_fft_restore.py
voxel 를 입력으로 받는 것)voxel_to_mips · voxel_to_mesh · signed_distance_field · to_points · sobel3d · hessian3d · curvature_maps · edt_jfa
restoration)*Provenance: volrestore.py — 3D 연산자 레지스트리. 이 op 노트는 tools/opdocs.py md 가 자동 생성합니다(직접 편집하지 마세요).*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.