Hznnn

How to add Poisson noise to projection data？

12

``` dataset = dicom.read_file('./L067_FD_1_1.CT.0001.0001.2015.12.22.18.09.40.840353.358074219.IMA') img = dataset.pixel_array.astype(np.float32) RescaleSlope = dataset.RescaleSlope RescaleIntercept = dataset.RescaleIntercept CT_img = img * RescaleSlope + RescaleIntercept reco_space = odl.uniform_discr(min_pt=[-128, -128], max_pt=[128, 128], shape=[512, 512], dtype='float32') angle_partition...

Read data in mayo_rec

1

``` if validation: fi = DATA_FOLDER + 'L286_FD_3_1.CT.0002.0201.2015.12.22.18.22.49.651226.358225786.npy' else: fi = file_loader.next_file() data = np.load(fi) phantom = space.element(np.rot90(data, -1)) phantom /= 1000.0 # convert go g/cm^3 data = operator(phantom) data...