python| 求在2°与10°视场下颜色的三刺激值_python 相对光谱功率分布-CSDN博客

本文链接：https://blog.csdn.net/Wangwenshuaicsdn/article/details/106268677
本程序提前导入了不同波长在不同光源不同视场中的反射率，相对光谱功率分布以及光谱三刺激值。学识浅薄，希望大佬能不吝赐教。
#        构建光谱反射率元组
tupleR =(0.0231,0.0206,0.018,0.0177,0.0181,0.0179,
0.0179,0.0182,0.0188,0.0187,0.0208,0.0306,
0.0519,0.067,0.0623,0.0549,0.0515,0.0468,
0.0396,0.0353,0.0466,0.1389,0.3418,
0.5255,0.6282,0.6757,0.6998,0.7153,0.731,
0.7422,0.7492,0.7517,0.752,0.7527,0.7578,0.7698,0,0,0,0,0)



#         分色光构建相对光谱功率分布
tupleD50 = (24.49,29.87,49.31,56.51,60.03,57.82,
74.82,87.25,90.61,91.37,95.11,91.96,
95.72,96.61,97.13,102.10,100.75,102.32,100.0,
97.74,98.92,93.50,97.69,99.27,99.04,95.72,
98.86,95.67,98.19,103.00,99.13,87.38,91.60,
92.89,76.85,86.51,92.58,78.23,57.69,82.92,78.27)

tupleD65 = (49.98,54.65,82.75,91.49,93.43,86.68,104.87,
100.00,117.01,117.81,114.86,115.92,108.81,
109.35,107.80,104.79,107.69,104.41,104.05,
96.33,95.79,88.69,90.01,89.60,87.70,83.29,
83.70,80.03,80.21,82.28,78.28,69.72,71.61,
74.35,61.60,69.89,75.09,63.59,46.42,66.81,63.38,)

tupleC = (33.00,47.40,63.30,80.60,98.10,112.40,121.50,
124.00,123.10,123.80,123.90,120.70,112.10,
102.30,96.90,98.00,102.10,105.20,105.30,102.30,
97.80,93.20,89.70,88.40,88.10,88.00,87.80,88.20,
87.90,86.30,84.00,80.20,76.30,72.40,68.30,64.40,
61.50,59.20,58.10,58.20,59.10,)

tupleA = (9.80,12.09,14.71,17.68,21.00,24.67,
28.70,33.09,37.81,42.87,48.24,53.91,
59.86,66.06,72.50,79.13,85.95,92.91,
100,107.18,114.44,121.73,129.04,136.35,
143.62,150.84,157.98,165.03,171.96,
178.77,185.43,191.93,198.26,204.41,
210.37,216.12,221.67,227.00,232.12,237.01,241.68)
#        构建波长列表
list1 = (380,390,400,410,420,430,440,
450,460,470,480,490,500,510,
520,530,540,550,560,570,580,
590,600,610,620,630,640,650,
660,670,680,690,700,710,720,
730,740,750,760,770,780)

#         光谱三刺激值
tupleX2=(0.001368,0.004243,0.01431,0.04351,0.13438,
0.2829,0.34828,0.3362,0.2908,0.19536,0.09564,
0.03201,0.0049,0.0093,0.06327,0.1655,0.2904,
0.43345,0.5945,0.7621,0.9163,1.0263,1.0622,1.0026,
0.85445,0.6424,0.4479,0.2835,0.1649,0.0874,0.04677,
0.0227,0.011359,0.00579,0.002899,0.00144,0.00069,
0.000332,0.000166,0.0000831,0.0000415)


tupleY2 = (0.000039,0.00012,0.000396,0.00121,0.004,
0.0116,0.023,0.038,0.06,0.09098,0.13902,
0.20802,0.323,0.503,0.71,0.862,0.954,0.99495,
0.995,0.952,0.87,0.757,0.631,0.503,0.381,
0.265,0.175,0.107,0.061,0.032,0.017,0.00821,
0.004102,0.002091,0.001047,0.00052,0.000249,
0.00012,0.00006,0.00003,0.000015)


tupleZ2 = (0.00645,0.02005,0.06785,0.2074,0.6456,
1.3856,1.74706,1.77211,1.6692,1.28764,
0.81295,0.46518,0.272,0.1582,
0.07825,0.04216,0.0203,0.00875,0.0039,
0.0021,0.00165,0.0011,0.0008,0.00034,
0.00019,0.00005,0.00002,0,0,0,0,0,0,0,
0,0,0,0,0,0,0)


tupleX10= (0.000160,0.002362,0.019110,0.084736,0.204492,
0.314679,0.383734,0.370702,0.302273,0.195618,
0.080507,0.016172,0.003816,0.037465,0.117749,
0.236491,0.376772,0.529826,0.705224,0.878655,
1.014160,1.118520,1.123990,1.030480,0.856297,
0.647467,0.431567,0.268329,0.152568,0.081261,
0.040851,0.019941,0.009577,0.004553,0.002175,
0.001045,0.000508,0.000251,0.000126,0.000065,0.000033,)

tupleY10 = (0.000017,0.000253,0.002004,0.008756,0.021391,
0.038676,0.062077,0.089456,0.128201,0.185190,
0.253589,0.339133,0.460777,0.606741,0.761757,
0.875211,0.961988,0.991761,0.997340,0.955552,
0.868934,0.777405,0.658341,0.527963,0.398057,
0.283493,0.179828,0.107633,0.060281,0.031800,
0.015905,0.007749,0.003718,0.001768,0.000846,
0.000407,0.000199,0.000098,0.000050,0.000025,0.000013,)

