平行ビームの順投影(サイノグラム)と、フィルタ補正逆投影による再構成、リングアーチファクトやビームハードニングの付与と補正を行います。再構成した体積はそのまま等値面としてメッシュ化できます。
Forward parallel-beam projection to a sinogram, filtered back-projection to a slice or a volume, and the injection and correction of ring artefacts and beam hardening. The reconstructed volume can be turned straight into a mesh.
向く: 投影数・雑音・アーチファクトが結果をどこで壊すかを、真値つきで詰めること。
向かない: ★閾値(iso-value)を 1 段変えるだけでボイド体積や肉厚の合否が反転します。1 つの数字で合否を出す前に、poc_ct_void_morphology が測っている「合否 1 個の数字は寿命に効く形に盲目」を読んでください。
import fullseye as fs
import numpy as np
sl = np.zeros((64, 64)); sl[24:40, 24:40] = 1.0
sino = fs.ledger.radon_transform(sl, np.linspace(0, 180, 120, endpoint=False))
print('サイノグラム', np.asarray(sino).shape)
radon_transform / fbp_volume / ring_artifact_remove / marching_cubespoc_ct_fidelity、poc_ct_void_morphology