#----------------------------------------------------------------------- #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 . # #----------------------------------------------------------------------- import astra import numpy as np vol_geom = astra.create_vol_geom(256, 256) proj_geom = astra.create_proj_geom('parallel', 1.0, 384, np.linspace(0,np.pi,180,False)) # Create volumes # initialized to zero v0 = astra.data2d.create('-vol', vol_geom) # initialized to 3.0 v1 = astra.data2d.create('-vol', vol_geom, 3.0) # initialized to a matrix. A may be a single, double or logical (0/1) array. import scipy.io A = scipy.io.loadmat('phantom.mat')['phantom256'] v2 = astra.data2d.create('-vol', vol_geom, A) # Projection data s0 = astra.data2d.create('-sino', proj_geom) # Initialization to a scalar or a matrix also works, exactly as with a volume. # Update data # set to zero astra.data2d.store(v0, 0) # set to a matrix astra.data2d.store(v2, A) # Retrieve data R = astra.data2d.get(v2) import pylab pylab.gray() pylab.imshow(R) pylab.show() # Free memory astra.data2d.delete(v0) astra.data2d.delete(v1) astra.data2d.delete(v2) astra.data2d.delete(s0)