嘿,你能帮我,我怎么能让我的泽尼克时刻代码,旋转不变量。

2视图(30天)

显示旧的评论

tamanna 2013年4月10日

0
链接

这个问题直接联系

//www.tatmou.com/matlabcentral/answers/71507-hey-can-u-help-me-that-how-can-i-make-my-zernike-moment-code-invariant-to-rotation

回答: 阿伦anoop米2019年9月9日

4评论
显示隐藏 3年长的评论

tamanna 2013年4月11日

编辑:沃尔特·罗伯森 2013年4月11日

                               n = 4;
                              
                               m = 2
                              
                               p = rgb2gray (imread (“C: \ \惠普桌面\ \用户2. jpg));
                              
                               N =大小(p, 1);
                              
                               (ZOH AOH PhiOH] = Zernikmoment (p, n, m);
                              
                               函数φ[Z] = Zernikmoment (p, n,米)
                              
                               N =大小(p, 1);
                              
                               x = 1: N;y = x;
                              
                               (X, Y) = meshgrid (X, Y);
                              
                               R =√(2 * X-N-1)。^ 2 + (2 * Y-N-1)。^ 2) / N;
                              
                               θ=量化((N-1-2。* Y + 2), (2。* x n + 1 - 2));
                              
                               R = (R < = 1)。* R;
                              
                               Rad = radialpoly (R n m);%得到径向多项式
                              
                               产品= p (x, y)。* Rad。* exp(θ1我* *);
                              
                               Z =总和(产品(:));%计算时刻
                              
                               问= nnz (R) + 1;%计算像素的数量在单位圆
                              
                               Z = (n + 1) * Z /问;%正常化的振幅的时刻
                              
                               一个= abs (Z);%计算时刻的振幅
                              
                               φ=角(Z) * 180 /π;%计算mement的阶段
                              
                               函数rad = radialpoly (r n m)
                              
                               rad = 0(大小(r));% Initilization
                              
                               为s = 0: (n-abs (m)) / 2
                              
                               c = (1) ^ * ! (n) / (! (s) * ! ((n + abs (m)) / 2 s) * ! ((n-abs (m)) / 2 s));
                              
                               rad = rad + c * r。^ (n * s);
                              
                               结束