mesoscope_motion_correction_final_ploy.py

import matplotlib.figure as mplt_figure
from numpy.core.fromnumeric import size
import pandas as pd
import imageio

#######################################################
#      sst = 5
fig = mplt_figure.Figure(figsize=(60,10))
axes = [fig.add_subplot(1, 6, ii) for ii in range(1,7,1)]
# fig.suptitle('Multiple experiments, smooth_sigma_time = 5', size=35)

#######################################################
#      Row 1
row_one = pd.read_csv("./final/1076633167_1076808563_suite2p_rigid_motion_transform.csv")
row_one_old = pd.read_csv("./final/1076808563_suite2p_rigid_motion_transform.csv")

axes[0].imshow(imageio.imread('./final/1076808563_suite2p_average_projection.png'), cmap='gray')
axes[1].imshow(imageio.imread('./final/1076808563_suite2p_maximum_projection.png'), cmap='gray')

axes[2].imshow(imageio.imread('./final/1076633167_1076808563_suite2p_average_projection.png'), cmap='gray')
axes[3].imshow(imageio.imread('./final/1076633167_1076808563_suite2p_maximum_projection.png'), cmap='gray')

axes[4].plot(row_one_old['x'], alpha=0.5, color='blue', label='smooth_sigma_time = 0.0')
axes[4].plot(row_one['x'], alpha=0.5, color='red', label='smooth_sigma_time = 3.0')
axes[4].set_ylabel('x', size=30)
axes[4].set_xlabel('Frame Number', size=20)
axes[4].tick_params(axis='y', labelsize=30)

axes[5].plot(row_one_old['y'], alpha=0.5, color='blue', label='smooth_sigma_time = 0.0')
axes[5].plot(row_one['y'], alpha=0.5, color='red', label='smooth_sigma_time = 3.0')
axes[5].set_ylabel('y', size=30)
axes[5].set_xlabel('Frame Number', size=20)
axes[5].tick_params(axis='y', labelsize=30)

axes[5].legend(loc='upper right', prop={'size': 20})

axes[0].set_title("Default Avg. Proj.", size=40)
axes[1].set_title("Default Max. Proj.", size=40)
axes[2].set_title("Avg. Proj.", size=40)
axes[3].set_title("Max. Proj.", size=40)
axes[4].set_title("Session 1076633167 Experiment 1076808563", size=40)

fig.tight_layout()
fig.savefig('./final/waynes_experiment.png')


# #######################################################
# #      sst = 5
# fig = mplt_figure.Figure(figsize=(60,45))
# axes = [fig.add_subplot(6, 6, ii) for ii in range(1,37,1)]
# # fig.suptitle('Multiple experiments, smooth_sigma_time = 5', size=35)
# #######################################################
# #      Row 1
# row_one = pd.read_csv("./final/1110070263_1110190110_suite2p_rigid_motion_transform.csv")
# row_one_old = pd.read_csv("./final/1110190110_suite2p_rigid_motion_transform.csv")

# axes[0].imshow(imageio.imread('./final/1110190110_suite2p_average_projection.png'), cmap='gray')
# axes[1].imshow(imageio.imread('./final/1110190110_suite2p_maximum_projection.png'), cmap='gray')

# axes[2].imshow(imageio.imread('./final/1110070263_1110190110_suite2p_average_projection.png'), cmap='gray')
# axes[3].imshow(imageio.imread('./final/1110070263_1110190110_suite2p_maximum_projection.png'), cmap='gray')

# axes[4].plot(row_one_old['x'], alpha=0.5, color='blue', label='smooth_sigma_time = 0.0')
# axes[4].plot(row_one['x'], alpha=0.5, color='red', label='smooth_sigma_time = 3.0')
# axes[4].set_ylabel('x', size=30)
# axes[4].set_xlabel('Frame Number', size=20)
# axes[4].tick_params(axis='y', labelsize=30)

# axes[5].plot(row_one_old['y'], alpha=0.5, color='blue', label='smooth_sigma_time = 0.0')
# axes[5].plot(row_one['y'], alpha=0.5, color='red', label='smooth_sigma_time = 3.0')
# axes[5].set_ylabel('y', size=30)
# axes[5].set_xlabel('Frame Number', size=20)
# axes[5].tick_params(axis='y', labelsize=30)

# axes[5].legend(loc='upper right', prop={'size': 20})

