Various IPAC related modelling and POF Videos etc

Author: eddrial

animation/change_gap_model.py

@@ -208,7 +208,7 @@ def __init__(self, drawobj, my_wrdobj, count, mod_params):
 
         self.show()
         scshframe = _screenshot()
-        image.write_png('d:\Profile\oqb\Desktop\presentations\POF2025\Animation\gap_apple\ivue32_{}_g{}_s{}.png'.format(mod_params.shiftmode,mod_params.gap,mod_params.rowshift), scshframe)
+        image.write_png('d:\Profile\oqb\Desktop\presentations\POF2025\Animation\gap_plain_apple\ivue32_{}_g{}_s{}.png'.format(mod_params.shiftmode,mod_params.gap,mod_params.rowshift), scshframe)
         self.close()
 
 
@@ -268,8 +268,8 @@ def on_draw(self, event):
     #create ID
 
     #shift loop
-    for mode in ['circular', 'linear']:
-        for shift in range(32,33,16):
+    for mode in ['circular']:
+        for shift in range(0,17,16):
             for gap in range (6, 326,1):
                 ggap = (gap-6)/8+6
                 rd.UtiDelAll()
@@ -286,7 +286,7 @@ def on_draw(self, event):
 
                 #a = ArbAPPLE(model_parameters = a_param)
 
-                a = id1.compensatedAPPLEv2(model_parameters = a_param)
+                a = id1.plainAPPLE(model_parameters = a_param)
 
                 #rd.ObjDrwOpenGL(a.cont.radobj)
 

animation/create_vid.py

@@ -7,19 +7,20 @@
 import os
 import numpy as np
 
-image_folder = 'd:\Profile\oqb\Desktop\presentations\POF2025\Animation\shift_plain_apple'
-video_name = 'circ_shift_plain.avi'
+image_folder = 'd:\Profile\oqb\Desktop\presentations\POF2025\Animation\gap_plain_apple'
+video_name = 'v_gap_plain.avi'
 
 
-image_ref = '{}\ivue32_circular_0.png'.format(image_folder)
+image_ref = '{}\ivue32_circular_g6.0_s16.png'.format(image_folder)
 frame = cv2.imread(image_ref)
 height, width, layers = frame.shape
 
 video = cv2.VideoWriter(os.path.join(image_folder,video_name), 0, 40    , (width,height))
 
-for i in np.arange(362):
+for i in np.arange(326):
     print(i)
-    image = 'ivue32_circular_{}.png'.format(i)
+    gap = (i-6)/8+6
+    image = 'ivue32_circular_g{}_s16.png'.format(gap)
     video.write(cv2.imread(os.path.join(image_folder, image)))
 
 cv2.destroyAllWindows()

animation/create_vid_force_vectors.py

@@ -7,61 +7,33 @@
 import os
 import numpy as np
 
-image_folder = 'd:\Profile\oqb\Desktop\presentations\POF2025\Schaefer\pics'
+image_folder = 'd:\Profile\oqb\Desktop\presentations\POF2025\Forces\Plots'
 
-vids = 'interpolate'
 
-if vids == 'non_interpolate':
+
+for result in ['forcevectorsapple_h',
+               'forcevectorsapple_v',
+               'forcevectorscompapple_h',
+               'forcevectorscompapple_v']:
 
-    for result in ['bfield',
-                   'flux_density',
-                   'trajectory']:
-        
-        
-        
-        image_ref = '{}/{}_gap_6.00_shift_0.00.png'.format(image_folder, result)
-        frame = cv2.imread(image_ref)
-        height, width, layers = frame.shape
-        
-        
-        for s in ['0.00','-16.00','16.00','32.00','-32.00']:
-            video_name = '{}_s{}_gap.avi'.format(result,s)
-            video = cv2.VideoWriter(os.path.join(image_folder,'vids',video_name), 0, 0.625, (width,height))
-        
-            for i in ['6.00','15.26','24.53','33.79','43.05']:
-                image = '{}/{}_gap_{}_shift_{}.png'.format(image_folder,result,i, s)
-                frame = cv2.imread(image)
-                resized_frame = cv2.resize(frame,(width,height))
-                #cv2.imwrite('')
-                print('{} gives frame size {}, originally {}'.format(i,resized_frame.shape, frame.shape))
-                video.write(resized_frame)
 
