cpos for I_c

hrt_pipe.py

@@ -93,7 +93,7 @@ def phihrt_pipe(data_f,dark_f,flat_f,norm_f = True, clean_f = False, sigma = 59,
     12. normalise to quiet sun
     13. calibration
         a) ghost correction - not implemented yet
-        b) cross talk correction - not implemented yet
+        b) cross talk correction
     14. rte inversion with cmilos
         a) output rte data products to fits file
 
@@ -291,7 +291,7 @@ def phihrt_pipe(data_f,dark_f,flat_f,norm_f = True, clean_f = False, sigma = 59,
     # TODO: Could check data dimensions? As an extra fail safe before progressing?
     #-----------------
     
-    print(voltagesData_arr)
+    #print(voltagesData_arr)
 
     #-----------------
     # READ FLAT FIELDS
@@ -637,7 +637,7 @@ def phihrt_pipe(data_f,dark_f,flat_f,norm_f = True, clean_f = False, sigma = 59,
 
         for scan in range(data_shape[-1]):
             
-            I_c = np.mean(data[512:1536,512:1536,0,0,int(scan)]) #mean of central 1k x 1k of continuum stokes I - all data cpos shifted to 0
+            I_c = np.mean(data[512:1536,512:1536,0,cpos_arr[0],int(scan)]) #mean of central 1k x 1k of continuum stokes I
             data[:,:,:,:,scan] = data[:,:,:,:,scan]/I_c
        
         printc('--------------------------------------------------------------',bcolors.OKGREEN)

run.py

 from hrt_pipe import phihrt_pipe
 import numpy as np
 
-sciencedata_fits_filenames = ['solo_L0_phi-hrt-ilam_0667414905_V202103221851C_0142230602.fits', 'solo_L0_phi-hrt-ilam_0667415054_V202103221851C_0142230603.fits', 'solo_L0_phi-hrt-ilam_0667415205_V202103221851C_0142230604.fits', 'solo_L0_phi-hrt-ilam_0667415354_V202103221851C_0142230605.fits', 'solo_L0_phi-hrt-ilam_0667415505_V202103221851C_0142230606.fits', 'solo_L0_phi-hrt-ilam_0667415654_V202103221851C_0142230607.fits', 'solo_L0_phi-hrt-ilam_0667415805_V202103221851C_0142230608.fits']#['../fits_files/solo_L0_phi-hrt-ilam_0667414748_V202103221851C_0142230201.fits']
+sciencedata_fits_filenames = ['solo_L0_phi-hrt-ilam_0667414748_V202103221851C_0142230201.fits']#['solo_L0_phi-hrt-ilam_0667414905_V202103221851C_0142230602.fits', 'solo_L0_phi-hrt-ilam_0667415054_V202103221851C_0142230603.fits', 'solo_L0_phi-hrt-ilam_0667415205_V202103221851C_0142230604.fits', 'solo_L0_phi-hrt-ilam_0667415354_V202103221851C_0142230605.fits', 'solo_L0_phi-hrt-ilam_0667415505_V202103221851C_0142230606.fits', 'solo_L0_phi-hrt-ilam_0667415654_V202103221851C_0142230607.fits', 'solo_L0_phi-hrt-ilam_0667415805_V202103221851C_0142230608.fits']#['../fits_files/solo_L0_phi-hrt-ilam_0667414748_V202103221851C_0142230201.fits']
 flatfield_fits_filename = '../fits_files/solo_L0_phi-hrt-flat_0667134081_V202103221851C_0162201100.fits'
 darkfield_fits_filename = '../fits_files/solo_L0_phi-fdt-ilam_20200228T155100_V202002281636_0022210004_000.fits'
 
@@ -27,5 +27,5 @@ c_talk_params[1,0] = q_int
 c_talk_params[1,1] = u_int
 c_talk_params[1,2] = v_int
 
-data = phihrt_pipe(sciencedata_fits_filenames[1], darkfield_fits_filename, flatfield_fits_filename, norm_stokes = True, prefilter_f = prefilter_f, 
-                    clean_f = False, ctalk_params = c_talk_params, ItoQUV = False, out_demod_file = False, out_dir = '/data/slam/home/sinjan/hrt_pipe_results/stp-136_ctalk/', rte = 'False')
\ No newline at end of file
+phihrt_pipe(sciencedata_fits_filenames, darkfield_fits_filename, flatfield_fits_filename, norm_stokes = True, prefilter_f = None, 
+                    clean_f = True, ctalk_params = c_talk_params, ItoQUV = True, out_demod_file = True, out_dir = '/data/slam/home/sinjan/hrt_pipe_results/', rte = 'False')
\ No newline at end of file