source: trunk/src/tup_moy_par_zone.py @ 54

#amsub
ialt34_abch1=(ambch1[13,:]>=3000)&(ambch1[13,:]<=4000)
ialt34_abch2=(ambch2[13,:]>=3000)&(ambch2[13,:]<=4000)
ialt34_abch3=(ambch3[13,:]>=3000)&(ambch3[13,:]<=4000)
ialt34_abch4=(ambch4[13,:]>=3000)&(ambch4[13,:]<=4000)
ialt34_abch5=(ambch5[13,:]>=3000)&(ambch5[13,:]<=4000)

tup_ab1_alt34=mean(ambch1[18,ialt34_abch1])
tup_ab2_alt34=mean(ambch2[18,ialt34_abch2])
tup_ab3_alt34=mean(ambch3[18,ialt34_abch3])
tup_ab4_alt34=mean(ambch4[18,ialt34_abch4])
tup_ab5_alt34=mean(ambch5[18,ialt34_abch5])

#amsub angles faibles
ialt34_abch1n1=ialt34_abch1 & (ambch1[7,:]<=13) | ialt34_abch1 & (ambch1[7,:]>=78)
ialt34_abch2n1=ialt34_abch2 & (ambch2[7,:]<=13) | ialt34_abch2 & (ambch2[7,:]>=78)
ialt34_abch3n1=ialt34_abch3 & (ambch3[7,:]<=13) | ialt34_abch3 & (ambch3[7,:]>=78)
ialt34_abch4n1=ialt34_abch4 & (ambch4[7,:]<=13) | ialt34_abch4 & (ambch4[7,:]>=78)
ialt34_abch5n1=ialt34_abch5 & (ambch5[7,:]<=13) | ialt34_abch5 & (ambch5[7,:]>=78)

tup_ab_alt34n1=np.zeros([5],float)
tup_ab_alt34n1[0]=mean(ambch1[18,ialt34_abch1n1])
tup_ab_alt34n1[1]=mean(ambch2[18,ialt34_abch2n1])
tup_ab_alt34n1[2]=mean(ambch3[18,ialt34_abch3n1])
tup_ab_alt34n1[3]=mean(ambch4[18,ialt34_abch4n1])
tup_ab_alt34n1[4]=mean(ambch5[18,ialt34_abch5n1])

#ecart type
tupet_ab_alt34n1=np.zeros([5],float)
tupet_ab_alt34n1[0]=sqrt(sum(ambch1[18,ialt34_abch1n1]**2)/size(ambch1[18,ialt34_abch1n1])-mean(ambch1[18,ialt34_abch1n1])**2)
tupet_ab_alt34n1[1]=sqrt(sum(ambch2[18,ialt34_abch2n1]**2)/size(ambch2[18,ialt34_abch2n1])-mean(ambch2[18,ialt34_abch2n1])**2)
tupet_ab_alt34n1[2]=sqrt(sum(ambch3[18,ialt34_abch3n1]**2)/size(ambch3[18,ialt34_abch3n1])-mean(ambch3[18,ialt34_abch3n1])**2)
tupet_ab_alt34n1[3]=sqrt(sum(ambch4[18,ialt34_abch4n1]**2)/size(ambch4[18,ialt34_abch4n1])-mean(ambch4[18,ialt34_abch4n1])**2)
tupet_ab_alt34n1[4]=sqrt(sum(ambch5[18,ialt34_abch5n1]**2)/size(ambch5[18,ialt34_abch5n1])-mean(ambch5[18,ialt34_abch5n1])**2)


#amsua
ialt34_ach1=(amch1[13,:]>=3000)&(amch1[13,:]<=4000)
ialt34_ach2=(amch2[13,:]>=3000)&(amch2[13,:]<=4000)
ialt34_ach3=(amch3[13,:]>=3000)&(amch3[13,:]<=4000)
ialt34_ach4=(amch4[13,:]>=3000)&(amch4[13,:]<=4000)
ialt34_ach5=(amch5[13,:]>=3000)&(amch5[13,:]<=4000)
ialt34_ach6=(amch6[13,:]>=3000)&(amch6[13,:]<=4000)
ialt34_ach7=(amch7[13,:]>=3000)&(amch7[13,:]<=4000)
ialt34_ach15=(amch15[13,:]>=3000)&(amch15[13,:]<=4000)

