Mercurial > repos > rv43 > tomo

diff tomo.py @ 24:817a8101d2e4 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
"planemo upload for repository https://github.com/rolfverberg/galaxytools commit 3d7910473008ec94bc0b4580d62e61310e69917d"
author rv43
date Mon, 18 Apr 2022 20:36:53 +0000
parents 3e336f5b6402
children 7791b6e7d4e6
line wrap: on
line diff
--- a/tomo.py	Mon Apr 18 20:02:08 2022 +0000
+++ b/tomo.py	Mon Apr 18 20:36:53 2022 +0000
@@ -1183,7 +1183,7 @@
         center = sinogram.shape[1]/2
 
         # try automatic center finding routines for initial value
-        print('OK2')
+        print(f'OK2 numcore = {num_core}')
         tomo_center = tomopy.find_center_vo(sinogram, ncore=num_core)
         center_offset_vo = tomo_center-center
         print(f'center_offset_vo = {center_offset_vo}')
@@ -1198,7 +1198,7 @@
             title = f'edges row{row} center offset{center_offset_vo:.2f} Vo'
             self._plotEdgesOnePlane(recon_plane, title, path='find_center_pngs')
             del recon_plane
-            print(f'center_type_selector = {center_type_selector}')
+            print(f'center_type_selector = {galaxy_param["center_type_selector"]}')
             if not galaxy_param['center_type_selector']:
                 del sinogram_T
                 return float(center_offset_vo)
@@ -1800,7 +1800,7 @@
                 if self.save_plots_only:
                     msnc.clearFig(f'theta={theta_start}')
             # center_stack order: row,theta,column
-            center_offset = self._findOenterOnePlane(center_stack[row,:,:], row, thetas_deg,
+            center_offset = self._findCenterOnePlane(center_stack[row,:,:], row, thetas_deg,
                     eff_pixel_size, cross_sectional_dim, num_core=num_core,
                     galaxy_param=galaxy_param)
         logging.info(f'upper_center_offset = {center_offset:.2f}')