小编que*_*ang的帖子

功率谱由numpy.fft.fft

无论我如何更改数据,通过以下代码绘制的图形只是ZERO周围的峰值.我的数据只是一列,记录了某种信号的每个时间点.是time_step一个值,我应该根据两个相邻点在我的数据的时间间隔定义?

data=np.loadtxt("timesequence",delimiter=",",usecols=(0,),unpack=True)

ps = np.abs(np.fft.fft(data))**2
time_step = 1

freqs = np.fft.fftfreq(data.size, time_step)
idx   = np.argsort(freqs)

pl.plot(freqs[idx], ps[idx])
pl.show()
Run Code Online (Sandbox Code Playgroud)

python signal-processing numpy fft

8
推荐指数
1
解决办法
1万
查看次数

如何在python中使用selenium定位四个元素

我试图将几个参数发布到此[url] [1]并按"submit"下载生成的csv文件.

我认为至少需要5个步骤.

python selenium web-scraping

8
推荐指数
1
解决办法
340
查看次数

如何使用scipy.interpolate中的griddata

我有一个三列(x像素,y像素,z值)数据,有一百万行.数据来自图像,并且存在重复的z值.现在我需要制作表面图.这个图像是一个很好的例子.但现在输出图像为空.有人可以查一下代码吗?

import numpy as np
from enthought.mayavi import mlab
from scipy.interpolate import griddata
x,y,z = np.loadtxt('test.csv',delimiter=',',usecols=(0,1,2),unpack=True) 
xi,yi = np.mgrid[0:3000:3000j, 0:3000:3000j]
zi = griddata((x, y), z, (xi, yi),method='linear')
mlab.surf(xi,yi,zi)
mlab.show()
Run Code Online (Sandbox Code Playgroud)

python interpolation numpy scipy mayavi

6
推荐指数
1
解决办法
1万
查看次数

如何使用scipy.stats.kstest /关于Kolmogorov-Smirnov测试的基本问题

帮助链接是http://docs.scipy.org/doc/scipy-0.7.x/reference/generated/scipy.stats.kstest.html 我现在可以计算ks-test值,但我不明白.代码如下.

from scipy import stats
import numpy as np
sample =np.loadtxt('mydata',delimiter=",",usecols=(2,),unpack=True)
print stats.kstest(sample, 'poisson', args=(1,))
Run Code Online (Sandbox Code Playgroud)

Q1
如果参考分布是常数,那么什么词可以代替上面的"泊松"?
Q2
是什么意思args=(1,)
Q3
如果有人对ks-test感兴趣,这里是wiki链接.
http://en.wikipedia.org/wiki/Kolmogorov%E2%80%93Smirnov_test
我们可以编写自己的python代码来练习吗?我们可以轻松获得max(D),但是如何在链接中获得Pr(k <= x)max(D)Pr(k <= x)之间的关系是什么?

python statistics scipy

4
推荐指数
1
解决办法
6610
查看次数

如何在bash中计算arccos()?

我需要arccos()在bash脚本中计算.
gawk可以在linux中计算cos(theta)sin(theta)
如何计算arccos()

linux bash

4
推荐指数
2
解决办法
3566
查看次数

numpy.random.poisson可以用来为图像添加泊松噪声吗?

对于一个图像的每个像素,我都有其x(int),y(int)和像素值(浮点数)。

现在我需要为图像添加噪点。numpy.random.poisson是否合适?

我很担心,因为它不像新像素值=原始值+噪声,而是

新像素值= numpy.random.poisson(原始值,1)。新值都是整数。

我的问题是标题。

我的目的是获得恒星的光度测量误差,但是我只有一个图像,所以我通过添加泊松噪声进行了模拟。请检查以下ccd图像中的数字。来源是红色功能。

在此处输入图片说明

numpy image image-processing noise

3
推荐指数
1
解决办法
2107
查看次数

如何通过matplotlib缩放直方图

你可以看到下面有直方图.
它像是 pl.hist(data1,bins=20,color='green',histtype="step",cumulative=-1)
如何缩放直方图?
例如,让直方图的高度为现在的三分之一.

此外,它是一种删除左侧垂直线的方法吗?

在此输入图像描述

python matplotlib histogram

2
推荐指数
1
解决办法
5156
查看次数

如何正确启动Selenium中的Chrome

我无法在Selenium中启动Chrome.

driver=webdriver.Chrome()
Run Code Online (Sandbox Code Playgroud)

Traceback(最近一次调用最后一次):

  File "<stdin>", line 1, in <module>

  File "/usr/anaconda/lib/python2.7/site-packages/selenium/webdriver/chrome/webdriver.py", line 67, in __init__
    desired_capabilities=desired_capabilities)

  File "/usr/anaconda/lib/python2.7/site-packages/selenium/webdriver/remote/webdriver.py", line 87, in __init__
    self.start_session(desired_capabilities, browser_profile)

  File "/usr/anaconda/lib/python2.7/site-packages/selenium/webdriver/remote/webdriver.py", line 141, in start_session
    'desiredCapabilities': desired_capabilities,
  File "/usr/anaconda/lib/python2.7/site-packages/selenium/webdriver/remote/webdriver.py", line 201, in execute
    self.error_handler.check_response(response)

  File "/usr/anaconda/lib/python2.7/site-packages/selenium/webdriver/remote/errorhandler.py", line 181, in check_response
    raise exception_class(message, screen, stacktrace)
selenium.common.exceptions.WebDriverException: Message: unknown error: Chrome failed to start: crashed
  (Driver info: chromedriver=2.33.506092 (733a02544d189eeb751fe0d7ddca79a0ee28cce4),platform=Linux 3.19.8-100.fc20.x86_64 x86_64)
Run Code Online (Sandbox Code Playgroud)

我从网址使用"chromedriver_linux64.zip 2017-10-03 21:09:52 3.90MB" .

python selenium google-chrome

2
推荐指数
1
解决办法
1683
查看次数

如何在EPD升级scipy或numpy?

他们的版本在EPD中很老.
我试过pip,但似乎它们被封装在EPD中.
有没有办法升级它们?

python numpy scipy epd-python

1
推荐指数
1
解决办法
183
查看次数