tup_a1_alt34=mean(amch1[18,ialt34_ach1])
tup_a2_alt34=mean(amch2[18,ialt34_ach2])
tup_a3_alt34=mean(amch3[18,ialt34_ach3])
tup_a4_alt34=mean(amch4[18,ialt34_ach4])
tup_a5_alt34=mean(amch5[18,ialt34_ach5])
tup_a6_alt34=mean(amch6[18,ialt34_ach6])
tup_a7_alt34=mean(amch7[18,ialt34_ach7])
tup_a15_alt34=mean(amch15[18,ialt34_ach15])

#amsua angles faibles
ialt34_ach1n1=ialt34_ach1 & (amch1[7,:]<=5) | ialt34_ach1 & (amch1[7,:]>=26)
ialt34_ach2n1=ialt34_ach2 & (amch2[7,:]<=5) | ialt34_ach2 & (amch2[7,:]>=26)
ialt34_ach3n1=ialt34_ach3 & (amch3[7,:]<=5) | ialt34_ach3 & (amch3[7,:]>=26)
ialt34_ach4n1=ialt34_ach4 & (amch4[7,:]<=5) | ialt34_ach4 & (amch4[7,:]>=26)
ialt34_ach5n1=ialt34_ach5 & (amch5[7,:]<=5) | ialt34_ach5 & (amch5[7,:]>=26)
ialt34_ach6n1=ialt34_ach6 & (amch6[7,:]<=5) | ialt34_ach6 & (amch6[7,:]>=26)
ialt34_ach7n1=ialt34_ach7 & (amch7[7,:]<=5) | ialt34_ach7 & (amch7[7,:]>=26)
ialt34_ach15n1=ialt34_ach15 & (amch15[7,:]<=5) | ialt34_ach15 & (amch15[7,:]>=26)

tup_aa_alt34n1=np.zeros([8],float)
tup_aa_alt34n1[0]=mean(amch1[18,ialt34_ach1n1])
tup_aa_alt34n1[1]=mean(amch2[18,ialt34_ach2n1])
tup_aa_alt34n1[2]=mean(amch3[18,ialt34_ach3n1])
tup_aa_alt34n1[3]=mean(amch4[18,ialt34_ach4n1])
tup_aa_alt34n1[4]=mean(amch5[18,ialt34_ach5n1])
tup_aa_alt34n1[5]=mean(amch6[18,ialt34_ach6n1])
tup_aa_alt34n1[6]=mean(amch7[18,ialt34_ach7n1])
tup_aa_alt34n1[7]=mean(amch15[18,ialt34_ach15n1])

#ecart type
tupet_aa_alt34n1=np.zeros([8],float)
tupet_aa_alt34n1[0]=sqrt(sum(amch1[18,ialt34_ach1n1]**2)/size(amch1[18,ialt34_ach1n1])-mean(amch1[18,ialt34_ach1n1])**2)
tupet_aa_alt34n1[1]=sqrt(sum(amch2[18,ialt34_ach2n1]**2)/size(amch2[18,ialt34_ach2n1])-mean(amch2[18,ialt34_ach2n1])**2)
tupet_aa_alt34n1[2]=sqrt(sum(amch3[18,ialt34_ach3n1]**2)/size(amch3[18,ialt34_ach3n1])-mean(amch3[18,ialt34_ach3n1])**2)
tupet_aa_alt34n1[3]=sqrt(sum(amch4[18,ialt34_ach4n1]**2)/size(amch4[18,ialt34_ach4n1])-mean(amch4[18,ialt34_ach4n1])**2)
tupet_aa_alt34n1[4]=sqrt(sum(amch5[18,ialt34_ach5n1]**2)/size(amch5[18,ialt34_ach5n1])-mean(amch5[18,ialt34_ach5n1])**2)
tupet_aa_alt34n1[5]=sqrt(sum(amch6[18,ialt34_ach6n1]**2)/size(amch6[18,ialt34_ach6n1])-mean(amch6[18,ialt34_ach6n1])**2)
tupet_aa_alt34n1[6]=sqrt(sum(amch7[18,ialt34_ach7n1]**2)/size(amch7[18,ialt34_ach7n1])-mean(amch7[18,ialt34_ach7n1])**2)
tupet_aa_alt34n1[7]=sqrt(sum(amch15[18,ialt34_ach15n1]**2)/size(amch15[18,ialt34_ach15n1])-mean(amch15[18,ialt34_ach15n1])**2)

