From 9bb7b68b15fa3d8ee68223398ea6d0859ee9c10b Mon Sep 17 00:00:00 2001 From: Michael-T-McCann Date: Wed, 11 Sep 2024 09:11:40 -0600 Subject: [PATCH] Remove unfinished example --- examples/scripts/ct_large_projection.py | 42 ------------------------- 1 file changed, 42 deletions(-) delete mode 100644 examples/scripts/ct_large_projection.py diff --git a/examples/scripts/ct_large_projection.py b/examples/scripts/ct_large_projection.py deleted file mode 100644 index 5067698b..00000000 --- a/examples/scripts/ct_large_projection.py +++ /dev/null @@ -1,42 +0,0 @@ -#!/usr/bin/env python -# -*- coding: utf-8 -*- -# This file is part of the SCICO package. Details of the copyright -# and user license can be found in the 'LICENSE.txt' file distributed -# with the package. - -r""" -Large-scale CT Projection -========================= - -This example demonstrates using SCICO's X-ray projector on a large-scale -volume. - -""" - -import numpy as np - -import jax - -from scico.examples import create_block_phantom -from scico.linop.xray import XRayTransform3D - -N = 1000 -num_views = 10 - -in_shape = (N, N, N) -x = create_block_phantom(in_shape) - -det_shape = (N, N) - -rot_X = 90.0 - 16.0 -rot_Y = np.linspace(0, 180, num_views, endpoint=False) -angles = np.stack(np.broadcast_arrays(rot_X, rot_Y), axis=-1) -matrices = XRayTransform3D.matrices_from_euler_angles( - in_shape, det_shape, "XY", angles, degrees=True -) - - -H = XRayTransform3D(in_shape, matrices, det_shape) - -proj = H @ x -jax.block_until_ready(proj)