平行ビームの順投影(サイノグラム)と、フィルタ補正逆投影による再構成、リングアーチファクトやビームハードニングの付与と補正を行います。再構成した体積はそのまま等値面としてメッシュ化できます。
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 → ring_artifact_remove → beam_hardening_correct → fbp_volume → vol_rle_components → vol_rle_volume
radon_transform で順投影(真値つきの検証用)→ ring_artifact_remove / beam_hardening_correct でサイノグラムを補正 → fbp_volume で再構成 → vol_rle_components → vol_rle_volume で領域の体積。
等値面は marching_cubes、3-D の順投影は radon_volume、回転中心のずれは sinogram_center_of_rotation。
iso-value を 1 段変えるだけでボイド体積・肉厚の合否が反転する(poc_ct_void_morphology)。投影数を減らすとフィルタ無し逆投影と 24 本で並び 12 本で逆転(poc_ct_fidelity)。
体積はボクセルピッチの 3 乗。ピッチは既知寸法のファントムを再構成して実測する(公称ピッチで 1 % ずれると体積 3 %)。
radon_transform / fbp_volume / ring_artifact_remove / marching_cubespoc_ct_fidelity、poc_ct_void_morphology