-    cv2.destroyAllWindows()
-    video.release()
 
-else:
-    for result in ['flux_density_distribution',
-                   'power_distro']:
-        
-        
-        
-        image_ref = '{}/{}_gap_6.00_shift_0.00_interpolated.png'.format(image_folder, result)
-        frame = cv2.imread(image_ref)
-        height, width, layers = frame.shape
-        
-        
-        for s in ['0.00','-16.00','16.00','32.00','-32.00']:
-            video_name = '{}_s{}_gap_interpolated.avi'.format(result,s)
-            video = cv2.VideoWriter(os.path.join(image_folder,'vids',video_name), 0, 0.625, (width,height))
-        
-            for i in ['6.00','15.26','24.53','33.79','43.05']:
-                image = '{}/{}_gap_{}_shift_{}_interpolated.png'.format(image_folder,result,i, s)
-                frame = cv2.imread(image)
-                resized_frame = cv2.resize(frame,(width,height))
-                #cv2.imwrite('')
-                print('{} gives frame size {}, originally {}'.format(i,resized_frame.shape, frame.shape))
-                video.write(resized_frame)
+    image_ref = '{}\\forcevectorsapple_v6.png'.format(image_folder, result)
+    frame = cv2.imread(image_ref)
+    height, width, layers = frame.shape
 
-    cv2.destroyAllWindows()
-    video.release()
-    
+
+    video_name = '{}gap.avi'.format(result)
+    video = cv2.VideoWriter(os.path.join(image_folder,video_name), 0, 5, (width,height))
+
+    for i in range(41):
+        image = '{}\\{}{}.png'.format(image_folder,result,i+6)
+        frame = cv2.imread(image)
+#        resized_frame = cv2.resize(frame,(width,height))
+        #cv2.imwrite('')
+ #       print('{} gives frame size {}, originally {}'.format(i,resized_frame.shape, frame.shape))
+        video.write(frame)
+
+cv2.destroyAllWindows()
+video.release()
+

animation/create_vid_stefan.py

@@ -0,0 +1,43 @@
+'''
+Created on Apr 7, 2025
+
+@author: oqb
+'''
+import cv2
+import os
+import numpy as np
+
+image_folder = 'd:\\Profile\\oqb\\Desktop\\presentations\\POF2025\\Schaefer\\final_plots'
+
+fin_plots = {
+    'bfields' : 'bfield',
+    'flux_densities' : 'fd',
+    'flux_density_distributions' : 'fdd',
+    'power_distributions' : 'pd',
+    'trajectories': 'traj'
+    }
+
+for result in fin_plots:
+    
+    
+    
+    image_ref = '{}\\{}\\{}_gap_7.00_shift_0.00.png'.format(image_folder, result, fin_plots[result])
+    frame = cv2.imread(image_ref)
+    height, width, layers = frame.shape
+    
+    
+    for s in ['0.00','-8.00','8.00','-16.00','16.00']:
+        video_name = '{}_s{}_gap.avi'.format(result,s)
+        video = cv2.VideoWriter(os.path.join(image_folder,'final_vids',video_name), 0, 5, (width,height))
+    
+        for i in range(7,48):
+            image = '{}\\{}\\{}_gap_{}.00_shift_{}.png'.format(image_folder,result,fin_plots[result],i, s)
+            frame = cv2.imread(image)
+            resized_frame = cv2.resize(frame,(width,height))
+            #cv2.imwrite('')
+            print('{} gives frame size {}, originally {}'.format(i,resized_frame.shape, frame.shape))
+            video.write(resized_frame)
+
+cv2.destroyAllWindows()
+
+video.release()