def main():
# Initialize recon parameters
nr_projections = 320
detector_rows = 192
detector_columns = 256
recon_voxel_sizes = [64, 128, 256]
#Start pset stuff
extraTime = time()
pset,data,jobs,results = initPsets()
extraTime = time() - extraTime
for voxels in recon_voxel_sizes:
print voxels
x_voxels = y_voxels = z_voxels = voxels
# Load input data
(projections,
combined_matrix,
z_voxel_coords,
transform_matrix,
volume_weight) = fdkload(nr_projections,
detector_rows,
detector_columns,
voxels)
# Initialize 3D volume data
recon_volume = zeros((z_voxels, y_voxels, x_voxels), dtype=float32)
# Reconstruct 3D Volume from recorded 2D images
start = time()
#result = fdkcore(nr_projections, projections, combined_matrix,
# z_voxel_coords, transform_matrix, z_voxels,
# detector_rows, detector_columns, recon_volume,
# volume_weight, count_out=False)
result = fdkcoreserver(nr_projections, projections, combined_matrix,
z_voxel_coords, transform_matrix, z_voxels,
detector_rows, detector_columns, recon_volume,
volume_weight,pset,data,jobs,results)
stop = time()
print x_voxels, 'Time taken:', stop - start + extraTime, 'cpus:', result[0],'workers:',result[1]
# print 'System size', (x_voxels, y_voxels, z_voxels), 'Time taken', stop - start + extraTime
pset.halt()
def main():
# Initialize recon parameters
nr_projections = 320
detector_rows = 192
detector_columns = 256
recon_voxel_sizes = [64, 128, 256]
for voxels in recon_voxel_sizes:
x_voxels = y_voxels = z_voxels = voxels
# Load input data
(projections, combined_matrix, z_voxel_coords, transform_matrix,
volume_weight) = fdkload(nr_projections, detector_rows,
detector_columns, voxels)
# Initialize 3D volume data
recon_volume = zeros((z_voxels, y_voxels, x_voxels), dtype=float32)
# Reconstruct 3D Volume from recorded 2D images
start = time()
result = fdkcore(nr_projections,
projections,
combined_matrix,
z_voxel_coords,
transform_matrix,
z_voxels,
detector_rows,
detector_columns,
recon_volume,
volume_weight,
count_out=False)
stop = time()
print 'System size', (x_voxels, y_voxels,
z_voxels), 'Time taken', stop - start
print result
def main():
# Initialize recon parameters
nr_projections = 320
detector_rows = 192
detector_columns = 256
voxels = x_voxels = y_voxels = z_voxels =256
# Load input data
(projections,
combined_matrix,
z_voxel_coords,
transform_matrix,
volume_weight) = fdkload(nr_projections,
detector_rows,
detector_columns,
voxels)
# Initialize 3D volume data
recon_volume = zeros((z_voxels, y_voxels, x_voxels), dtype=float32)
print "Showing filtered projections"
#show(projections)
# Reconstruct 3D Volume from recorded 2D images
start = time()
result = fdkcore(nr_projections, projections, combined_matrix,
z_voxel_coords, transform_matrix, z_voxels,
detector_rows, detector_columns, recon_volume,
volume_weight, count_out=True)
stop = time()
print 'System size', (x_voxels, y_voxels, z_voxels), 'Time taken', stop - start
print "Showing reconstructed volume slices"
show(result)
def main():
# Initialize recon parameters
nr_projections = 320
detector_rows = 192
detector_columns = 256
recon_voxel_sizes = [64,128,256]
for voxels in recon_voxel_sizes:
x_voxels = y_voxels = z_voxels = voxels
# Load input data
(projections,
combined_matrix,
z_voxel_coords,
transform_matrix,
volume_weight) = fdkload(nr_projections,
detector_rows,
detector_columns,
voxels)
# Initialize 3D volume data
recon_volume = zeros((z_voxels, y_voxels, x_voxels), dtype=float32)
# Reconstruct 3D Volume from recorded 2D images
start = time()
result = fdkcore(nr_projections, projections, combined_matrix,
z_voxel_coords, transform_matrix, z_voxels,
detector_rows, detector_columns, recon_volume,
volume_weight, count_out=False)
stop = time()
print 'System size', (x_voxels, y_voxels, z_voxels), 'Time taken', stop - start
print result