# axes[0].set_title("Default Avg. Proj.", size=40)
# axes[1].set_title("Default Max. Proj.", size=40)
# axes[2].set_title("Avg. Proj.", size=40)
# axes[3].set_title("Max. Proj.", size=40)
# axes[4].set_title("Session 1110070263 Experiment 1110190110", size=40)

# #######################################################
# #      Row 2
# row_one = pd.read_csv("./final/1110609895_1110738903_suite2p_rigid_motion_transform.csv")
# row_one_old = pd.read_csv("./final/1110738903_suite2p_rigid_motion_transform.csv")

# axes[6].imshow(imageio.imread('./final/1110738903_suite2p_average_projection.png'), cmap='gray')
# axes[7].imshow(imageio.imread('./final/1110738903_suite2p_maximum_projection.png'), cmap='gray')

# axes[8].imshow(imageio.imread('./final/1110609895_1110738903_suite2p_average_projection.png'), cmap='gray')
# axes[9].imshow(imageio.imread('./final/1110609895_1110738903_suite2p_maximum_projection.png'), cmap='gray')

# axes[10].plot(row_one_old['x'], alpha=0.5, color='blue', label='smooth_sigma_time = 0.0')
# axes[10].plot(row_one['x'], alpha=0.5, color='red', label='smooth_sigma_time = 3.0')
# axes[10].set_ylabel('x', size=30)
# axes[10].set_xlabel('Frame Number', size=20)
# axes[10].tick_params(axis='y', labelsize=30)

# axes[11].plot(row_one_old['y'], alpha=0.5, color='blue', label='smooth_sigma_time = 0.0')
# axes[11].plot(row_one['y'], alpha=0.5, color='red', label='smooth_sigma_time = 3.0')
# axes[11].set_ylabel('y', size=30)
# axes[11].set_xlabel('Frame Number', size=20)
# axes[11].tick_params(axis='y', labelsize=30)

# axes[11].legend(loc='upper right', prop={'size': 20})

# axes[10].set_title("Session 1110609895 Experiment 1110738903", size=40)

# #######################################################
# #      Row 3
# row_one = pd.read_csv("./final/1110835610_1110928086_suite2p_rigid_motion_transform.csv")
# row_one_old = pd.read_csv("./final/1110928086_suite2p_rigid_motion_transform.csv")

# axes[12].imshow(imageio.imread('./final/1110928086_suite2p_average_projection.png'), cmap='gray')
# axes[13].imshow(imageio.imread('./final/1110928086_suite2p_maximum_projection.png'), cmap='gray')

# axes[14].imshow(imageio.imread('./final/1110835610_1110928086_suite2p_average_projection.png'), cmap='gray')
# axes[15].imshow(imageio.imread('./final/1110835610_1110928086_suite2p_maximum_projection.png'), cmap='gray')

# axes[16].plot(row_one_old['x'], alpha=0.5, color='blue', label='smooth_sigma_time = 0.0')
# axes[16].plot(row_one['x'], alpha=0.5, color='red', label='smooth_sigma_time = 3.0')
# axes[16].set_ylabel('x', size=30)
# axes[16].set_xlabel('Frame Number', size=20)
# axes[16].tick_params(axis='y', labelsize=30)

# axes[17].plot(row_one_old['y'], alpha=0.5, color='blue', label='smooth_sigma_time = 0.0')
# axes[17].plot(row_one['y'], alpha=0.5, color='red', label='smooth_sigma_time = 3.0')
# axes[17].set_ylabel('y', size=30)
# axes[17].set_xlabel('Frame Number', size=20)
# axes[17].tick_params(axis='y', labelsize=30)

# axes[17].legend(loc='upper right', prop={'size': 20})

# axes[16].set_title("Session 1110835610 Experiment 1110928086", size=40)

# #######################################################
# #      Row 4
# row_one = pd.read_csv("./final/1111018150_1111121124_suite2p_rigid_motion_transform.csv")
# row_one_old = pd.read_csv("./final/1111121124_suite2p_rigid_motion_transform.csv")

# axes[18].imshow(imageio.imread('./final/1111121124_suite2p_average_projection.png'), cmap='gray')
# axes[19].imshow(imageio.imread('./final/1111121124_suite2p_maximum_projection.png'), cmap='gray')

# axes[20].imshow(imageio.imread('./final/1111018150_1111121124_suite2p_average_projection.png'), cmap='gray')
# axes[21].imshow(imageio.imread('./final/1111018150_1111121124_suite2p_maximum_projection.png'), cmap='gray')

