summaryrefslogtreecommitdiffstats
path: root/demos
diff options
context:
space:
mode:
authorEdoardo Pasca <edo.paskino@gmail.com>2017-08-23 16:54:01 +0100
committerEdoardo Pasca <edo.paskino@gmail.com>2017-08-23 16:54:01 +0100
commitb4ea57f337172dad7352b4ccc09f2ba0f17dbe35 (patch)
tree7a68f204b448f8e66c70f59ba6086c485a23bec0 /demos
parentbd4c6178258212bf075f19ff3638c12d46e0f819 (diff)
parentcc3a464ec587e95ddfd421cd3836a7677dfb9744 (diff)
downloadregularization-b4ea57f337172dad7352b4ccc09f2ba0f17dbe35.tar.gz
regularization-b4ea57f337172dad7352b4ccc09f2ba0f17dbe35.tar.bz2
regularization-b4ea57f337172dad7352b4ccc09f2ba0f17dbe35.tar.xz
regularization-b4ea57f337172dad7352b4ccc09f2ba0f17dbe35.zip
Merge branch 'pythonize' of https://github.com/vais-ral/CCPi-FISTA_Reconstruction into pythonize
Diffstat (limited to 'demos')
-rw-r--r--demos/exportDemoRD2Data.m35
1 files changed, 35 insertions, 0 deletions
diff --git a/demos/exportDemoRD2Data.m b/demos/exportDemoRD2Data.m
new file mode 100644
index 0000000..028353b
--- /dev/null
+++ b/demos/exportDemoRD2Data.m
@@ -0,0 +1,35 @@
+clear all
+close all
+%%
+% % adding paths
+addpath('../data/');
+addpath('../main_func/'); addpath('../main_func/regularizers_CPU/');
+addpath('../supp/');
+
+load('DendrRawData.mat') % load raw data of 3D dendritic set
+angles_rad = angles*(pi/180); % conversion to radians
+size_det = size(data_raw3D,1); % detectors dim
+angSize = size(data_raw3D, 2); % angles dim
+slices_tot = size(data_raw3D, 3); % no of slices
+recon_size = 950; % reconstruction size
+
+Sino3D = zeros(size_det, angSize, slices_tot, 'single'); % log-corrected sino
+% normalizing the data
+for jj = 1:slices_tot
+ sino = data_raw3D(:,:,jj);
+ for ii = 1:angSize
+ Sino3D(:,ii,jj) = log((flats_ar(:,jj)-darks_ar(:,jj))./(single(sino(:,ii)) - darks_ar(:,jj)));
+ end
+end
+
+Sino3D = Sino3D.*1000;
+Weights3D = single(data_raw3D); % weights for PW model
+clear data_raw3D
+
+hdf5write('DendrData.h5', '/Weights3D', Weights3D)
+hdf5write('DendrData.h5', '/Sino3D', Sino3D, 'WriteMode', 'append')
+hdf5write('DendrData.h5', '/angles_rad', angles_rad, 'WriteMode', 'append')
+hdf5write('DendrData.h5', '/size_det', size_det, 'WriteMode', 'append')
+hdf5write('DendrData.h5', '/angSize', angSize, 'WriteMode', 'append')
+hdf5write('DendrData.h5', '/slices_tot', slices_tot, 'WriteMode', 'append')
+hdf5write('DendrData.h5', '/recon_size', recon_size, 'WriteMode', 'append') \ No newline at end of file