擅长:python、mysql、java
<p>最好的方法可能是使用np.Array.sum(axis=1)。假设g(Z[i,j])调用是有效的。实际上,您甚至不需要任何循环:</p>
<pre><code>import numpy as np
Z = np.random.randint(0,255,(10,15))
g=np.random.randint(0,10,(256))
lEks = np.random.rand((15))
def genHDR(Z,g,lEks):
numerator = (Z*(g[Z]-lEks.view().reshape((1,)+ lEks.shape))).sum(axis=1)
denominator = Z.sum(axis=1)
HDRimage = numerator/denominator
return HDRimage
</code></pre>