2

这是我正在处理的代码的一部分。看最后一行。

def xyz2srgb(xyz):
    """Convert from XYZ colorspace to sRGB"""
    xyz = xyz.T
    srgb = N.zeros_like(xyz)
    srgb[0] = 3.241*xyz[0] - 1.5374*xyz[1] - 0.4986*xyz[2]
    srgb[1] = -0.9692*xyz[0] + 1.876*xyz[1] + 0.0416*xyz[2]
    srgb[2] = 0.0556*xyz[0] - 0.204*xyz[1] + 1.057*xyz[2]

    srgb /= float(srgb.max())

    a = 0.055
    srgb = N.where(srgb <= 0.0031308,
                12.92*srgb,
                (1+a)*srgb**(1/2.4)-a)

它抛出一个错误:

 ./example.py:29: RuntimeWarning: invalid value encountered in power
  (1+a)*srgb**(1/2.4)-a)

这是一些范围限制吗?我在 SO 上找到了一个类似的线程:Numpy 错误:在电源中遇到无效值

有什么线索吗?谢谢。

4

2 回答 2

4

您应该打印出srgb. 可能不是范围误差,是这里的负系数

srgb[0] = 3.241*xyz[0] - 1.5374*xyz[1] - 0.4986*xyz[2]
srgb[1] = -0.9692*xyz[0] + 1.876*xyz[1] + 0.0416*xyz[2]
srgb[2] = 0.0556*xyz[0] - 0.204*xyz[1] + 1.057*xyz[2]

使元素的元素完全有可能srgb是负数,除非你已经在处理复数,否则当你掌权时你会发出警告:

>>> import numpy
>>> numpy.array([-2.0])
array([-2.])
>>> numpy.array([-2.0])**(1/2.4)
__main__:1: RuntimeWarning: invalid value encountered in power
array([ nan])
>>> numpy.array([-2.0+0j])**(1/2.4)
array([ 0.34548198+1.28935629j])

参数where将在传递之前进行评估,因此将在srgb <= 0.0031308条件分支划分值之前获取权力,并触发警告。

于 2012-05-14T11:29:25.037 回答
0

对于这些错误,您可以使用指数中的浮点数。例如 :

>>> test = np.linspace(0,10,1000);)

>>> expos = 2*np.power(test,**6.**)+2*np.power(test,**4.**)

此致

于 2016-11-18T13:04:37.700 回答