def test_project():
assert_allclose(project(read_file("obj.npy"), read_file("angle.npy")), read_file("proj.npy"), rtol=1e-2)
except OSError as exc: # Python >2.5
if exc.errno == errno.EEXIST and os.path.isdir(path):
pass
else: raise
'''
print('Creating Object...')
obj = shepp3d((128, 128, 128))
numpy.save('./data/largeTestObj.npy', obj)
print('Creating Angle...')
theta = angles(200)
numpy.save('./data/largeTestTheta.npy', theta)
print('Creating Projection...')
data = project(obj, theta)
numpy.save('./data/largeTestData.npy', data)
print('Creating ospml_hybrid...')
pmlrecon = recon(data, theta, algorithm='ospml_hybrid', num_iter=2)
numpy.save('./data/pml/largeTestRecon_OH.npy', pmlrecon)
print('Creating ospml_quad...')
pmlrecon = recon(data, theta, algorithm='ospml_quad', num_iter=2)
numpy.save('./data/pml/largeTestRecon_OQ.npy', pmlrecon)
print('Creating pml_hybrid...')
pmlrecon = recon(data, theta, algorithm='pml_hybrid', num_iter=2)
numpy.save('./data/pml/largeTestRecon_PH.npy', pmlrecon)
print('Creating pml_quad...')
def test_project(self):
assert_allclose(project(read_file('obj.npy'), read_file('angle.npy')),
read_file('proj.npy'),
rtol=1e-2)
def test_project():
assert_allclose(
project(read_file('obj.npy'), read_file('angle.npy')),
read_file('proj.npy'), rtol=1e-2)
if exc.errno == errno.EEXIST and os.path.isdir(path):
pass
else: raise
'''
print('Creating Object...')
obj=shepp3d((128,128,128))
numpy.save('./data/largeTestObj.npy',obj)
print('Creating Angle...')
theta=angles(200)
numpy.save('./data/largeTestTheta.npy',theta)
print('Creating Projection...')
data=project(obj,theta)
numpy.save('./data/largeTestData.npy',data)
print('Creating ospml_hybrid...')
pmlrecon=recon(data, theta, algorithm='ospml_hybrid', num_iter=2)
numpy.save('./data/pml/largeTestRecon_OH.npy',pmlrecon)
print('Creating ospml_quad...')
pmlrecon=recon(data, theta, algorithm='ospml_quad', num_iter=2)
numpy.save('./data/pml/largeTestRecon_OQ.npy',pmlrecon)
print('Creating pml_hybrid...')
pmlrecon=recon(data, theta, algorithm='pml_hybrid', num_iter=2)
numpy.save('./data/pml/largeTestRecon_PH.npy',pmlrecon)
print('Creating pml_quad...')
