Trajectory optimisation for cone-beam computed tomography (CT) determines which information sparse-view scans acquire. Fixed candidate pools prevent off-grid refinement and require new object-specific precomputation for each acquisition manifold. We make every source pose an individual continuous variable and move all poses jointly by gradient ascent on the scanner's kinematic manifold. The objective combines soft-Tuy plane coverage, continuous View Covariance Loss, and an analytic attenuation-aware ray-bundle penalty. The same optimiser handles circular, limited C-arm, two-axis, and freesphere parametrisations. On a Defrise flange, continuous selection recovers laminar defects invisible to a circular orbit, matches discrete swap search on the free sphere at the sparser budget, and leads at the denser one, with the same objective evaluated in every arm. A moderate elevation band already recovers most of the free-sphere gain at the defects, so the same optimiser transfers to bounded scanner envelopes. Photon noise preserves the ordering on the flange and compresses it on a dense fuel nozzle. Sparseprescan planning benefits from matching prescan and planned acquisition manifolds. Selection takes seconds rather than minutes without an object-specific reconstruction basis. Prescan-planned poses were executed on a robot CT bench and reconstructed in a common frame, demonstrating feasibility but no consistent metric gain over uniform band sampling. Continuous pose optimisation incorporates attenuation and scanner constraints directly into sparse-view acquisition design.
Figures & tables
Fig. 1: Overview of the gradient-based trajectory-optimisation pipeline. (a) Attenuation model μ . (b) The viewing sphere carries two families of symbols. Red dots are the source positions si , and blue arrows from the centre are the sampled Radon normals nj , one for each plane whose coverage is scored. Thick red arrows are the coverage gradient, which is tangential to the sphere by Proposition 2 . Every source carries such a gradient, and three are drawn so that the panel stays readable. (c) The objective combines the soft-Tuy coverage surrogate, the attenuation-aware bundle penalty, and the continuous VCL information score; Adam updates the kinematic parameters on the manifold. (d) Two outputs from separate runs. Above, red dots are the optimised source positions V⋆ of a coverage-only free-sphere run at k=24 , thin red lines are their rays to the region of interest, and the black cross marks its centre. Below, SART reconstructions at k=80 for the two-axis gantry and the free sphere, named in the image.
Specification
Defrise flange
Fuel nozzle
Digital camera
Material
aluminium, steel
316L stainless steel
plastics, circuit boards, metal
Max. photon energy (keV)
monochromatic
180
175
Detector dimensions (pixels)
256×256
1024×1024
3072×3072
Pixel pitch in detector (mm)
0.5
0.2
0.139
Source-to-detector distance (mm)
900
808.508
1995.6
Source-to-isocentre distance (mm)
500
243.307
997.1
TABLE I: Configuration of the three test objects. The flange is simulated under a monochromatic model, and the camera detector is mean-pooled by a factor of four before reconstruction.
Fig. 2: The Defrise flange phantom. (a) Three-dimensional view of the aluminium body in translucent blue with two press-fit steel inserts in red. (b) Central sagittal slice with the two LoF stacks far off the source plane and the two tilted control cracks in the column. (c) Axial slice through the upper flange disk with the LoF layer, the steel insert, and the pore cluster in the attenuation shadow of the insert.
Fig. 3: Directional absorption of the Defrise flange, shown as the bundle path integral τˉ towards three target points as a function of the source azimuth and elevation on the acquisition sphere. Dark is transparent. Towards the isocentre the transparent region is the equatorial band, since only the column lies in the way there. Towards an off-plane defect the minimum moves into the hemisphere of that defect and leaves a narrow equatorial notch. We aim the term at the isocentre because that choice needs no assumption about which defect matters.
Fig. 4: The scanned object mounted on the rotation stage (a) and the custom robot-based cone-beam CT test bench with the two-axis source and detector manipulators (b).
k=40
k=80
Object
Selector
PSNR
SSIM
HFEN
PSNR
SSIM
HFEN
Defrise flange, 3843
VCLS, discrete
45.82±0.02
0.981±0.000
1.81±0.02
49.76±0.09
0.992±0.000
0.91±0.02
soft-Tuy greedy, discrete
41.32±0.56
0.951±0.003
3.46±0.37
45.25±0.39
0.980±0.000
1.87±0.15
cov+bundle
44.66±0.60
0.963±0.004
2.00±0.15
47.07±0.39
0.981±0.000
1.43±0.10
cov+VCL
43.91±0.53
0.962±0.002
2.36±0.21
49.36±0.15
0.989±0.000
0.95±0.02
cov+VCL+bundle
45.72±0.25
0.970±0.002
1.78±0.07
50.32±0.10
0.988±0.000
0.81±0.00
TABLE II: Sphere-based selectors on both virtual objects, cold-started, mean ± population standard deviation. PSNR in dB and SSIM are higher-is-better, HFEN is lower-is-better. The soft-Tuy greedy row is the discrete fixed-grid selector on the same 720 -point sphere and the initialisation of every continuous row, so the three arms below it measure what leaving the candidate grid adds. Bold marks the best value per column within each object.
Fig. 5: Noise-free reconstructions at k=40 , the nozzle rendered at 2563 for legibility. The top row is a central axial slice of the fuel nozzle with the cooling-hole ring in the inset, the bottom row a central sagittal slice of the Defrise flange with the upper lack-of-fusion stack. The nozzle selections are visually equivalent, as the sub-decibel margins of Table II predict, whereas the flange separates them in the streak background.
k=40
k=80
Reachable set
PSNR
SSIM
HFEN
ROI-PSNR
PSNR
SSIM
HFEN
ROI-PSNR
equatorial circle
42.56±0.00
0.972±0.000
2.91±0.00
25.36±0.00
44.67±0.00
0.985±0.000
2.14±0.00
26.00±0.00
band, ±30∘ elevation
43.93±0.21
0.958±0.002
2.26±0.09
30.61±1.01
48.55±0.06
0.987±0.000
1.11±0.02
36.83±0.84
C-arm, ±110∘/±45∘
42.08±0.51
0.946±0.004
2.88±0.20
30.45±1.34
46.98±0.66
0.980±0.002
1.48±0.17
36.91±1.70
full sphere
45.61±0.34
0.968±0.003
1.78±0.08
34.64±0.74
51.01±0.28
0.990±0.001
0.72±0.04
39.56±0.98
TABLE III: Kinematic study on the Defrise flange, noise-free, native 3843 , mean ± population standard deviation. ROI-PSNR is measured at the upper lack-of-fusion stack, the designed circle-blind defect family; the lower stack behaves alike. Every row uses the same composite objective, greedy initialisation and pose parametrisation; only the reachable set changes, entering the refinement stage as bounds. Bold marks the best value per column.
Fig. 6: Defrise flange at k=80 , native 3843 . Each geometry column pairs the sagittal reconstruction (top, with LoF inset) with its source poses (bottom). The equatorial circle loses laminar contrast; the bench band already resolves the layers. Open blue circles show Tuy-greedy poses, filled red dots the optimised poses, grey lines their displacement, and dashed lines the reachable limits. Trajectory panels share the elevation scale.
ns prescan views
Selector
Prescan
8
16
32
64
cov+bundle
circle
−2.67
−1.05
−1.69
−2.10
band
−0.04
−0.81
−1.24
−0.02
cov+VCL+bundle
circle
−5.62
−2.37
−3.21
−4.86
band
−3.20
−2.35
−2.75
−1.40
TABLE IV: Prescan sensitivity on the Defrise flange at native 3843 , k=40 , Kmax=360 . Each entry is the paired PSNR difference in dB between planning on a prescan reconstruction from ns views and planning on the true volume, so zero is the oracle plan and negative is worse, mean over five seeds. Only our own selectors are listed.
VCLS
Ours
Phantom
Grid
k
Precomp.
Swap
Total
Total
Gain
Defrise flange
3843
40
98.1
0.6
98.7
5.1
19×
80
98.1
5.0
103.1
10.1
10×
Fuel nozzle
5123
40
262.1
1.2
263.3
5.6
47×
80
262.1
3.6
265.8
13.0
20×
TABLE V: Selection wall-clock in seconds at Kmax=720 on one Linux workstation with the Torch/CUDA backend, reconstruction excluded. VCLS is split into the object-specific candidate-basis precompute and the swap loop that uses it. Ours is the greedy-initialised coverage plus bundle arm. Best per row in bold, second best underlined.
Fig. 7: Measured central axial slices at k=100 , showing the quantitative 1200 -view FDK reference, the reference-coupled circular subset, same-manifold uniform sampling, and the two acquired planned arms. Gold boxes and the strip above show the complete pin region. Cyan boxes and the strip below show a high-attenuation component above the lens together with the streaks around it. All panels use the common reference grey-value window. Band acquisitions are reconstructed natively in the reference frame, and all panels use the residual-selected TV weight.
Trajectory
k
PSNR
SSIM
HFEN
circular subset
100
50.02
0.9785
16.93
uniform on band
100
45.53
0.9552
27.21
planned bundle
100
44.22
0.9495
31.31
planned full composite
100
44.09
0.9510
32.95
circular subset
400
51.49
0.9795
12.11
uniform on band
400
45.09
0.9523
27.58
TABLE VI: Measured camera experiment against the quantitative 1200 -view FDK reference. The circular subsets share the reference acquisition. Every metric is computed on the same reference-frame full-volume pair. Best values within each budget are bold and second-best values are underlined.
X-ray computed tomography (CT) reconstructs volumetric representations of objects from projection images obtained by transmitting X-rays through a target. Recent splat-based tomography, which represents a volume as a continuous distribution of 3D Gaussians, has demonstrated both high reconstruction quality and fast convergence in cone-beam sparse-view CT. However, when deployed in real CT systems with limited and non-uniform view distributions, we observe distinctive streak and strip artifacts that are far more pronounced than in conventional reconstruction methods. Through detailed analysis, we show that these artifacts primarily originate from pose inaccuracies in the acquisition geometry rather than from view sparsity itself. We revisit pose sensitivity in the splatting formulation and derive a stable gradient-based framework that jointly refines geometric parameters during reconstruction. Our study not only identifies how pose perturbations propagate through the differentiable projection operator but also reveals why splat-based CT is particularly vulnerable to geometric misalignment. The resulting formulation remains lightweight and easily integrable into existing pipelines while substantially improving reconstruction fidelity under real-world sparse-view conditions.
This work introduces a continuous soft near-orthogonality score and a resolution-aware saturated coverage objective for projection selection in region-of-interest focused cone-beam CT, grounded in Tuy's completeness theory. Replacing the binary hit-or-miss model of classical Tuy completeness with a graded, differentiable formulation preserves a direct link to achievable feature sizes while enabling both efficient approximate and exact optimisation. We establish that the underlying discrete decision problems are NP-complete via polynomial-time reductions from Set Cover, motivating a submodular greedy algorithm with proven (1−1/e) approximation guarantees and a mixed-integer linear program (MILP) that provides certified optimality bounds. The MILP serves as a quality certificate for the greedy solution rather than a competing optimiser. The primary empirical finding confirms this relationship: across a systematic benchmark spanning six target regions, multiple projection budgets, and four controlled occlusion conditions, the pooled median greedy-to-MILP objective ratio was 0.998, with a substantial fraction of cases certified globally optimal. A binary formulation is included as a diagnostic baseline; it strengthens hard directional completeness but is weaker on the continuous coverage scale. We additionally introduce Effective Spatial Resolution (ESR), a physically interpretable trajectory-level diagnostic that maps directional sampling gaps to achievable feature sizes. ESR correlates reliably with matched reconstruction quality across projection budgets and occlusion levels, providing a practical bridge between the selection stage and the image domain without requiring reconstruction.
Sparse-view computed tomography (CT) is critical for reducing radiation exposure to patients. Recent advances in radiative 3D Gaussian Splatting (3DGS) have enabled fast and accurate sparse-view CT reconstruction. Despite these algorithmic advancements, practical reconstruction fidelity remains fundamentally bounded by the quality of the captured data, raising the crucial yet underexplored problem of X-ray active view selection. Existing active view selection methods are primarily designed for natural-light scenes and fail to capture the unique geometric ambiguities and physical attenuation properties inherent in X-ray imaging. In this paper, we present Perturbed Gaussian Ensemble, an active view selection framework that integrates uncertainty modeling with sequential decision-making, tailored for X-ray Gaussian Splatting. Specifically, we identify low-density Gaussian primitives that are likely to be uncertain and apply stochastic density scaling to construct an ensemble of plausible Gaussian density fields. For each candidate projection, we measure the structural variance of the ensemble predictions and select the one with the highest variance as the next best view. Extensive experimental results on arbitrary-trajectory CT benchmarks demonstrate that our density-guided perturbation strategy effectively eliminates geometric artifacts and consistently outperforms existing baselines in progressive tomographic reconstruction under unified view selection protocols.
Yulun Wu, Ruyi Zha, Wei Cao +3
University of Illinois Urbana-Champaign · Australian National University · Johns Hopkins University