无论我如何更改数据,通过以下代码绘制的图形只是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) 我试图将几个参数发布到此[url] [1]并按"submit"下载生成的csv文件.
我认为至少需要5个步骤.
我有一个三列(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) 帮助链接是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)之间的关系是什么?
我需要arccos()在bash脚本中计算.
gawk可以在linux中计算cos(theta)和sin(theta)
如何计算arccos()?
对于一个图像的每个像素,我都有其x(int),y(int)和像素值(浮点数)。
现在我需要为图像添加噪点。numpy.random.poisson是否合适?
我很担心,因为它不像新像素值=原始值+噪声,而是
新像素值= numpy.random.poisson(原始值,1)。新值都是整数。
我的问题是标题。
我的目的是获得恒星的光度测量误差,但是我只有一个图像,所以我通过添加泊松噪声进行了模拟。请检查以下ccd图像中的数字。来源是红色功能。

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

我无法在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" .
他们的版本在EPD中很老.
我试过pip,但似乎它们被封装在EPD中.
有没有办法升级它们?
python ×7
numpy ×4
scipy ×3
selenium ×2
bash ×1
epd-python ×1
fft ×1
histogram ×1
image ×1
linux ×1
matplotlib ×1
mayavi ×1
noise ×1
statistics ×1
web-scraping ×1