#ssmis
ialt34_ssch1=(ssch1[11,:]>=3000)&(ssch1[11,:]<=4000)
ialt34_ssch2=(ssch2[11,:]>=3000)&(ssch2[11,:]<=4000)
ialt34_ssch3=(ssch3[11,:]>=3000)&(ssch3[11,:]<=4000)
ialt34_ssch4=(ssch4[11,:]>=3000)&(ssch4[11,:]<=4000)
ialt34_ssch5=(ssch5[11,:]>=3000)&(ssch5[11,:]<=4000)
ialt34_ssch8=(ssch8[11,:]>=3000)&(ssch8[11,:]<=4000)
ialt34_ssch9=(ssch9[11,:]>=3000)&(ssch9[11,:]<=4000)
ialt34_ssch10=(ssch10[11,:]>=3000)&(ssch10[11,:]<=4000)
ialt34_ssch11=(ssch11[11,:]>=3000)&(ssch11[11,:]<=4000)
ialt34_ssch12=(ssch12[11,:]>=3000)&(ssch12[11,:]<=4000)
ialt34_ssch13=(ssch13[11,:]>=3000)&(ssch13[11,:]<=4000)
ialt34_ssch14=(ssch14[11,:]>=3000)&(ssch14[11,:]<=4000)
ialt34_ssch15=(ssch15[11,:]>=3000)&(ssch15[11,:]<=4000)
ialt34_ssch16=(ssch16[11,:]>=3000)&(ssch16[11,:]<=4000)
ialt34_ssch17=(ssch17[11,:]>=3000)&(ssch17[11,:]<=4000)
ialt34_ssch18=(ssch18[11,:]>=3000)&(ssch18[11,:]<=4000)

tup_ss1_alt34=mean(ssch1[16,ialt34_ssch1])
tup_ss2_alt34=mean(ssch2[16,ialt34_ssch2])
tup_ss3_alt34=mean(ssch3[16,ialt34_ssch3])
tup_ss4_alt34=mean(ssch4[16,ialt34_ssch4])
tup_ss5_alt34=mean(ssch5[16,ialt34_ssch5])
tup_ss8_alt34=mean(ssch8[16,ialt34_ssch8])
tup_ss9_alt34=mean(ssch9[16,ialt34_ssch9])
tup_ss10_alt34=mean(ssch10[16,ialt34_ssch10])
tup_ss11_alt34=mean(ssch11[16,ialt34_ssch11])
tup_ss12_alt34=mean(ssch12[16,ialt34_ssch12])
tup_ss13_alt34=mean(ssch13[16,ialt34_ssch13])
tup_ss14_alt34=mean(ssch14[16,ialt34_ssch14])
tup_ss15_alt34=mean(ssch15[16,ialt34_ssch15])
tup_ss16_alt34=mean(ssch16[16,ialt34_ssch16])
tup_ss17_alt34=mean(ssch17[16,ialt34_ssch17])
tup_ss18_alt34=mean(ssch18[16,ialt34_ssch18])

