examples/python/Vision_files_preprocess.py: add View offset to remove any rotational shifts between Siemens and STIR Vision reconstructions

Author: NicoleJurjew

examples/python/Vision_files_preprocess.py

@@ -149,6 +149,10 @@ def process_files(prompts_header_filename,
     ## the crystal depth of interaction (DOI) from 7mm to 10mm to minimize the differences.
     if apply_DOI_adaption: DOI_adaption(prompts_from_e7, 10)
 
+    ###################### View Offset ADAPTION ############################
+
+    view_offset_adaption(prompts_from_e7, 0.0324)
+
     ## write to file so you can use it later
     prompts_from_e7.write_to_file(os.path.join(STIR_output_folder, prompts_filename_STIR_corr_DOI))
 
@@ -484,6 +488,15 @@ def DOI_adaption(projdata, DOI_new):
     DOI = proj_info.get_scanner().get_average_depth_of_interaction()
     print('New Depth of interaction:', DOI)
 
+def view_offset_adaption(projdata, view_offset):
+    proj_info = projdata.get_proj_data_info()
+
+    VO = proj_info.get_scanner().get_intrinsic_azimuthal_tilt()
+    print('Current view offset (rad):', VO)
+    proj_info.get_scanner().set_intrinsic_azimuthal_tilt(view_offset)
+    VO = proj_info.get_scanner().get_average_depth_of_interaction()
+    print('New Depth of interaction(rad):', VO)
+
 def check_if_compressed(header_filename):
     with open(header_filename) as f:
         data = f.read()