- android - 多次调用 OnPrimaryClipChangedListener
- android - 无法更新 RecyclerView 中的 TextView 字段
- android.database.CursorIndexOutOfBoundsException : Index 0 requested, 光标大小为 0
- android - 使用 AppCompat 时,我们是否需要明确指定其 UI 组件(Spinner、EditText)颜色
我正在尝试设置以下非线性成本函数 (1) 和 (2)
通过使用 numpy einsum 函数。
我尝试将它们转换为 python,其中 R
是形状 (14,3,3) 和 v
、vt
(表示等式 (1) 中的 v') 的形状为 (14,3)。 old_points
和 new_points
分别代表 p
和 p'
,其形状为 (6890,3)。
def costfunc(x, old_points, new_points, weights, n_joints):
"""
Set up non-linear cost functions by using equations from LBS:
(1) p'_i = sum{j}(w_ji (R_j p_i + v_j))
(2) p_i - sum{j}(w_ji (R_j p'_i + v'_j))
where Rt denotes the transpose of R.
:param old_points: original vertex positions
:param new_points: transformed vertex positions
:param weights: weight matrix obtained from spectral clustering
:param n_joints: number of joints
:return: non-linear cost functions as in (1), (2) to find the root of
"""
# Extract rotations R, Rt and offsets v, v' from rv
R = np.array([(np.array(x[j * 15:j * 15 + 9]).reshape(3, 3)) for j in range(n_joints)])
Rt = np.array([R[j].T for j in range(n_joints)])
v = np.array([(np.array(x[j * 15 + 9:j * 15 + 12])) for j in range(n_joints)])
vt = np.array([(np.array(x[j * 15 + 12:j * 15 + 15])) for j in range(n_joints)])
## Use equations (1) and (2) for the non-linear pass.
# R_j p_i
Rp = np.einsum('jkl,il', R, old_points)
Rtv = np.einsum('jkl,il', Rt, v)
# Rt_j p'_i
Rtp = np.einsum('jkl,il', Rt, new_points)
Rvt = np.einsum('jkl,il', R, vt)
# w_ji (Rp_ij - Rtv_j)
wRpv = np.einsum('ji,ijk->ik', weights, Rp - Rvt)
# w_ji (Rtp'_ij - Rv'_j)
wRtpv = np.einsum('ji,ijk->ik', weights, Rtp - Rtv)
# Set up a non-linear cost function, then compute the squared norm.
d = new_points - wRpv
dt = old_points - wRtpv
norm = np.linalg.norm(d, axis=1)
normt = np.linalg.norm(dt, axis=1)
result = np.concatenate([norm, normt])
return np.power(result, 2)
现在,计算 wRpv
和 wRtpv
的行中存在错误ValueError:操作数无法与形状 (6890,14,3) (14,14,3) 一起广播
。我该如何解决这个问题?非常感谢任何帮助!
最佳答案
我现在明白了。这是解决方案:
def costfunc(x, old_points, new_points, weights, n_joints):
"""
Set up non-linear cost functions by using equations from LBS:
(1) p'_i = sum{j}(w_ji (R_j p_i + v_j))
(2) p_i - sum{j}(w_ji (R_j p'_i + v'_j))
where Rt denotes the transpose of R.
:param old_points: original vertex positions
:param new_points: transformed vertex positions
:param weights: weight matrix obtained from spectral clustering
:param n_joints: number of joints
:return: non-linear cost functions as in (1), (2) to find the root of
"""
# Extract rotations R, Rt and offsets v, v' from rv
R = np.array([(np.array(x[j * 15:j * 15 + 9]).reshape(3, 3)) for j in range(n_joints)])
Rt = np.array([R[j].T for j in range(n_joints)])
v = np.array([(np.array(x[j * 15 + 9:j * 15 + 12])) for j in range(n_joints)])
vt = np.array([(np.array(x[j * 15 + 12:j * 15 + 15])) for j in range(n_joints)])
## Use equations (1) and (2) for the non-linear pass.
# R_j p_i
Rp = np.einsum('jkl,il', R, old_points)
# Rt_j p'_i
Rtp = np.einsum('jkl,il', Rt, new_points)
# R_j v'_j
Rvt = np.array([R[i] @ vt[i] for i in range(n_joints)])
# Rt_j v_j
Rtv = np.array([Rt[i] @ v[i] for i in range(n_joints)])
# w_ji (Rp_ij - Rtv_j)
wRpv = np.einsum('ji,ijk->ik', weights, Rp - Rvt)
# w_ji (Rtp'_ij - Rv'_j)
wRtpv = np.einsum('ji,ijk->ik', weights, Rtp - Rtv)
# Set up a non-linear cost function, then compute the squared norm.
d = new_points - wRpv
dt = old_points - wRtpv
norm = np.linalg.norm(d, axis=1)
normt = np.linalg.norm(dt, axis=1)
result = np.concatenate([norm, normt])
return np.power(result, 2)
关于python - 非线性函数的 numpy 爱因斯坦求和,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/57017071/
我正在尝试在 OCaml 中创建一个函数,该函数在数学中执行求和函数。 我试过这个: sum n m f = if n = 0 then 0 else if n > m then f
我正在尝试找到一个可以帮助我解决问题的公式。 这个公式应该对每个靠近(总是在左边)具有相同名称的单元格的单元格求和(或工作)。如下所示: 将每个大写字母视为 “食谱”并且每个小写字母为 “成分” .在
让它成为以下 python pandas DataFrame,其中每一行代表一个人在酒店的住宿。 | entry_date | exit_date | days | other_columns
我有显示客户来电的数据。我有客户号码、电话号码(1 个客户可以有多个)、每个语音调用的日期记录以及调用持续时间的列。表看起来如下示例。 CusID | PhoneNum | Date
让它成为以下 python pandas DataFrame,其中每一行代表一个人在酒店的住宿。 | entry_date | exit_date | days | other_columns
我得到了两列数据; 答: 2013年12月31日 2013年12月30日 2013年12月29日 2013年12月28日 2013年12月27日 2012年12月26日 B: 10 10 10 10
我对 double 格式的精度有疑问。 示例: double K=0, L=0, M=0; scanf("%lf %lf %lf", &K, &L, &M); if((K+L) 我的测试输入: K
我有以下数组: int[,] myArray1 = new int[2, 3] { { 1, 2, 3 }, { 4, 6, 8 } }; int[,] myArray2 = new int[2, 3
我需要有关报告查询的帮助。我在该方案的底部有一个发票表,需要一种方法来获取总计费金额,同时在此数据库方案中的较高点进行条件过滤。我需要加入其他表,这会导致 SUM 函数返回不需要的结果。 这是我正在使
我有一个使用innodb作为存储引擎的MySQL数据库,并且我有许多采用基本形式的查询: SELECT bd.billing, SUM(CASE WHEN tc.transaction_class
尝试创建一个查询来给出总胜、平和负。我有以下查询 SELECT CASE WHEN m.home_team = '192' AND m.home_full_time_score
我正在尝试生成一份报告,显示排名靠前的推荐人以及他们推荐的人产生了多少收入。 这是我的表格的缩写版本: Users Table ------------------ id referral_user_
我有以下查询,并得到了预期的结果: SELECT IF (a1>b1,'1','0') AS a1r, IF (a2>b2,'1','0') AS a2r,
我尝试了几种不同的解决方案,但都没有成功。我给出的表格是一个示例,其设计和功能与我实际使用的表格类似: PK | Color | Count -------------------
我正在尝试构建一个查询来检查我的库存。 SELECT COUNT(*) AS item_count, reseller_id, sum(sold) as sold_count, sum(refunde
我试图解决一个看起来像下面编写的代码的问题,但由于缺乏知识和阅读 sqlalchemy 文档,我还没有真正找到解决问题的方法。 目标: 如果 year_column 中的年份相同,则获取 sales_
我有一个包含一周中多天的表格。一周中的每一天都有独特的属性,例如冰淇淋是否在这一天成功送达: ID DAY_WEEK ICE_CREAM 1 Monday
首先,我有一个名为store_00的表 id | ref | item | qty | cost | sell 1 22 x1 5 10 15 2 22
我正在编写一个程序,计算每个数字的总和,直到 1000。例如,1+2+3+4+5....+100。首先,我将求和作业分配给 10 个处理器:处理器 0 得到 1-100,处理器 1 得到 101-20
我想在一个循环中一次对多个属性求和: class Some(object): def __init__(self, acounter, bcounter): self.acou
我是一名优秀的程序员,十分优秀!