Merge pull request #47 from FZJ-INM1-BDA/bugfix_typecastMap64fTo32f

bugfix: type cast float64 maps to float32

Author: xgui3783

app/service/request_utils.py

@@ -507,6 +507,8 @@ def get_path_to_regional_map(query_id, roi, space_of_interest):
     # cache fails, fetch from source
     def save_new_nii(cached_fullpath):
         import nibabel as nib
+        import numpy as np
+
         # fix regional_map if necessary
         regional_map.image.header.set_xyzt_units('mm', 'sec')
 
@@ -515,7 +517,22 @@ def save_new_nii(cached_fullpath):
 
         # num channel
         regional_map.image.header['dim'][5] = 1
-        nib.save(regional_map.image, cached_fullpath)
+
+        # cast type float64 to float32
+        if regional_map.image.header.get_data_dtype() == np.float64:
+            fdata=regional_map.image.get_fdata()
+            new_data=fdata.astype(np.float32)
+            regional_map.image.set_data_dtype(np.float32)
+
+            if regional_map.image.header['sizeof_hdr'] == 348:
+                new_image=nib.Nifti1Image(new_data, regional_map.image.affine, regional_map.image.header)
+            elif regional_map.image.header['sizeof_hdr'] == 540:
+                new_image=nib.Nifti2Image(new_data, regional_map.image.affine, regional_map.image.header)
+            else:
+                raise IOError('regional map has incorrect sizeof_hdr')
+            nib.save(new_image, cached_fullpath)
+        else:
+            nib.save(regional_map.image, cached_fullpath)
 
     return get_cached_file(cached_filename, save_new_nii)
 

requirements/dev.txt

@@ -1,3 +1,2 @@
 -r common.txt
-git+https://github.com/FZJ-INM1-BDA/siibra-python@develop
-# git+https://github.com/FZJ-INM1-BDA/siibra-python@refactor/dataset
+git+https://github.com/FZJ-INM1-BDA/siibra-python@main

requirements/prod.txt

@@ -1,2 +1,2 @@
 -r common.txt
-siibra==0.1a8
+siibra==0.3a0