我想删除numpy.array中的选定列.这就是我做的:
n [397]: a = array([[ NaN, 2., 3., NaN],
.....: [ 1., 2., 3., 9]])
In [398]: print a
[[ NaN 2. 3. NaN]
[ 1. 2. 3. 9.]]
In [399]: z = any(isnan(a), axis=0)
In [400]: print z
[ True False False True]
In [401]: delete(a, z, axis = 1)
Out[401]:
array([[ 3., NaN],
[ 3., 9.]])
在此示例中,我的目标是删除包含NaN的所有列.我希望最后一个命令导致:
array([[2., 3.],
[2., 3.]])
我怎样才能做到这一点？
假设我在matplotlib中有两个数字,每个数字有一个数字:
import matplotlib.pyplot as plt
f1 = plt.figure()
plt.plot(range(0,10))
f2 = plt.figure()
plt.plot(range(10,20))
然后我一次性展示
plt.show()
有没有办法单独显示它们,即只显示f1？
或者更好:我如何单独管理数字,如下面的"如意"代码(不起作用):
f1 = plt.figure()
f1.plot(range(0,10))
f1.show()
我在Numpy/Scipy/Matplotlib中编写科学软件.在家用计算机上开发应用程序后,我现在对编写简单的Web应用程序感兴趣.示例:用户上传图像或音频文件,我的程序使用Numpy/Scipy处理它,输出使用Matplotlib显示在浏览器上,或者用户可以下载已处理的文件.
我已经支付了安装了Python 2.4.3的托管,但没有Numpy/Scipy.我也没有通过命令行访问shell.只需拖放FTP即可.非常有限,但我可以使用简单的Python/CGI脚本.
令人惊讶的是,网络搜索显示很少有适合网络托管的选项,内置了这些功能.(如果我错了,请指导我.)我正在学习Google App Engine,但我仍然没有完全了解它工具和限制.什么是网络没告诉我的是,别人也有类似的担忧.
希望找到解决方案,我想我会向这个令人敬畏的SO社区提出这些简单的问题:
是否有一种简单的方法可以将numpy(或任何第三方软件包/库)安装到我已经托管的空间中？我知道托管空间上的Python路径,我知道家用计算机上的相关Python/Numpy目录.我可以简单地复制文件并让它工作吗？本地和远程系统都运行Ubuntu.
存在Numpy/Matplotlib的托管站点(免费或付费),或者如果没有安装,是否可以安装它？是否有任何文档化的站点可以使用工作应用程序引用,无论多么简单？
Google App Engine能以任何方式帮助我吗？还是完全是为了别的什么？你或其他人用它来编写Python/Numpy中的科学应用程序吗？如果是这样,你能参考一下吗？
谢谢您的帮助.
编辑:在下面的有用答案之后,我在Slicehost购买了20美元的计划,到目前为止我喜欢它!(我首先尝试使用Amazon EC2.我必须是愚蠢的,但我无法让它工作.)使用Apache设置Ubuntu服务器需要花费几个小时(我是Apache新手).它允许我完成我想要的Python以及更多.我现在也拥有自己的版本控制远程存储库.再次感谢!
编辑2:近两年后,我尝试了Linode和EC2(再次).Linode很棒.这次EC2似乎更容易 - 也许它只是增加了经验,或者可能是亚马逊对AWS管理控制台所做的改进.对于那些对Numpy/Scipy/Matplotlib/Audiolab感兴趣的人,每当我启动EC2实例时,这是我的Ubuntu备忘单:
ec2:~$ sudo aptitude install build-essential python-scipy ipython
python-matplotlib python-dev python-setuptools libsndfile-dev
libasound2-dev mysql-server python-mysqldb
Upload scikits.audiolab-0.11.0
ec2:~/scikits.audiolab-0.11.0$ sudo python setup.py install
ec2:~$ sudo rm -rf scikits.audiolab-0.11.0
ec2:~$ nano .ipython/ipy_user_conf.py
ip.ex('import matplotlib; matplotlib.use("Agg"); import scipy, pylab,
scipy.signal as sig, scipy.linalg as lin, scipy.sparse as spar,
os, sys, MySQLdb, boto; from scikits import audiolab')
import ipy_greedycompleter
import ipy_autoreload
我发现自己经常编写这样的代码:
k = 0
for i in mylist:
# y[k] = some function of i
k += 1
相反,我能做到
for k in range(K):
# y[k] = some function of mylist[k]
但这似乎不是"pythonic".(你知道......索引.Ick!)是否有一些语法允许我使用循环,列表推导或生成器同时提取索引(k)和元素(i)？任务是在科学计算中,因此循环体中有很多东西,我认为列表理解可能不够强大.
我也欢迎有关相关概念的提示,即使我可能也没有.谢谢.
我想用Pylab绘制一个混淆矩阵.沿水平轴的类标签很长,所以我想绘制它们垂直旋转.但是,我还想将它们绘制在轴的顶部,而不是下面.
此命令可以在底部绘制垂直标签:
pylab.imshow(confusion_matrix)
pylab.xticks(..., rotation='vertical')
并且此命令可以在顶部绘制水平标签而不旋转:
pylab.matshow(confusion_matrix)
但我找不到两者兼而有之.以下命令不起作用.
pylab.matshow(confusion_matrix)
pylab.xticks(..., rotation='vertical')
你能建议一种方法用垂直旋转在轴顶部用xticks绘制混淆矩阵吗？谢谢.
编辑
马克,谢谢你的帮助.通过更仔细地检查刻度属性,它让我走上正轨.与你的答案和我想要的答案的唯一区别是将这个想法应用于AxesImage,而不是绘图.经过调查,这里是答案:
im = pylab.matshow(confusion_matrix)
for label in im.axes.xaxis.get_ticklabels():
label.set_rotation(90)
im.figure.show()
对于那些阅读...不要忘记show()!我忘了我需要刷新数字.见下面的输出.
当我尝试在一些高维输入上应用一些ML算法(分类,更具体地说,特别是SVM)时,我应该如何处理,我得到的结果不太令人满意？
可以显示1,2或3维数据以及算法的结果,这样您就可以了解正在发生的事情,并了解如何解决问题.一旦数据超过3个维度,除了直观地使用参数,我不确定如何攻击它？
我有一个矩阵P与形状MxN和三维张量T与形状KxNxR.我想P与每个NxR矩阵相乘T,产生一个KxMxR3d张量.
P.dot(T).transpose(1,0,2)给出了期望的结果.有没有更好的解决方案(即摆脱transpose)这个问题？这必须是一个非常常见的操作,所以我假设,其他人已经找到了不同的方法,例如使用tensordot(我试过但未能获得所需的结果).意见/观点将受到高度赞赏!
什么是最快(或最"Pythonic")的转换方式
x = [False, False, True, True]
进12？(如果有这样的方式.)
如果x是一个numpy.arraybools怎么办？对此有特殊的命令吗？
我有一个大的m-by-n布尔数组,其中每个n元素行代表一个高维特征向量的低维散列.(在上面的例子中,n = 4.)我想知道答案,以便尽可能地压缩我的数据.谢谢.
编辑:谢谢你的回复!使用以下测试代码,
t = 0
for iter in range(500):
B = scipy.signbit(scipy.randn(1000,20))
for b in B:
t0 = time.clock()
# test code here
t1 = time.clock()
t += (t1-t0)
print t
...这是我的Thinkpad笔记本电脑上的运行时间:
当然,我欢迎任何可能确认或反驳我的数据的独立测试!
编辑:在我的答案下面,int(j)改为简单j仍然有效,但运行速度慢六倍!然后,如果使用bool进行铸造,其他答案可能会更快int.但我懒得再测试一切.
我希望将图片准确地放在文本中的特定位置.我使用下面的命令
\begin{enumerate}
\item T.D. Lee(1957 Physics Nobel Laureate)
\begin{figure}
\begin{center}
\includegraphics[scale=0.5]{TdLee.eps}
\end{center}
\end{figure}
我有很多像这样的数字.会发生什么是我看到不同顺序的数字和项目.你可以让我知道我应该使用哪种乳胶命令,以便将图片放在每个项目的下方或右侧？
谢谢
库尔特
我想要实现的目标如下:我需要声音文件(.wav)的频率值进行分析.我知道很多程序会给出值的可视图(谱图),但我需要原始数据.我知道这可以通过FFT完成,并且应该在python中相当容易编写脚本但不确定如何完全执行它.因此,假设文件中的信号长度为.4s,那么我希望多次测量,为程序测量的每个时间点以及它找到的值(频率)(以及可能的功率(dB))提供输出作为数组.复杂的是我想分析鸟歌,它们经常有谐波或信号超出频率范围(例如1000-2000赫兹).我希望程序也能输出这些信息,因为这对于我想对数据做的分析非常重要:)
现在有一段看起来非常像我想要的代码,但我认为它并没有给我所有我想要的价值....(感谢Justin Peel将这个问题发布到另一个问题:))所以我聚集在一起我需要numpy和pyaudio但不幸的是我不熟悉python所以我希望Python专家可以帮助我吗？
源代码:
# Read in a WAV and find the freq's
import pyaudio
import wave
import numpy as np
chunk = 2048
# open up a wave
wf = wave.open('test-tones/440hz.wav', 'rb')
swidth = wf.getsampwidth()
RATE = wf.getframerate()
# use a Blackman window
window = np.blackman(chunk)
# open stream
p = pyaudio.PyAudio()
stream = p.open(format =
p.get_format_from_width(wf.getsampwidth()),
channels = wf.getnchannels(),
rate = RATE,
output = True)
# read some data
data = wf.readframes(chunk)
# play stream and find the …