#-----------------------------------------------------------------------
#Copyright 2013 Centrum Wiskunde & Informatica, Amsterdam
#
#Author: Daniel M. Pelt
#Contact: D.M.Pelt@cwi.nl
#Website: http://dmpelt.github.io/pyastratoolbox/
#
#
#This file is part of the Python interface to the
#All Scale Tomographic Reconstruction Antwerp Toolbox ("ASTRA Toolbox").
#
#The Python interface to the ASTRA Toolbox is free software: you can redistribute it and/or modify
#it under the terms of the GNU General Public License as published by
#the Free Software Foundation, either version 3 of the License, or
#(at your option) any later version.
#
#The Python interface to the ASTRA Toolbox is distributed in the hope that it will be useful,
#but WITHOUT ANY WARRANTY; without even the implied warranty of
#MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
#GNU General Public License for more details.
#
#You should have received a copy of the GNU General Public License
#along with the Python interface to the ASTRA Toolbox. If not, see <http://www.gnu.org/licenses/>.
#
#-----------------------------------------------------------------------
from . import data3d_c as d
import numpy as np
[docs]def create(datatype,geometry,data=None):
"""Create a 3D object.
:param datatype: Data object type, '-vol' or '-sino'.
:type datatype: :class:`string`
:param geometry: Volume or projection geometry.
:type geometry: :class:`dict`
:param data: Data to fill the constructed object with, either a scalar or array.
:type data: :class:`float` or :class:`numpy.ndarray`
:returns: :class:`int` -- the ID of the constructed object.
"""
return d.create(datatype,geometry,data)
[docs]def link(datatype, geometry, data):
"""Link a 3D numpy array with the toolbox.
:param datatype: Data object type, '-vol' or '-sino'.
:type datatype: :class:`string`
:param geometry: Volume or projection geometry.
:type geometry: :class:`dict`
:param data: Numpy array to link
:type data: :class:`numpy.ndarray`
:returns: :class:`int` -- the ID of the constructed object.
"""
if not isinstance(data,np.ndarray):
raise ValueError("Input should be a numpy array")
if not data.dtype==np.float32:
raise ValueError("Numpy array should be float32")
if not (data.flags['C_CONTIGUOUS'] and data.flags['ALIGNED']):
raise ValueError("Numpy array should be C_CONTIGUOUS and ALIGNED")
return d.create(datatype,geometry,data,True)
[docs]def get(i):
"""Get a 3D object.
:param i: ID of object to get.
:type i: :class:`int`
:returns: :class:`numpy.ndarray` -- The object data.
"""
return d.get(i)
[docs]def get_shared(i):
"""Get a 3D object with memory shared between the ASTRA toolbox and numpy array.
:param i: ID of object to get.
:type i: :class:`int`
:returns: :class:`numpy.ndarray` -- The object data.
"""
return d.get_shared(i)
[docs]def get_single(i):
"""Get a 3D object in single precision.
:param i: ID of object to get.
:type i: :class:`int`
:returns: :class:`numpy.ndarray` -- The object data.
"""
return g.get_single(i)
[docs]def store(i,data):
"""Fill existing 3D object with data.
:param i: ID of object to fill.
:type i: :class:`int`
:param data: Data to fill the object with, either a scalar or array.
:type data: :class:`float` or :class:`numpy.ndarray`
"""
return d.store(i,data)
[docs]def get_geometry(i):
"""Get the geometry of a 3D object.
:param i: ID of object.
:type i: :class:`int`
:returns: :class:`dict` -- The geometry of object with ID ``i``.
"""
return d.get_geometry(i)
[docs]def change_geometry(i, geometry):
"""Change the geometry of a 3D object.
:param i: ID of object.
:type i: :class:`int`
:param geometry: Volume or projection geometry.
:type geometry: :class:`dict`
"""
return d.change_geometry(i, geometry)
[docs]def dimensions(i):
"""Get dimensions of a 3D object.
:param i: ID of object.
:type i: :class:`int`
:returns: :class:`tuple` -- dimensions of object with ID ``i``.
"""
return d.dimensions(i)
[docs]def delete(ids):
"""Delete a 2D object.
:param ids: ID or list of ID's to delete.
:type ids: :class:`int` or :class:`list`
"""
return d.delete(ids)
[docs]def clear():
"""Clear all 3D data objects."""
return d.clear()
[docs]def info():
"""Print info on 3D objects in memory."""
return d.info()