diff options
| author | Edoardo Pasca <edo.paskino@gmail.com> | 2017-08-23 16:51:18 +0100 |
|---|---|---|
| committer | Edoardo Pasca <edo.paskino@gmail.com> | 2017-10-11 15:48:34 +0100 |
| commit | a203949c84484fe2641e39451f033d20d445b1f3 (patch) | |
| tree | e3f03559387a3ba4f0ddfd3ee8013bdb0cbde28c /demos/exportDemoRD2Data.m | |
| parent | 9f8fb57e1e89c1ad200d9c7eada5c653be34db66 (diff) | |
| download | regularization-a203949c84484fe2641e39451f033d20d445b1f3.tar.gz regularization-a203949c84484fe2641e39451f033d20d445b1f3.tar.bz2 regularization-a203949c84484fe2641e39451f033d20d445b1f3.tar.xz regularization-a203949c84484fe2641e39451f033d20d445b1f3.zip | |
export/import data from hdf5
Added file to export the data from DemoRD2.m to HDF5 to pass it to Python.
Added file to import the data from DemoRD2.m from HDF5.
Diffstat (limited to 'demos/exportDemoRD2Data.m')
| -rw-r--r-- | demos/exportDemoRD2Data.m | 35 |
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 |