- c - 在位数组中找到第一个零
- linux - Unix 显示有关匹配两种模式之一的文件的信息
- 正则表达式替换多个文件
- linux - 隐藏来自 xtrace 的命令
我年轻的时候会编程一点,但我从来都不是很好。我发现 Python 非常适合我想做的事情。
我有一个 Excel 文件,其中包含我使用 numpy genfromtxt 函数读取的数据(64 列,18496 行)。我想将所有内容都放在一个名为 H 的 3D 矩阵中。我使用了三个循环来这样做,但我知道这不是最有效的。
data = np.genfromtxt(filename, delimiter = ";",skiprows = 11)
H = np.zeros((N,N,Nt))
for k in np.arange(N):
for l in np.arange(N):
for m in np.arange(Nt):
H[k,l,m] = data[m+Nt*k,l]
是否有切肉刀(计算速度更快)可以这样做。我想使用 numpy shape,但我做不到。
谢谢
最佳答案
您可以使用 np.reshape
reshape 形状& 然后用 np.transpose
重新排列维度, 像这样 -
H = data.reshape(N,Nt,N).transpose(0,2,1)
代替np.transpose
,也可以使用np.swapaxes
基本上我们在那里交换 axes 1,2
,就像这样 -
H = data.reshape(N,Nt,N).swapaxes(1,2)
sample 运行-
In [300]: N = 2
...: Nt = 3
...: data = np.random.randint(0,9,(N*Nt,N))
...:
In [301]: data
Out[301]:
array([[3, 6],
[7, 4],
[8, 1],
[8, 7],
[4, 8],
[2, 3]])
In [302]: H = np.zeros((N,N,Nt),dtype=data.dtype)
...: for k in np.arange(N):
...: for l in np.arange(N):
...: for m in np.arange(Nt):
...: H[k,l,m] = data[m+Nt*k,l]
...:
In [303]: H
Out[303]:
array([[[3, 7, 8],
[6, 4, 1]],
[[8, 4, 2],
[7, 8, 3]]])
In [304]: data.reshape(N,Nt,N).transpose(0,2,1)
Out[304]:
array([[[3, 7, 8],
[6, 4, 1]],
[[8, 4, 2],
[7, 8, 3]]])
运行时测试-
In [8]: # Input
...: N = 10
...: Nt = 10*50
...: data = np.random.randint(0,9,(N*Nt,N))
...:
...: def original_app(data):
...: H = np.zeros((N,N,Nt),dtype=data.dtype)
...: for k in np.arange(N):
...: for l in np.arange(N):
...: for m in np.arange(Nt):
...: H[k,l,m] = data[m+Nt*k,l]
...: return H
...:
In [9]: np.allclose(original_app(data),data.reshape(N,Nt,N).transpose(0,2,1))
Out[9]: True
In [10]: %timeit original_app(data)
10 loops, best of 3: 56.1 ms per loop
In [11]: %timeit data.reshape(N,Nt,N).transpose(0,2,1)
1000000 loops, best of 3: 1.25 µs per loop
关于python - 将 2D 数组复制到 3D 数组 - Python/NumPy,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/34576433/
我正在编写一个应用程序,允许用户创建一个“问卷”,然后向其中添加问题。我正在使用核心数据来存储信息。我创建了一个问卷实体,并与问题实体建立了“一对多”关系。我的问题是,如果要允许用户复制(复制)整个调
有没有办法复制或复制 SharedPreference?或者我需要从一个变量中获取每个变量,然后将它们放入另一个变量中吗? 最佳答案 尝试这样的事情: //sp1 is the shared pref
下面的(A)和(B)有区别吗? (假设 NON ARC,如果重要的话) // --- (A) --- @interface Zoo : NSObject{} @property (copy) Dog
我正在尝试将 mysql SELECT 查询保存到文件中,如下所示: $result = mysqli_query($db,$sql); $out = fopen('tmp/csv.csv', 'w'
我需要创建一个 CVPixelBufferRef 的副本,以便能够使用副本中的值以按位方式操作原始像素缓冲区。我似乎无法使用 CVPixelBufferCreate 或 CVPixelBufferCr
我在 Source 文件夹中有一个 Active wave 录音 wave-file.wav。我需要使用新名称 wave-file-copy.wav 将此文件复制到 Destination 文件夹。
在使用 GNU Autotools 构建的项目中,我有一个脚本需要通过 make 修改以包含安装路径。这是一个小例子: configure.ac: AC_INIT(foobar, 1.0) AC_PR
我想将 SQL 的行复制到同一个表中。但是在我的表中,我有一个“文本”列。 使用此 SQL: CREATE TEMPORARY TABLE produit2 ENGINE=MEMORY SELECT
谁能给我解释一下 df2 = df1 df2 = df1.copy() df3 = df1.copy(deep=False) 我已经尝试了所有选项并执行了以下操作: df1 = pd.DataFram
Hazelcast 是否具有类似于 Ehcache 的复制? http://www.ehcache.org/generated/2.9.0/pdf/Ehcache_Replication_Guide.
我有以下拓扑。一个 Ubuntu 16.04。运行我的全局 MySQL 服务器的 Amazon AWS 上的实例。我想将此服务器用作许多本地主服务器(Windows 机器 MySQL 服务器)的从服务
使用 SQLyog,我正在测试表中是否设置了正确的值。我尝试过 SELECT type_service FROM service WHERE email='test@gmail.com' 因此,只输出
有人可以提供一些关于如何配置 ElasticSearch 进行复制的说明。我在 Windows 中运行 ES,并且了解如果我在同一台服务器上多次运行 bat 文件,则会启动一个单独的 ES 实例,并且
一 点睛 ThreadGroup 复制线程的两个方法。 public int enumerate(Thread list[]) // 会将 ThreadGroup 中的 active 线程全部复制到
一 点睛 ThreadGroup 复制线程组的两个方法。 public int enumerate(ThreadGroup list[]) // 相对于 enumerate(list,true) pu
官方documentation Cassandra 说: Configure the keyspace and create the new datacenter: Use ALTER KEYSPAC
This question already has answers here: How to weight smoothing by arbitrary factor in ggplot2? (2个答
我们有一个表格来表明对各种俱乐部的兴趣。输出将数据记录在 Excel 电子表格中,其中列有他们的首选姓名、姓氏、电子邮件、代词,以及他们感兴趣的俱乐部的相应列中的“1”(下面的模型)。 我们希望为俱乐
This question already has answers here: Closed 8 years ago. Possible Duplicate: In vim, how do I get
如何复制形状及其所在的单元格?当我手动复制时,形状会跟随单元格,但是当我使用宏进行复制时,我会得到除形状之外的所有其他内容。 Cells(sourceRow, sourceColumn).Copy C
我是一名优秀的程序员,十分优秀!