- iOS/Objective-C 元类和类别
- objective-c - -1001 错误,当 NSURLSession 通过 httpproxy 和/etc/hosts
- java - 使用网络类获取 url 地址
- ios - 推送通知中不播放声音
我正在编写一个程序来查找闭合平面曲线的椭圆傅里叶系数 (EFC),但在系数归一化方面遇到了问题。
闭合平面折线 p 由点集 m_points 描述:m_points[i][0] 保持 xi 坐标,m_points i保持 yi 坐标。我从 0 开始到 m_num_points-1。
折线的 EFC 由该算法计算(系数在 EFD 数组中)。
// calc accumulated curve length, delta x, and delta y
dt[0] := 0;
dx[0] := 0;
dy[0] := 0;
tp[0] := 0;
T := 0;
for i:=1 to p.m_num_points-1 do begin
va := VectorAffineSubtract(p.m_points[i], p.m_points[i-1]);
dt[i] := VectorLength( va );
dx[i] := p.m_points[i][0] - p.m_points[i-1][0];
dy[i] := p.m_points[i][1] - p.m_points[i-1][1];
tp[i] := tp[i-1] + dt[i];
T := tp[i];
end;
Tpi := T / (2*PI*PI);
piT := PI*2 / T;
// find elliptic fourier coefficients
// first
An := 0;
Cn := 0;
for i:=0 to p.m_num_points-1 do begin
An := An + (p.m_points[i][0] + p.m_points[i-1][0]) * dt[i] / 2;
Cn := Cn + (p.m_points[i][1] + p.m_points[i-1][1]) * dt[i] / 2;
end;
EFD[0][0] := An / T;
EFD[0][1] := 0;
EFD[0][2] := Cn / T;
EFD[0][3] := 0;
// next
for n:=1 to m_num_efd do begin
Tn := Tpi / (n*n);
piTn := piT * n;
An := 0;
Bn := 0;
Cn := 0;
Dn := 0;
for i:=1 to p.m_num_points-1 do begin
An := An + dx[i]/dt[i]*( cos(piTn*tp[i]) - cos(piTn*tp[i-1]) );
Bn := Bn + dx[i]/dt[i]*( sin(piTn*tp[i]) - sin(piTn*tp[i-1]) );
Cn := Cn + dy[i]/dt[i]*( cos(piTn*tp[i]) - cos(piTn*tp[i-1]) );
Dn := Dn + dy[i]/dt[i]*( sin(piTn*tp[i]) - sin(piTn*tp[i-1]) );
end;
EFD[n][0] := An * Tn;
EFD[n][1] := Bn * Tn;
EFD[n][2] := Cn * Tn;
EFD[n][3] := Dn * Tn;
end;
为了从谐波集中恢复折线,我使用了该算法。
// restore outline
resP := TYPolyline.create();
for i:=0 to p.m_num_points-1 do begin
xn := EFD[0][0];
yn := EFD[0][2];
for n:=1 to m_num_efd do begin
xn := xn + EFD[n][0] * cos(piT*n*tp[i]) + EFD[n][1] * sin(piT*n*tp[i]);
yn := yn + EFD[n][2] * cos(piT*n*tp[i]) + EFD[n][3] * sin(piT*n*tp[i]);
end;
resP.add_point( xn, yn );
end;
现在我需要规范化 EFD 并将它们放在一个新数组 NEFD 中。我这样做:
// Normalize EFD
An := EFD[0][0];
Bn := EFD[0][1];
Cn := EFD[0][2];
Dn := EFD[0][3];
for n:=0 to m_num_efd do begin
NEFD[n][0] := EFD[n][0];
NEFD[n][1] := EFD[n][1];
NEFD[n][2] := EFD[n][2];
NEFD[n][3] := EFD[n][3];
end;
// rotate starting point
angl_o := 0.5 * ArcTan2( 2*(An*Bn + Cn*Dn), An*An + Cn*Cn - Bn*Bn - Dn*Dn );
for n:=1 to m_num_efd+1 do begin
NEFD[n-1][0]:= EFD[n-1][0]*cos((n)*angl_o) + EFD[n-1][1]*sin((n)*angl_o);
NEFD[n-1][1]:=-EFD[n-1][0]*sin((n)*angl_o) + EFD[n-1][1]*cos((n)*angl_o);
NEFD[n-1][2]:= EFD[n-1][2]*cos((n)*angl_o) + EFD[n-1][3]*sin((n)*angl_o);
NEFD[n-1][3]:=-EFD[n-1][3]*sin((n)*angl_o) + EFD[n-1][3]*cos((n)*angl_o);
end;
// make the semi-major axis parallel to the x-axis
angl_w := ArcTan( NEFD[1][2] / NEFD[1][0] );
cs := cos(angl_w);
ss := sin(angl_w);
for n:=1 to m_num_efd do begin
NEFD[n][0] := cs*NEFD[n][0] + ss*NEFD[n][2];
NEFD[n][1] := cs*NEFD[n][1] + ss*NEFD[n][3];
NEFD[n][2] :=-ss*NEFD[n][0] + cs*NEFD[n][2];
NEFD[n][3] :=-ss*NEFD[n][1] + cs*NEFD[n][3];
end;
// size invariant
R := sqrt( NEFD[1][0]*NEFD[1][0] );
for n:=0 to m_num_efd do begin
NEFD[n][0] := NEFD[n][0] / R;
NEFD[n][1] := NEFD[n][1] / R;
NEFD[n][2] := NEFD[n][2] / R;
NEFD[n][3] := NEFD[n][3] / R;
end;
当我尝试使半长轴平行于 X 并旋转系数时,我得到了一个丑陋的结果。看起来恢复后的形状是围绕 z 轴旋转的。 (左侧为源形状,右侧为还原形状。)
我的代码有什么问题?
更新
@MBo 成功回答后,需要对代码进行以下修改:
// MODIFIED: change n-1 to n
// rotate starting point
angl_o := 0.5 * ArcTan2( 2*(An*Bn + Cn*Dn), An*An + Cn*Cn - Bn*Bn - Dn*Dn );
for n:=1 to m_num_efd+1 do begin
NEFD[n][0]:= EFD[n][0]*cos((n)*angl_o) + EFD[n][1]*sin((n)*angl_o);
NEFD[n][1]:=-EFD[n][0]*sin((n)*angl_o) + EFD[n][1]*cos((n)*angl_o);
NEFD[n][2]:= EFD[n][2]*cos((n)*angl_o) + EFD[n][3]*sin((n)*angl_o);
NEFD[n][3]:=-EFD[n][3]*sin((n)*angl_o) + EFD[n][3]*cos((n)*angl_o);
end;
// MODIFIED: change left NEFD to EFD
// make the semi-major axis parallel to the x-axis
angl_w := ArcTan( NEFD[1][2] / NEFD[1][0] );
cs := cos(angl_w);
ss := sin(angl_w);
for n:=1 to m_num_efd do begin
EFD[n][0] := cs*NEFD[n][0] + ss*NEFD[n][2];
EFD[n][1] := cs*NEFD[n][1] + ss*NEFD[n][3];
EFD[n][2] :=-ss*NEFD[n][0] + cs*NEFD[n][2];
EFD[n][3] :=-ss*NEFD[n][1] + cs*NEFD[n][3];
end;
// ADDED: place normalized EFD into NEFD
for n:=0 to m_num_efd do begin
NEFD[n][0] := EFD[n][0];
NEFD[n][1] := EFD[n][1];
NEFD[n][2] := EFD[n][2];
NEFD[n][3] := EFD[n][3];
end;
最佳答案
可能原因:
在此循环中,您使用NEFD[n][0]
的修改 值来计算NEFD[n][2]
的新值( NEFD[n][1]
相同)。看来您必须存储和使用未修改的值。
for n:=1 to m_num_efd do begin
NEFD[n][0] := cs*NEFD[n][0] + ss*NEFD[n][2];
NEFD[n][1] := cs*NEFD[n][1] + ss*NEFD[n][3];
vvvvvvvvv
NEFD[n][2] :=-ss*NEFD[n][0] + cs*NEFD[n][2];
vvvvvvvvv
NEFD[n][3] :=-ss*NEFD[n][1] + cs*NEFD[n][3];
end;
关于algorithm - 如何归一化椭圆傅里叶系数?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/36710928/
我在一本书(Interview Question)中读到这个问题,想在这里详细讨论这个问题。请点亮它。 问题如下:- 隐私和匿名化 马萨诸塞州集团保险委员会早在 1990 年代中期就有一个绝妙的主意
我最近接受了一次面试,面试官给了我一些伪代码并提出了相关问题。不幸的是,由于准备不足,我无法回答他的问题。由于时间关系,我无法向他请教该问题的解决方案。如果有人可以指导我并帮助我理解问题,以便我可以改
这是我的代码 public int getDist(Node root, int value) { if (root == null && value !=0) return
就效率而言,Strassen 算法应该停止递归并应用乘法的最佳交叉点是多少? 我知道这与具体的实现和硬件密切相关,但对于一般情况应该有某种指南或某人的一些实验结果。 在网上搜索了一下,问了一些他们认为
我想学习一些关于分布式算法的知识,所以我正在寻找任何书籍推荐。我对理论书籍更感兴趣,因为实现只是个人喜好问题(我可能会使用 erlang(或 c#))。但另一方面,我不想对算法进行原始的数学分析。只是
我想知道你们中有多少人实现了计算机科学的“ classical algorithms ”,例如 Dijkstra's algorithm或现实世界中的数据结构(例如二叉搜索树),而不是学术项目? 当有
我正在解决旧编程竞赛中的一些示例问题。在这个问题中,我们得到了我们有多少调酒师以及他们知道哪些食谱的信息。制作每杯鸡尾酒需要 1 分钟,我们需要使用所有调酒师计算是否可以在 5 分钟内完成订单。 解决
关闭。这个问题是opinion-based .它目前不接受答案。 想要改进这个问题? 更新问题,以便 editing this post 可以用事实和引用来回答它. 关闭 8 年前。 Improve
我开始学习 Nodejs,但我被困在中间的某个地方。我从 npm 安装了一个新库,它是 express -jwt ,它在运行后显示某种错误。附上代码和错误日志,请帮助我! const jwt = re
我有一个证书,其中签名算法显示“sha256rsa”,但指纹算法显示“sha1”。我的证书 SHA1/SHA2 的标识是什么? 谢谢! 最佳答案 TL;TR:签名和指纹是完全不同的东西。对于证书的强度
我目前在我的大学学习数据结构类(class),并且在之前的类(class)中做过一些算法分析,但这是我在之前的类(class)中遇到的最困难的部分。我们现在将在我的数据结构类(class)中学习算法分
有一个由 N 个 1x1 方格组成的区域,并且该区域的所有部分都是相连的(没有任何方格无法到达的方格)。 下面是一些面积的例子。 我想在这个区域中选择一些方块,并且两个相邻的方块不能一起选择(对角接触
我有一些多边形形状的点列表,我想将其包含在我页面上的 Google map 中。 我已经从原始数据中删除了尽可能多的不必要的多边形,现在我剩下大约 12 个,但它们非常详细以至于导致了问题。现在我的文
我目前正在实现 Marching Squares用于计算等高线曲线,我对此处提到的位移位的使用有疑问 Compose the 4 bits at the corners of the cell to
我正在尝试针对给定算法的约束满足问题实现此递归回溯函数: function BACKTRACKING-SEARCH(csp) returns solution/failure return R
是否有包含反函数的库? 作为项目的一部分,我目前正在研究测向算法。我正在使用巴特利特相关性。在 Bartlett 相关性中,我需要将已经是 3 次矩阵乘法(包括 Hermitian 转置)的分子除以作
关闭。这个问题不符合Stack Overflow guidelines .它目前不接受答案。 这个问题似乎与 help center 中定义的范围内的编程无关。 . 关闭 8 年前。 Improve
问题的链接是UVA - 1394 : And There Was One . 朴素的算法是扫描整个数组并在每次迭代中标记第 k 个元素并在最后停止:这需要 O(n^2) 时间。 我搜索了一种替代算法并
COM 中创建 GUID 的函数 (CoCreateGUID) 使用“分散唯一性算法”,但我的问题是,它是什么? 谁能解释一下? 最佳答案 一种生成 ID 的方法,该 ID 具有一定的唯一性保证,而不
在做一个项目时我遇到了这个问题,我将在这个问题的实际领域之外重新措辞(我想我可以谈论烟花的口径和形状,但这会使理解更加复杂).我正在寻找一种(可能是近似的)算法来解决它。 我有 n 个不同大小的容器,
我是一名优秀的程序员,十分优秀!