#ecart type
tupet_ss_alt34=np.zeros([16],float)
tupet_ss_alt34[0]=sqrt(sum(ssch1[16,ialt34_ssch1]**2)/size(ssch1[16,ialt34_ssch1])-mean(ssch1[16,ialt34_ssch1])**2)
tupet_ss_alt34[1]=sqrt(sum(ssch2[16,ialt34_ssch2]**2)/size(ssch2[16,ialt34_ssch2])-mean(ssch2[16,ialt34_ssch2])**2)
tupet_ss_alt34[2]=sqrt(sum(ssch3[16,ialt34_ssch3]**2)/size(ssch3[16,ialt34_ssch3])-mean(ssch3[16,ialt34_ssch3])**2)
tupet_ss_alt34[3]=sqrt(sum(ssch4[16,ialt34_ssch4]**2)/size(ssch4[16,ialt34_ssch4])-mean(ssch4[16,ialt34_ssch4])**2)
tupet_ss_alt34[4]=sqrt(sum(ssch5[16,ialt34_ssch5]**2)/size(ssch5[16,ialt34_ssch5])-mean(ssch5[16,ialt34_ssch5])**2)
tupet_ss_alt34[5]=sqrt(sum(ssch8[16,ialt34_ssch8]**2)/size(ssch8[16,ialt34_ssch8])-mean(ssch8[16,ialt34_ssch8])**2)
tupet_ss_alt34[6]=sqrt(sum(ssch9[16,ialt34_ssch9]**2)/size(ssch9[16,ialt34_ssch9])-mean(ssch9[16,ialt34_ssch9])**2)
tupet_ss_alt34[7]=sqrt(sum(ssch10[16,ialt34_ssch10]**2)/size(ssch10[16,ialt34_ssch10])-mean(ssch10[16,ialt34_ssch10])**2)
tupet_ss_alt34[8]=sqrt(sum(ssch11[16,ialt34_ssch11]**2)/size(ssch11[16,ialt34_ssch11])-mean(ssch11[16,ialt34_ssch11])**2)
tupet_ss_alt34[9]=sqrt(sum(ssch12[16,ialt34_ssch12]**2)/size(ssch12[16,ialt34_ssch12])-mean(ssch12[16,ialt34_ssch12])**2)
tupet_ss_alt34[10]=sqrt(sum(ssch13[16,ialt34_ssch13]**2)/size(ssch13[16,ialt34_ssch13])-mean(ssch13[16,ialt34_ssch13])**2)
tupet_ss_alt34[11]=sqrt(sum(ssch14[16,ialt34_ssch14]**2)/size(ssch14[16,ialt34_ssch14])-mean(ssch14[16,ialt34_ssch14])**2)
tupet_ss_alt34[12]=sqrt(sum(ssch15[16,ialt34_ssch15]**2)/size(ssch15[16,ialt34_ssch15])-mean(ssch15[16,ialt34_ssch15])**2)
tupet_ss_alt34[13]=sqrt(sum(ssch16[16,ialt34_ssch16]**2)/size(ssch16[16,ialt34_ssch16])-mean(ssch16[16,ialt34_ssch16])**2)
tupet_ss_alt34[14]=sqrt(sum(ssch17[16,ialt34_ssch17]**2)/size(ssch17[16,ialt34_ssch17])-mean(ssch17[16,ialt34_ssch17])**2)
tupet_ss_alt34[15]=sqrt(sum(ssch18[16,ialt34_ssch18]**2)/size(ssch18[16,ialt34_ssch18])-mean(ssch18[16,ialt34_ssch18])**2)

tupnd_ss_alt34=np.zeros([16],float)
tupet_ss_alt34[0]=sum(ialt34_ssch1)
tupet_ss_alt34[1]=sum(ialt34_ssch2)
tupet_ss_alt34[2]=sum(ialt34_ssch3)
tupet_ss_alt34[3]=sum(ialt34_ssch4)
tupet_ss_alt34[4]=sum(ialt34_ssch5)
tupet_ss_alt34[5]=sum(ialt34_ssch8)
tupet_ss_alt34[6]=sum(ialt34_ssch9)
tupet_ss_alt34[7]=sum(ialt34_ssch10)
tupet_ss_alt34[8]=sum(ialt34_ssch11)
tupet_ss_alt34[9]=sum(ialt34_ssch12)
tupet_ss_alt34[10]=sum(ialt34_ssch13)
tupet_ss_alt34[11]=sum(ialt34_ssch14)
tupet_ss_alt34[12]=sum(ialt34_ssch15)
tupet_ss_alt34[13]=sum(ialt34_ssch16)
tupet_ss_alt34[14]=sum(ialt34_ssch17)
tupet_ss_alt34[15]=sum(ialt34_ssch18)




for i in range(0,16):
    ti=str(ssmisch[i])
    t='tupet_ss_alt34['+str(i)+']=sum(ialt34_ssch'+ti+')'
    print t



ssmisch=np.array([1,2,3,4,5,8,9,10,11,12,13,14,15,16,17,18])
for i in range(0,16):
    ti=str(ssmisch[i])
    t='tup_ss_alt34['+str(i)+']=sqrt(sum(ssch'+ti+'[16,ialt34_ssch'+ti+'])/size(ssch'+ti+'[16,ialt34_ssch'+ti+'])+mean(ssch'+ti+'+[16,ialt34_ssch'+ti+']))'
    print t