1 files changed, 60 insertions, 0 deletions
diff --git a/src/operators/UfoStackedProjector.py b/src/operators/UfoStackedProjector.py
new file mode 100644
index 0000000..e2e65fd
--- /dev/null
+++ b/src/operators/UfoStackedProjector.py
@@ -0,0 +1,60 @@
+# ========================================================================
+# Copyright 2019 Science Technology Facilities Council
+# Copyright 2019 University of Manchester
+#
+# This work is part of the Core Imaging Library developed by Science Technology
+# Facilities Council and University of Manchester
+#
+# Licensed under the Apache License, Version 2.0 (the "License");
+# you may not use this file except in compliance with the License.
+# You may obtain a copy of the License at
+#
+#         http://www.apache.org/licenses/LICENSE-2.0.txt
+#
+# Unless required by applicable law or agreed to in writing, software
+# distributed under the License is distributed on an "AS IS" BASIS,
+# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+# See the License for the specific language governing permissions and
+# limitations under the License.
+#
+# =========================================================================
+
+from cil.optimisation.operators import LinearOperator
+from src.processors import UfoStackedForwardProjector, UfoStackedBackProjector
+
+class UfoStackedProjector(LinearOperator):
+    """UFO projector modified to use DataSet and geometry."""
+
+    def __init__(self, geomv, geomp, precision_mode='single', stack_num=2):
+
+        super(UfoStackedProjector, self).__init__(domain_geometry=geomv, range_geometry=geomp)
+
+        self.sinogram_geometry = geomp
+        self.volume_geometry = geomv
+
+        self.fp = UfoStackedForwardProjector(volume_geometry=geomv, sinogram_geometry=geomp,
+                                             precision_mode=precision_mode, stack_num=stack_num)
+        self.bp = UfoStackedBackProjector(volume_geometry=geomv, sinogram_geometry=geomp,
+                                          precision_mode=precision_mode, stack_num=stack_num)
+
+    def direct(self, IM, out=None):
+        self.fp.set_input(IM)
+
+        if out is None:
+            return self.fp.get_output()
+        else:
+            out.fill(self.fp.get_output())
+
+    def adjoint(self, DATA, out=None):
+        self.bp.set_input(DATA)
+
+        if out is None:
+            return self.bp.get_output()
+        else:
+            out.fill(self.bp.get_output())
+
+    def domain_geometry(self):
+        return self.volume_geometry
+
+    def range_geometry(self):
+        return self.sinogram_geometry