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