# axes[22].plot(row_one_old['x'], alpha=0.5, color='blue', label='smooth_sigma_time = 0.0')
# axes[22].plot(row_one['x'], alpha=0.5, color='red', label='smooth_sigma_time = 3.0')
# axes[22].set_ylabel('x', size=30)
# axes[22].set_xlabel('Frame Number', size=20)
# axes[22].tick_params(axis='y', labelsize=30)

# axes[23].plot(row_one_old['y'], alpha=0.5, color='blue', label='smooth_sigma_time = 0.0')
# axes[23].plot(row_one['y'], alpha=0.5, color='red', label='smooth_sigma_time = 3.0')
# axes[23].set_ylabel('y', size=30)
# axes[23].set_xlabel('Frame Number', size=20)
# axes[23].tick_params(axis='y', labelsize=30)

# axes[23].legend(loc='upper right', prop={'size': 20})

# axes[22].set_title("Session 1111018150 Experiment 1111121124", size=40)

# #######################################################
# #      Row 5
# row_one = pd.read_csv("./final/1111197126_1111313373_suite2p_rigid_motion_transform.csv")
# row_one_old = pd.read_csv("./final/1111313373_suite2p_rigid_motion_transform.csv")

# axes[24].imshow(imageio.imread('./final/1111313373_suite2p_average_projection.png'), cmap='gray')
# axes[25].imshow(imageio.imread('./final/1111313373_suite2p_maximum_projection.png'), cmap='gray')

# axes[26].imshow(imageio.imread('./final/1111197126_1111313373_suite2p_average_projection.png'), cmap='gray')
# axes[27].imshow(imageio.imread('./final/1111197126_1111313373_suite2p_maximum_projection.png'), cmap='gray')

# axes[28].plot(row_one_old['x'], alpha=0.5, color='blue', label='smooth_sigma_time = 0.0')
# axes[28].plot(row_one['x'], alpha=0.5, color='red', label='smooth_sigma_time = 3.0')
# axes[28].set_ylabel('x', size=30)
# axes[28].set_xlabel('Frame Number', size=20)
# axes[28].tick_params(axis='y', labelsize=30)

# axes[29].plot(row_one_old['y'], alpha=0.5, color='blue', label='smooth_sigma_time = 0.0')
# axes[29].plot(row_one['y'], alpha=0.5, color='red', label='smooth_sigma_time = 3.0')
# axes[29].set_ylabel('y', size=30)
# axes[29].set_xlabel('Frame Number', size=20)
# axes[29].tick_params(axis='y', labelsize=30)

# axes[29].legend(loc='upper right', prop={'size': 20})

# axes[28].set_title("Session 1111197126 Experiment 1111313373", size=40)


# #######################################################
# #      Row 6
# row_one = pd.read_csv("./final/1111398401_1111493062_suite2p_rigid_motion_transform.csv")
# row_one_old = pd.read_csv("./final/1111493062_suite2p_rigid_motion_transform.csv")

# axes[30].imshow(imageio.imread('./final/1111493062_suite2p_average_projection.png'), cmap='gray')
# axes[31].imshow(imageio.imread('./final/1111493062_suite2p_maximum_projection.png'), cmap='gray')

# axes[32].imshow(imageio.imread('./final/1111398401_1111493062_suite2p_average_projection.png'), cmap='gray')
# axes[33].imshow(imageio.imread('./final/1111398401_1111493062_suite2p_maximum_projection.png'), cmap='gray')

# axes[34].plot(row_one_old['x'], alpha=0.5, color='blue', label='smooth_sigma_time = 0.0')
# axes[34].plot(row_one['x'], alpha=0.5, color='red', label='smooth_sigma_time = 3.0')
# axes[34].set_ylabel('x', size=30)
# axes[34].set_xlabel('Frame Number', size=20)
# axes[34].tick_params(axis='y', labelsize=30)

# axes[35].plot(row_one_old['y'], alpha=0.5, color='blue', label='smooth_sigma_time = 0.0')
# axes[35].plot(row_one['y'], alpha=0.5, color='red', label='smooth_sigma_time = 3.0')
# axes[35].set_ylabel('y', size=30)
# axes[35].set_xlabel('Frame Number', size=20)
# axes[35].tick_params(axis='y', labelsize=30)

# axes[35].legend(loc='upper right', prop={'size': 20})

# axes[34].set_title("Session 1111398401 Experiment 1111493062", size=40)

# #######################################################
# #      Save

# fig.tight_layout()
# fig.savefig('./final/multiple_bad_sst_3.png')