擅长:python、mysql、java
<p>好吧,多亏了你的帮助和你介绍给我的帖子,我终于明白了。在</p>
<p>计算我的全部数据(300000+)大约需要一秒钟!在</p>
<p>我使用了以下代码:</p>
<pre><code>def runningmean(x,N):
cumsum =np.cumsum(np.insert(x,0,0))
return (cumsum[N:] -cumsum[:-N])/N
photosignal = np.array(seg.analogsignal[0], ndmin =1)
photosignalaverage = runningmean(photosignal, 3000)
holder = np.zeros(2999)
photosignalaverage = np.append(holder,photosignalaverage)
detalfsignal = (photosignal-photosignalaverage)/abs(photosignalaverage)
</code></pre>
<p>Photosignal将我的原始信号存储在numpy数组中。
Photosignalaverage使用cumsum计算photosignal中每个数据点的运行平均值。然后,我将前2999个值作为0添加到maintian,其大小与我的photosignal相同。在</p>
<p>然后我用基本的numpy计算得到delta F/F信号。在</p>
<p>再次感谢你的反馈,真的很有帮助!在</p>