Zarr support

VERSION

@@ -1 +1 @@
-1.0.4
+1.0.6

recipe/meta.yaml

@@ -41,7 +41,8 @@ requirements:
     - opencv
     - opencv-contrib-python
     - pywavelets
-    - tifffile    
+    - tifffile
+    - zarr    
 
 about:
     home: https://github.com/tomography/tomocupy

setup.py

@@ -10,13 +10,14 @@
     packages=find_packages('src'),
     zip_safe=False,
     install_requires=[
-        'cupy',
+    	'cupy',
         'opencv-python',
         'h5py',
         'numexpr',
         'numpy',
         'pywavelets',
         'setuptools',  # for pkg_resources at runtime
         'tifffile',
+        'zarr',
     ]
 )

src/tomocupy/config.py

@@ -224,7 +224,7 @@ def default_parameter(func, param):
         'default': 'none',
         'type': str,
         'help': "Phase retrieval correction method",
-        'choices': ['none', 'paganin', 'Gpaganin']},
+        'choices': ['none', 'paganin', 'Gpaganin', 'FourierFilter', 'farago']},
     'energy': {
         'default': 0,
         'type': float,
@@ -240,7 +240,7 @@ def default_parameter(func, param):
     'retrieve-phase-delta-beta': {
         'default': 1500.0,
         'type': float,
-        'help': "delta/beta material for Generalized Paganin"},
+        'help': "delta/beta material for Generalized Paganin & Farago"},
     'retrieve-phase-W': {
         'default': 2e-4,
         'type': float,
@@ -249,6 +249,22 @@ def default_parameter(func, param):
         'type': utils.positive_int,
         'default': 1,
         'help': "Padding with extra slices in z for phase-retrieval filtering"},
+    'FFratio': {
+        'default': 100,
+        'type': float,
+        'help': "Shape of the Fourier Filter window"},  
+    'FFdim': {
+        'default': 2,
+        'type': float,
+        'help': "1 for sinograms, 2 for projections"},  
+    'FFlog': {
+        'default': 0,
+        'type': int,
+        'help': "0/1"},
+    'FFpad': {
+        'default': 150,
+        'type': float,
+        'help': "Padding size for the Fourier Filter"},              
 }
 
 SECTIONS['rotate-proj'] = {
@@ -395,21 +411,6 @@ def default_parameter(func, param):
         'type': int,
         'default': -1,
         'help': "End row to find the rotation center"},
-    'dtype': {
-        'default': 'float32',
-        'type': str,
-        'choices': ['float32', 'float16'],
-        'help': "Data type used for reconstruction. Note float16 works with power of 2 sizes.", },
-    'save-format': {
-        'default': 'tiff',
-        'type': str,
-        'help': "Output format",
-        'choices': ['tiff', 'h5', 'h5sino', 'h5nolinks']},
-    'clear-folder': {
-        'default': 'False',
-        'type': str,
-        'help': "Clear output folder before reconstruction",
-        'choices': ['True', 'False']},
     'fbp-filter': {
         'default': 'parzen',
         'type': str,
@@ -446,6 +447,41 @@ def default_parameter(func, param):
         'type': float,
         'help': "Pixel size [microns]"},
 }
+SECTIONS['output'] ={
+    'dtype': {
+        'default': 'float32',
+        'type': str,
+        'choices': ['float32', 'float16'],
+        'help': "Data type used for reconstruction. Note float16 works with power of 2 sizes.", },
+    'save-format': {
+        'default': 'tiff',
+        'type': str,
+        'help': "Output format",
+        'choices': ['tiff', 'h5', 'h5sino', 'h5nolinks', 'zarr']},
+    'zarr-compression': {
+        'default': 'blosclz',
+        'type': str,
+        'help': "ZARR compression format",
+        'choices': ['blosclz', 'lz4', 'zstd']},    
+    'zarr-chunk': {
+        'default': '8,64,64',
+        'type': str,
+        'help': "ZARR chunk size"},
+    'large-data': {
+        'default': False,
+        'type': bool,
+        'help': "If Active it computes ldchunk chunks of angular projections"},
+    'ldchunk': {
+        'default': 10,
+        'type': int,
+        'help': "Number of angular chunks for large-data"},                  
+    'clear-folder': {
+        'default': 'False',
+        'type': str,
+        'help': "Clear output folder before reconstruction",
+        'choices': ['True', 'False']}
+}
+
 
 SECTIONS['beam-hardening'] = {
     'beam-hardening-method': {
@@ -564,9 +600,9 @@ def default_parameter(func, param):
 
 
 RECON_PARAMS = ('file-reading', 'remove-stripe',
-                'reconstruction', 'fw', 'ti', 'vo-all', 'lamino', 'reconstruction-types', 'beam-hardening')
+                'reconstruction', 'fw', 'ti', 'vo-all', 'lamino', 'reconstruction-types', 'beam-hardening', 'output')
 RECON_STEPS_PARAMS = ('file-reading', 'remove-stripe', 'reconstruction',
-                      'retrieve-phase', 'fw', 'ti', 'vo-all', 'lamino', 'reconstruction-steps-types', 'rotate-proj', 'beam-hardening')
+                      'retrieve-phase', 'fw', 'ti', 'vo-all', 'lamino', 'reconstruction-steps-types', 'rotate-proj', 'beam-hardening', 'output')
 
 NICE_NAMES = ('General', 'File reading', 'Remove stripe',
               'Remove stripe FW', 'Remove stripe Titarenko', 'Remove stripe Vo' 'Retrieve phase', 'Reconstruction')
@@ -773,7 +809,7 @@ def update_hdf_process(fname, args=None, sections=None):
                         hdf_file.require_dataset(dataset, shape=(1,), dtype=dt)
                         log.info(name + ': ' + str(value))
                         try:
-                            hdf_file[dataset][0] = np.string_(str(value))
+                            hdf_file[dataset][0] = np.bytes_(str(value))
                         except TypeError:
                             log.error(
                                 "Could not convert value {}".format(value))

src/tomocupy/dataio/reader.py

@@ -84,7 +84,7 @@ def init_sizes(self):
 
         # read data sizes and projection angles with a reader
         sizes = self.read_sizes()
-        theta = self.read_theta()
+        theta = self.read_theta(sizes['nproji'])
         nproji = sizes['nproji']
         nzi = sizes['nzi']
         ni = sizes['ni']
@@ -306,12 +306,15 @@ def read_sizes(self):
 
         return sizes
 
-    def read_theta(self):
+    def read_theta(self, projections):
         """Read projection angles (in radians)"""
 
         with h5py.File(args.file_name) as file_in:
-            theta = file_in['/exchange/theta'][:].astype('float32')/180*np.pi
-
+                if '/exchange/theta' in file_in:
+                        theta = file_in['/exchange/theta'][:].astype('float32') / 180 * np.pi
+                else:
+                        # If 'theta' doesn't exist, calculate it over the range [0, proj]
+                        theta = np.linspace(0, np.pi, projections, dtype='float32')
         return theta
 
     def read_data_chunk_to_queue(self, data_queue, ids_proj, st_z, end_z, st_n, end_n, id_z, in_dtype):