tupleZ10 = (0.000705,0.010482,0.086011,0.389366,
0.972542,1.553480,1.967280,1.994800,
1.745370,1.317560,0.772125,0.415254,
0.218502,0.112044,0.060709,0.030451,
0.013676,0.003988,0.000000,0.000000,
0.000000,0.000000,0.000000,0.000000,
0.000000,0.000000,0.000000,0.000000,
0.000000,0.000000,0.000000,0.000000,
0.000000,0.000000,0.000000,0.000000,
0.000000,0.000000,0.000000,0.000000,0.000000,)


#         初始化
k = 0.095
dict1 = dict(zip(list1,tupleR))

dictD50 = dict(zip(list1,tupleD50))
dictD65 = dict(zip(list1,tupleD65))
dictC = dict(zip(list1,tupleC))
dictA = dict(zip(list1,tupleA))

dict2X = dict(zip(list1,tupleX2))
dict2Y = dict(zip(list1,tupleY2))
dict2Z = dict(zip(list1,tupleZ2))

dict10X = dict(zip(list1,tupleX10))
dict10Y = dict(zip(list1,tupleY10))
dict10Z = dict(zip(list1,tupleZ10))




def d502():
    sum1 = 0
    sum2 = 0
    sum3 = 0
    for i in list1:
        sum1 = sum1 + k*dict1[i]*dictD50[i]*dict2X[i]
        sum2 = sum2 + k*dict1[i]*dictD50[i]*dict2Y[i]
        sum3 = sum3 + k*dict1[i]*dictD50[i]*dict2Z[i]
    print("X =%.1f"%sum1)
    print("Y =%.1f"%sum2)
    print("Z =%.1f"%sum3)

def  d5010():
    sum1 = 0
    sum2 = 0
    sum3 = 0
    for i in list1:
        sum1 = sum1 + k*dict1[i]*dictD50[i]*dict10X[i]
        sum2 = sum2 + k*dict1[i]*dictD50[i]*dict10Y[i]
        sum3 = sum3 + k*dict1[i]*dictD50[i]*dict10Z[i]
    print("X =%.1f"%sum1)
    print("Y =%.1f"%sum2)
    print("Z =%.1f"%sum3)


def d652():
    sum1 = 0
    sum2 = 0
    sum3 = 0
    for i in list1:
        sum1 = sum1 + k*dict1[i]*dictD65[i]*dict2X[i]
        sum2 = sum2 + k*dict1[i]*dictD65[i]*dict2Y[i]
        sum3 = sum3 + k*dict1[i]*dictD65[i]*dict2Z[i]
    print("X =%.1f"%sum1)
    print("Y =%.1f"%sum2)
    print("Z =%.1f"%sum3)

def  d6510():
    sum1 = 0
    sum2 = 0
    sum3 = 0
    for i in list1:
        sum1 = sum1 + k*dict1[i]*dictD65[i]*dict10X[i]
        sum2 = sum2 + k*dict1[i]*dictD65[i]*dict10Y[i]
        sum3 = sum3 + k*dict1[i]*dictD65[i]*dict10Z[i]
    print("X =%.1f"%sum1)
    print("Y =%.1f"%sum2)
    print("Z =%.1f"%sum3)

def  c2():
    sum1 = 0
    sum2 = 0
    sum3 = 0
    for i in list1:
        sum1 = sum1 + k*dict1[i]*dictC[i]*dict2X[i]
        sum2 = sum2 + k*dict1[i]*dictC[i]*dict2Y[i]
        sum3 = sum3 + k*dict1[i]*dictC[i]*dict2Z[i]
    print("X =%.1f"%sum1)
    print("Y =%.1f"%sum2)
    print("Z =%.1f"%sum3)

def  c10():
    sum1 = 0
    sum2 = 0
    sum3 = 0
    for i in list1:
        sum1 = sum1 + k*dict1[i]*dictC[i]*dict10X[i]
        sum2 = sum2 + k*dict1[i]*dictC[i]*dict10Y[i]
        sum3 = sum3 + k*dict1[i]*dictC[i]*dict10Z[i]
    print("X =%.1f"%sum1)
    print("Y =%.1f"%sum2)
    print("Z =%.1f"%sum3)

def  a2():
    sum1 = 0
    sum2 = 0
    sum3 = 0
    for i in list1:
        sum1 = sum1 + k*dict1[i]*dictA[i]*dict2X[i]
        sum2 = sum2 + k*dict1[i]*dictA[i]*dict2Y[i]
        sum3 = sum3 + k*dict1[i]*dictA[i]*dict2Z[i]
    print("X =%.1f"%sum1)
    print("Y =%.1f"%sum2)
    print("Z =%.1f"%sum3)

def  a10():
    sum1 = 0
    sum2 = 0
    sum3 = 0
    for i in list1:
        sum1 = sum1 + k*dict1[i]*dictA[i]*dict10X[i]
        sum2 = sum2 + k*dict1[i]*dictA[i]*dict10Y[i]
        sum3 = sum3 + k*dict1[i]*dictA[i]*dict10Z[i]
    print("X =%.1f"%sum1)
    print("Y =%.1f"%sum2)
    print("Z =%.1f"%sum3)


   

#    程序运行
sight = input('光源：')
print('视场：',end = '')
hor = input()
if (sight =='D50'):
    if hor == '2':
        d502()
    else:
        d5010()


if sight =='D65':
        if hor == '2':
            d652()
        else:
            d6510()
          


if sight =='C':
        if hor == '2':
            c2()
        else:
            c10()



if sight =='A':
        if hor == '2':
            a2()
        else:
            a10()