API tools faq
paste
Advertisement
Guest User

Untitled

a guest
Jul 16th, 2019
109
0
Never
Add comment
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Python 20.12 KB | None | 0 0
raw download clone embed print report
  1. def Plot_t_delta_incl(d_delta,d_delta_av,d_delta_comp):
  2.        
  3.     for n_centralina in d_delta:
  4.         # find temperature column name and tiltimeter columns names
  5.         l_col_incl=[]
  6.         l_n_incl=[]
  7.         for column_name in d_delta[n_centralina].columns:
  8.             if "TEMP"==column_name[:4]:
  9.                 T_col_name=column_name
  10.             elif "INCL"==column_name[:4]:
  11.                 l_col_incl.append(column_name)
  12.                 n_incl=int(''.join(x for x in column_name if x.isdigit()))
  13.                 if n_incl not in l_n_incl:
  14.                     l_n_incl.append(n_incl)
  15.                    
  16.         if l_n_incl:               # check if list is not empty
  17.             for n_incl in l_n_incl:
  18.                 l_data=[name for name in l_col_incl if str(n_incl) in name]
  19.                 l_data_pila=[name for name in l_data if "PILA" in name]
  20.                 l_data_spalla=[name for name in l_data if "SPALLA" in name]
  21.  
  22.                 if l_data_pila:
  23.                     t=d_delta[n_centralina][l_data_pila[0]].index
  24.                     t_av=d_delta_av[n_centralina][l_data_pila[0]].index
  25.                     data1 = d_delta[n_centralina][l_data_pila[0]]
  26.                     data2 = d_delta[n_centralina][l_data_pila[1]]
  27.                     data3=d_delta_av[n_centralina][l_data_pila[0]]
  28.                     data4=d_delta_av[n_centralina][l_data_pila[1]]
  29.                     #data5=d_delta_comp[n_centralina][l_data_pila[0]]
  30.                     #data6=d_delta_comp[n_centralina][l_data_pila[1]]
  31.                     dataT = d_delta[n_centralina][T_col_name].values
  32.                
  33.                     fig, ax1 = plt.subplots()
  34.                     color0 = 'tab:orange'
  35.                     color1 = 'tab:blue'
  36.                     ax1.set_xlabel('time')
  37.                     ax1.set_ylabel('rotazione [°]')
  38.                     ax1.plot(t, data1, color=color0, label=l_data_pila[0],linewidth=0.5)
  39.                     ax1.plot(t, data2, color=color1, label=l_data_pila[1],linewidth=0.5)
  40.                     ax1.plot(t_av, data3, color=color0, label=l_data_pila[0]+'av',linewidth=1.0)
  41.                     ax1.plot(t_av, data4, color=color1, label=l_data_pila[1]+'av',linewidth=1.0)
  42.                     #ax1.plot(t, data5, color=color0, label=l_data_pila[0]+'comp',linewidth=1.0,linestyle="--")
  43.                     #ax1.plot(t, data6, color=color1, label=l_data_pila[1]+'comp',linewidth=1.0,linestyle="--")
  44.                     ax1.tick_params(axis='y')
  45.  
  46.                     # set y axis limits            
  47.                     ax1.set_ylim([-0.2,0.2])
  48.                    
  49.                     plt.grid(True, which='major',axis='both', linestyle='--',dashes=[10, 10], linewidth=0.5)
  50.                     leg=plt.legend(bbox_to_anchor=(0., 1.02, 1., .102), loc=3,
  51.                                    ncol=2, mode="expand", borderaxespad=0.,fontsize=8)
  52.                
  53.                     for line in leg.get_lines():
  54.                         line.set_linewidth(2)
  55.            
  56.                     ax2 = ax1.twinx()  # instantiate a second axes that shares the same x-axis
  57.            
  58.                     color = 'tab:grey'
  59.                     ax2.set_ylabel('T [°C]', color=color)
  60.                     ax2.plot(t, dataT, color=color, label="TEMPERATURE", linewidth=0.5,alpha=0.4)
  61.                     ax2.tick_params(axis='y', labelcolor=color)
  62.                     ax2.set_ylim([-25,25])
  63.                    
  64.                     fig.tight_layout()
  65.                     fig.canvas.set_window_title(l_data_pila[0][:-1])
  66.                     plt.show()
  67.                    
  68.                 if l_data_spalla:
  69.                     t=d_delta[n_centralina][l_data_spalla[0]].index
  70.                     t_av=d_delta_av[n_centralina][l_data_spalla[0]].index
  71.                     data1 = d_delta[n_centralina][l_data_spalla[0]]
  72.                     data2 = d_delta[n_centralina][l_data_spalla[1]]
  73.                     data3=d_delta_av[n_centralina][l_data_spalla[0]]
  74.                     data4=d_delta_av[n_centralina][l_data_spalla[1]]
  75.                     #data5=d_delta_comp[n_centralina][l_data_spalla[0]]
  76.                     #data6=d_delta_comp[n_centralina][l_data_spalla[1]]
  77.                     dataT = d_delta[n_centralina][T_col_name].values
  78.                
  79.                     fig, ax1 = plt.subplots()
  80.                     color0 = 'tab:orange'
  81.                     color1 = 'tab:blue'
  82.                     ax1.set_xlabel('time')
  83.                     ax1.set_ylabel('rotazione [°]')
  84.                     ax1.plot(t, data1, color=color0, label=l_data_spalla[0],linewidth=0.5)
  85.                     ax1.plot(t, data2, color=color1, label=l_data_spalla[1],linewidth=0.5)
  86.                     ax1.plot(t_av, data3, color=color0, label=l_data_spalla[0]+'av',linewidth=1.0)
  87.                     ax1.plot(t_av, data4, color=color1, label=l_data_spalla[1]+'av',linewidth=1.0)
  88.                     #ax1.plot(t, data5, color=color0, label=l_data_spalla[0]+'comp',linewidth=1.0,linestyle="--")
  89.                     #ax1.plot(t, data6, color=color1, label=l_data_spalla[1]+'comp',linewidth=1.0,linestyle="--")
  90.                     ax1.tick_params(axis='y')
  91.  
  92.                     # set y axis limits            
  93.                     ax1.set_ylim([-0.2,0.2])
  94.                    
  95.                     plt.grid(True, which='major',axis='both', linestyle='--',dashes=[10, 10], linewidth=0.5)
  96.                     leg=plt.legend(bbox_to_anchor=(0., 1.02, 1., .102), loc=3,
  97.                                    ncol=2, mode="expand", borderaxespad=0.,fontsize=8)
  98.                
  99.                     for line in leg.get_lines():
  100.                         line.set_linewidth(2)
  101.            
  102.                     ax2 = ax1.twinx()  # instantiate a second axes that shares the same x-axis
  103.            
  104.                     color = 'tab:grey'
  105.                     ax2.set_ylabel('T [°C]', color=color)
  106.                     ax2.plot(t, dataT, color=color, label="TEMPERATURE", linewidth=0.5,alpha=0.4)
  107.                     ax2.tick_params(axis='y', labelcolor=color)
  108.                     ax2.set_ylim([-25,25])
  109.                    
  110.                     fig.tight_layout()
  111.                     fig.canvas.set_window_title(l_data_spalla[0][:-1])
  112.                     plt.show()
  113.                    
  114. def Plot_t_delta_comp_incl(d_delta_comp, d_delta_av_comp):
  115.        
  116.     for n_centralina in d_delta_comp:
  117.         # find temperature column name and tiltimeter columns names
  118.         l_col_incl=[]
  119.         l_n_incl=[]
  120.         for column_name in d_delta_comp[n_centralina].columns:
  121.             if "TEMP"==column_name[:4]:
  122.                 T_col_name=column_name
  123.             elif "INCL"==column_name[:4]:
  124.                 l_col_incl.append(column_name)
  125.                 n_incl=int(''.join(x for x in column_name if x.isdigit()))
  126.                 if n_incl not in l_n_incl:
  127.                     l_n_incl.append(n_incl)
  128.                    
  129.         if l_n_incl:                    # check if list is not empty
  130.             for n_incl in l_n_incl:
  131.                 l_data=[name for name in l_col_incl if str(n_incl) in name]
  132.                 l_data_pila=[name for name in l_data if "PILA" in name]
  133.                 l_data_spalla=[name for name in l_data if "SPALLA" in name]
  134.  
  135.                 if l_data_pila:
  136.                     t=d_delta_comp[n_centralina][l_data_pila[0]].index
  137.                     t_av=d_delta_av_comp[n_centralina][l_data_pila[0]].index
  138.                     data1 = d_delta_comp[n_centralina][l_data_pila[0]]
  139.                     data2 = d_delta_comp[n_centralina][l_data_pila[1]]
  140.                     data3=d_delta_av_comp[n_centralina][l_data_pila[0]]
  141.                     data4=d_delta_av_comp[n_centralina][l_data_pila[1]]
  142.                     dataT = d_delta_comp[n_centralina][T_col_name].values
  143.                
  144.                     fig, ax1 = plt.subplots()
  145.                     color0 = 'tab:orange'
  146.                     color1 = 'tab:blue'
  147.                     ax1.set_xlabel('time')
  148.                     ax1.set_ylabel('rotazione [°]')
  149.                     ax1.plot(t, data1, color=color0, label=l_data_pila[0]+'comp',linewidth=0.5)
  150.                     ax1.plot(t, data2, color=color1, label=l_data_pila[1]+'comp',linewidth=0.5)
  151.                     ax1.plot(t_av, data3, color=color0, label=l_data_pila[0]+'av'+'comp',linewidth=1.0)
  152.                     ax1.plot(t_av, data4, color=color1, label=l_data_pila[1]+'av'+'comp',linewidth=1.0)
  153.                     ax1.tick_params(axis='y')
  154.  
  155.                     # set y axis limits            
  156.                     ax1.set_ylim([-0.2,0.2])
  157.                    
  158.                     plt.grid(True, which='major',axis='both', linestyle='--',dashes=[10, 10], linewidth=0.5)
  159.                     leg=plt.legend(bbox_to_anchor=(0., 1.02, 1., .102), loc=3,
  160.                                    ncol=2, mode="expand", borderaxespad=0.,fontsize=8)
  161.                
  162.                     for line in leg.get_lines():
  163.                         line.set_linewidth(2)
  164.                     '''
  165.                    ax2 = ax1.twinx()  # instantiate a second axes that shares the same x-axis
  166.                    
  167.                    color = 'tab:grey'
  168.                    ax2.set_ylabel('T [°C]', color=color)
  169.                    ax2.plot(t, dataT, color=color, label="TEMPERATURE", linewidth=0.5,alpha=0.4)
  170.                    ax2.tick_params(axis='y', labelcolor=color)
  171.                    ax2.set_ylim([-25,25])
  172.                    '''
  173.                     fig.tight_layout()
  174.                     fig.canvas.set_window_title(l_data_pila[0][:-1]+"comp")
  175.                     plt.show()
  176.                    
  177.                 if l_data_spalla:
  178.                     t=d_delta_comp[n_centralina][l_data_spalla[0]].index
  179.                     t_av=d_delta_av_comp[n_centralina][l_data_spalla[0]].index
  180.                     data1 = d_delta_comp[n_centralina][l_data_spalla[0]]
  181.                     data2 = d_delta_comp[n_centralina][l_data_spalla[1]]
  182.                     data3=d_delta_av_comp[n_centralina][l_data_spalla[0]]
  183.                     data4=d_delta_av_comp[n_centralina][l_data_spalla[1]]
  184.                     dataT = d_delta_comp[n_centralina][T_col_name].values
  185.                
  186.                     fig, ax1 = plt.subplots()
  187.                     color0 = 'tab:orange'
  188.                     color1 = 'tab:blue'
  189.                     ax1.set_xlabel('time')
  190.                     ax1.set_ylabel('rotazione [°]')
  191.                     ax1.plot(t, data1, color=color0, label=l_data_spalla[0]+'comp',linewidth=0.5)
  192.                     ax1.plot(t, data2, color=color1, label=l_data_spalla[1]+'comp',linewidth=0.5)
  193.                     ax1.plot(t_av, data3, color=color0, label=l_data_spalla[0]+'av'+'comp',linewidth=1.0)
  194.                     ax1.plot(t_av, data4, color=color1, label=l_data_spalla[1]+'av'+'comp',linewidth=1.0)
  195.                     ax1.tick_params(axis='y')
  196.  
  197.                     # set y axis limits            
  198.                     ax1.set_ylim([-0.2,0.2])
  199.                    
  200.                     plt.grid(True, which='major',axis='both', linestyle='--',dashes=[10, 10], linewidth=0.5)
  201.                     leg=plt.legend(bbox_to_anchor=(0., 1.02, 1., .102), loc=3,
  202.                                    ncol=2, mode="expand", borderaxespad=0.,fontsize=8)
  203.                
  204.                     for line in leg.get_lines():
  205.                         line.set_linewidth(2)
  206.                     '''
  207.                    ax2 = ax1.twinx()  # instantiate a second axes that shares the same x-axis
  208.            
  209.                    color = 'tab:grey'
  210.                    ax2.set_ylabel('T [°C]', color=color)
  211.                    ax2.plot(t, dataT, color=color, label="TEMPERATURE", linewidth=0.5,alpha=0.4)
  212.                    ax2.tick_params(axis='y', labelcolor=color)
  213.                    ax2.set_ylim([-25,25])
  214.                    '''
  215.                     fig.canvas.set_window_title(l_data_spalla[0][:-1]+"comp")
  216.                     fig.tight_layout()    
  217.                     plt.show()
  218.                                      
  219.  
  220. def Plot_t_delta_fess(d_delta,d_delta_av,d_delta_comp,d_coeff_Pear):
  221.  
  222.     for n_centralina in d_delta:        
  223.         # find temperature column name and crack width gauge columns names
  224.         l_col_fess=[]
  225.         l_n_fess=[]
  226.         for column_name in d_delta[n_centralina].columns:
  227.             if "TEMP"==column_name[:4]:
  228.                 T_col_name=column_name
  229.             elif "FESS"==column_name[:4]:
  230.                 l_col_fess.append(column_name)
  231.                 n_fess=int(''.join(x for x in column_name if x.isdigit()))
  232.                 if n_fess not in l_n_fess:
  233.                     l_n_fess.append(n_fess)
  234.         if l_n_fess:     # check if list is not empty
  235.             for n_fess in l_n_fess:
  236.                 if n_fess != 100:
  237.                     l_data=[name for name in l_col_fess if str(n_fess) in name]
  238.                     t = d_delta[n_centralina][l_data[0]].index
  239.                     t_av=d_delta_av[n_centralina][l_data[0]].index
  240.                     data1 = d_delta[n_centralina][l_data[0]]
  241.                     data2 = d_delta[n_centralina][l_data[1]]
  242.                     data3 = d_delta[n_centralina][l_data[2]]
  243.                     data4 = d_delta_av[n_centralina][l_data[0]]
  244.                     data5 = d_delta_av[n_centralina][l_data[1]]
  245.                     data6 = d_delta_av[n_centralina][l_data[2]]
  246.                     dataT = d_delta[n_centralina][T_col_name].values
  247.              
  248.                     fig, ax1 = plt.subplots()
  249.  
  250.                     color = 'tab:orange'
  251.                     ax1.set_xlabel('time')
  252.                     ax1.set_ylabel('spostamento [mm]', color=color)
  253.                     ax1.plot(t, data1, color=color, label=l_data[0],linewidth=0.5)
  254.                     ax1.plot(t_av, data4, color=color, label=l_data[0]+'av',linewidth=0.5)
  255.                     ax1.tick_params(axis='y', labelcolor=color)
  256.  
  257.                     color = 'tab:blue'
  258.                     ax1.set_xlabel('time')
  259.                     ax1.set_ylabel('spostamento [mm]', color=color)
  260.                     ax1.plot(t, data2, color=color, label=l_data[1],linewidth=0.5)
  261.                     ax1.plot(t_av, data5, color=color, label=l_data[1]+'av',linewidth=0.5)
  262.                     ax1.tick_params(axis='y', labelcolor=color)
  263.            
  264.                     color = 'tab:green'
  265.                     ax1.set_xlabel('time')
  266.                     ax1.set_ylabel('spostamento [mm]', color=color)
  267.                     ax1.plot(t, data3, color=color, label=l_data[2],linewidth=0.5)
  268.                     ax1.plot(t_av, data6, color=color, label=l_data[2]+'av',linewidth=0.5)
  269.                     ax1.tick_params(axis='y', labelcolor=color)
  270.            
  271.                     ax1.set_ylim([-15,15])
  272.                    
  273.                     textstr1=l_data[0]+' '+'c_Pear='+str(round(d_coeff_Pear[n_centralina][l_data[0]][0],2))
  274.                     textstr2=l_data[1]+' '+'c_Pear='+str(round(d_coeff_Pear[n_centralina][l_data[1]][0],2))
  275.                     textstr3=l_data[2]+' '+'c_Pear='+str(round(d_coeff_Pear[n_centralina][l_data[2]][0],2))
  276.                     ax1.text(0.05, 0.95,'%s\n%s\n%s'%(textstr1,textstr2,textstr3),
  277.                              transform=ax1.transAxes,fontsize=6,
  278.                              verticalalignment='top')
  279.        
  280.                     plt.grid(True, which='major',axis='both', linestyle='--',dashes=[10, 10], linewidth=0.5)
  281.                     leg=plt.legend(bbox_to_anchor=(0., 1.02, 1., .102), loc=3,
  282.                                ncol=3, mode="expand", borderaxespad=0.,fontsize=8)
  283.                    
  284.                     for line in leg.get_lines():
  285.                         line.set_linewidth(2)
  286.        
  287.                     ax2 = ax1.twinx()  # instantiate a second axes that shares the same x-axis
  288.  
  289.                     color = 'tab:grey'
  290.                     ax2.set_ylabel('T [°C]', color=color)
  291.                     ax2.plot(t, dataT, color=color, label="TEMPERATURE", linewidth=0.5, alpha=0.4)
  292.                     ax2.tick_params(axis='y', labelcolor=color)
  293.                     ax2.set_ylim([-25,25])
  294.  
  295.                     fig.canvas.set_window_title(l_data[0][:-2]+l_data[0][-1])
  296.                     fig.tight_layout()    
  297.                     plt.show()
  298.  
  299. def Plot_T_delta(d_delta,d_delta_av,d_fit_par):
  300.    
  301.     for n_centralina in d_delta:
  302.        
  303.         # find temperature column name and tiltimeter columns names
  304.         l_col_incl=[]
  305.         l_n_incl=[]
  306.         for column_name in d_delta[n_centralina].columns:
  307.             if "TEMP"==column_name[:4]:
  308.                 T_col_name=column_name
  309.             elif "INCL"==column_name[:4]:
  310.                 l_col_incl.append(column_name)
  311.                 n_incl=int(''.join(x for x in column_name if x.isdigit()))
  312.                 if n_incl not in l_n_incl:
  313.                     l_n_incl.append(n_incl)
  314.                    
  315.         if l_n_incl:               # check if list is not empty
  316.             for n_incl in l_n_incl:
  317.                 l_data=[name for name in l_col_incl if str(n_incl) in name]
  318.                 l_data_pila=[name for name in l_data if "PILA" in name]
  319.                 l_data_spalla=[name for name in l_data if "SPALLA" in name]
  320.  
  321.                 if l_data_pila:
  322.                     T = d_delta[n_centralina][T_col_name]
  323.                     data1 = d_delta[n_centralina][l_data_pila[0]]
  324.                     data2= d_delta[n_centralina][l_data_pila[1]]
  325.    
  326.                     p1 = np.poly1d(d_fit_par[n_centralina][l_data_pila[0]])
  327.                     p2 = np.poly1d(d_fit_par[n_centralina][l_data_pila[1]])
  328.                     xp = np.linspace(min(T), max(T), 100)
  329.  
  330.                     fig, ax1 = plt.subplots()
  331.    
  332.                     color0 = 'tab:orange'
  333.                     color1 = 'tab:blue'
  334.                     ax1.set_xlabel('TEMPERATURE')
  335.                     ax1.set_ylabel('rotazione [°]', color=color0)
  336.                     ax1.scatter(T, data1, s=1, color=color0, label=l_data_pila[0])
  337.                     ax1.scatter(T, data2, s=1, color=color1, label=l_data_pila[1])
  338.                     ax1.plot(xp,p1(xp),color=color0)
  339.                     ax1.plot(xp,p2(xp), color=color1)
  340.                     ax1.tick_params(axis='y')
  341.                     #ax1.set_ylim()
  342.    
  343.                     plt.grid(True, which='major',axis='both', linestyle='--',dashes=[10, 10], linewidth=0.5)
  344.                     plt.legend(bbox_to_anchor=(0., 1.02, 1., .102), loc=3,
  345.                                ncol=2, mode="expand", borderaxespad=0.)
  346.                    
  347.                     fig.canvas.set_window_title(l_data_pila[0][:-1]+"T-delta")
  348.                     fig.tight_layout()      # otherwise the right y-label is slightly clipped
  349.                     plt.show()
  350.                    
  351.                 if l_data_spalla:
  352.                     T = d_delta[n_centralina][T_col_name]
  353.                     data1 = d_delta[n_centralina][l_data_spalla[0]]
  354.                     data2= d_delta[n_centralina][l_data_spalla[1]]
  355.    
  356.                     p1 = np.poly1d(d_fit_par[n_centralina][l_data_spalla[0]])
  357.                     p2 = np.poly1d(d_fit_par[n_centralina][l_data_spalla[1]])
  358.                     xp = np.linspace(min(T), max(T), 100)
  359.  
  360.                     fig, ax1 = plt.subplots()
  361.    
  362.                     color0 = 'tab:orange'
  363.                     color1 = 'tab:blue'
  364.                     ax1.set_xlabel('TEMPERATURE')
  365.                     ax1.set_ylabel('rotazione [°]', color=color0)
  366.                     ax1.scatter(T, data1, s=1, color=color0, label=l_data_spalla[0])
  367.                     ax1.scatter(T, data2, s=1, color=color1, label=l_data_spalla[1])
  368.                     ax1.plot(xp,p1(xp),color=color0)
  369.                     ax1.plot(xp,p2(xp), color=color1)
  370.                     ax1.tick_params(axis='y')
  371.                     #ax1.set_ylim()
  372.    
  373.                     plt.grid(True, which='major',axis='both', linestyle='--',dashes=[10, 10], linewidth=0.5)
  374.                     plt.legend(bbox_to_anchor=(0., 1.02, 1., .102), loc=3,
  375.                                ncol=2, mode="expand", borderaxespad=0.)
  376.                    
  377.                     fig.canvas.set_window_title(l_data_spalla[0][:-1]+" T-delta")
  378.                     fig.tight_layout()      # otherwise the right y-label is slightly clipped
  379.                     plt.show()
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement
Public Pastes
Advertisement
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy.  OK, I Understand
Not a member of Pastebin yet?
Sign Up, it unlocks many cool features!