- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我需要找到一种优化的方法来构建表示 3D 平面的 2D numpy 数组(即存储其 Z 值)给定三个点(以数学方式定义平面)和所需的数组大小。
给定三个点 p0
、p1
、p2
和一个 size
变量:
import numpy as np
import pylab as plt
p0 = [0.2, -0.4, -0.2]
p1 = [-0.6, 0.1, -0.8]
p2 = [-0.1, 0.3, -0.6]
size = (12000, 17000)
下面的函数展示了如何从这 3 个点构建一个平面的想法:
def plane(p0, p1, p2, size):
p0 = np.array(p0)
p1 = np.array(p1)
p2 = np.array(p2)
ux, uy, uz = [p1[0] - p0[0], p1[1] - p0[1], p1[2] - p0[2]]
vx, vy, vz = [p2[0] - p0[0], p2[1] - p0[1], p2[2] - p0[2]]
uv = [uy * vz - uz * vy,
uz * vx - ux * vz,
ux * vy - uy * vx]
norm = np.array(uv)
d = -np.array(p0).dot(norm)
xx, yy = np.meshgrid(range(size[1]), range(size[0]))
z = np.array((-norm[0] * xx - norm[1] * yy - d) * 1. / norm[2])
return(z)
生成的 z
如下所示:
这个实现在几秒钟内吃掉了我的内存:
您将如何优化它以尽可能以最低的内存使用量获得完全相同的结果(快速也不是可选的...)?
注意:size
在两个维度上可能大 10 倍 -> 这显然会使我的 Python 在 2016 年主流笔记本电脑上崩溃。
(使用最常见的科学图书馆就可以了)
我在 Ubuntu linux 18.04.5 上使用 Python 3.6.9(默认,2021 年 1 月 26 日,15:33:00)[GCC 8.4.0]。
最佳答案
直到计算d
,你的数组很小。之后的一切都是为操作的每个部分创建一个单独的 1.5GiB (12k * 17k * 8bytes/float) 临时数组。您可以完全摆脱这些数组中的大部分,并就地完成其余的工作。例如:
xx, yy = np.meshgrid(...)
您真的不需要 meshgrid:在 np.arange(size[0])[:, None]
上广播和 np.arange(size[1])
将制作您需要的所有阵列。对于那部分数据,这会将您的大小从 ~3GiB 减少到 ~230KiB:
xx = np.arange(size[1], dtype=float)
yy = np.arange(size[0], dtype=float)[:, None]
表达式z = ...
其中大约有七个临时数组,每个数组的大小为 1.5GiB:
t1 = -norm[0] * xx
t2 = norm[1] * yy
t3 = t1 - t2
t4 = t3 - d
t5 = t4 * 1
t6 = t5 / norm[2]
t7 = np.array(t6)
其中一些数组可能会在计算过程中被垃圾回收,但它们都需要在某个时候进行分配。您可以使用广播和 in-place operators 摆脱所有这些或 out
参数到适当的 ufunc
小号:
所以:
d /= norm[2]
norm /= norm[2]
xx *= -norm[0]
yy *= -norm[1]
yy -= d
z = xx + yy
即使在代码的紧凑部分,你也在做一堆不必要的操作。例如,您调用 np.array(p0)
两次。这是您的函数的更好版本:
def plane2(p0, p1, p2, size):
u = np.subtract(p1, p0)
v = np.subtract(p2, p0)
uxv = np.cross(u, v)
d = -np.dot(p0, uxv)
d /= uxv[2]
uxv /= uxv[2]
xx = np.arange(size[1], dtype=float)
xx *= -uxv[0]
yy = np.arange(size[0], dtype=float)[:, None]
yy *= -uxv[1]
yy -= d
return xx + yy
对于比建议尺寸小 10 倍的尺寸:
%timeit plane(p0, p1, p2, size)
70.1 ms ± 329 µs per loop (mean ± std. dev. of 7 runs, 10 loops each)
%timeit plane2(p0, p1, p2, size)
8.29 ms ± 302 µs per loop (mean ± std. dev. of 7 runs, 100 loops each)
当我将大小增加到您问题中的大小时,您的实现崩溃,而我的实现在合理的时间内完成:
%timeit plane(p0, p1, p2, size)
...
MemoryError: Unable to allocate 1.52 GiB for an array with shape (12000, 17000) and data type float64
%timeit plane2(p0, p1, p2, size)
819 ms ± 37.6 ms per loop (mean ± std. dev. of 7 runs, 1 loop each)
根据您拥有多少 RAM 以及 python 进程中空闲内存的设置,无论您执行操作的效率如何,您都可能无法将数组大小增加 100 倍至 152GiB。
关于python - 构建一个表示 3D 平面(存储其 Z 值)的 2D 数组,由 3 个点和所需的数组大小定义,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/68277692/
我的问题:非常具体。我正在尝试想出解析以下文本的最简单方法: ^^domain=domain_value^^version=version_value^^account_type=account_ty
好吧,这就是我的困境: 我正在为 Reddit 子版 block 开发常见问题解答机器人。我在 bool 逻辑方面遇到了麻烦,需要一双更有经验的眼睛(这是我在 Python 中的第一次冒险)。现在,该
它首先遍历所有 y 值,然后遍历所有 x 值。我需要 X 和 y 同时改变。 For x = 3 To lr + 1 For y = 2 To lr anyl.Cells(x, 1)
假设我有一个包含 2 列的 Excel 表格:单元格 A1 到 A10 中的日期和 B1 到 B10 中的值。 我想对五月日期的所有值求和。我有3种可能性: {=SUM((MONTH(A1:A10)=
如何转换 Z-score来自 Z-distribution (standard normal distribution, Gaussian distribution)到 p-value ?我还没有找到
我正在重写一些 Javascript 代码以在 Excel VBA 中工作。由于在这个网站上搜索,我已经设法翻译了几乎所有的 Javascript 代码!但是,有些代码我无法准确理解它在做什么。这是一
我遇到过包含日期格式的时间戳日期的情况。然后我想构建一个图表,显示“点击”项目的数量“每天”, //array declaration $array1 = array("Date" => 0); $a
我是scala的新手! 我的问题是,是否有包含成员的案例类 myItem:Option[String] 当我构造类时,我需要将字符串内容包装在: Option("some string") 要么 So
我正在用 PHP 创建一个登录系统。我需要用户使用他或她的用户名或电子邮件或电话号码登录然后使用密码。因为我知道在 Java 中我们会像 email==user^ username == user 这
我在 C++ 项目上使用 sqlite,但是当我在具有文本值的列上使用 WHERE 时出现问题 我创建了一个 sqlite 数据库: CREATE TABLE User( id INTEGER
当构造函数是显式时,它不用于隐式转换。在给定的代码片段中,构造函数被标记为 explicit。那为什么在 foo obj1(10.25); 情况下它可以工作,而在 foo obj2=10.25; 情况
我知道这是一个主观问题,所以如果需要关闭它,我深表歉意,但我觉得它经常出现,让我想知道是否普遍偏爱一种形式而不是另一种形式。 显然,最好的答案是“重构代码,这样你就不需要测试是否存在错误”,但有时没有
这两个 jQuery 选择器有什么区别? 以下是来自 w3schools.com 的定义: [attribute~=value] 选择器选择带有特定属性,其值包含特定字符串。 [attribute*=
为什么我们需要CSS [attribute|=value] Selector根本当 CSS3 [attribute*=value] Selector基本上完成相同的事情,浏览器兼容性几乎相似?是否存在
我正在解决 regx 问题。我已经有一个像这样的 regx [0-9]*([.][0-9]{2})。这是 amont 格式验证。现在,通过此验证,我想包括不应提供 0 金额。比如 10 是有效的,但
我正在研究计算机科学 A 考试的样题,但无法弄清楚为什么以下问题的正确答案是正确的。 考虑以下方法。 public static void mystery(List nums) { for (
好的,我正在编写一个 Perl 程序,它有一个我收集的值的哈希值(完全在一个完全独立的程序中)并提供给这个 Perl 脚本。这个散列是 (string,string) 的散列。 我想通过 3 种方式对
我有一个表数据如下,来自不同的表。仅当第三列具有值“债务”并且第一列(日期)具有最大值时,我才想从第四列中获取最大值。最终值基于 MAX(DATE) 而不是 MAX(PRICE)。所以用简单的语言来说
我有一个奇怪的情况,只有错误状态保存到数据库中。当“状态”应该为 true 时,我的查询仍然执行 false。 我有具有此功能的 Controller public function change_a
我有一个交易表(针对所需列进行了简化): id client_id value 1 1 200 2 2 150 3 1
我是一名优秀的程序员,十分优秀!