perform_warping.m

来自「该程序包为用MATLAB编写的一个实现小波变换的工具箱」· M 代码 · 共 17 行

function [t,u,v,ts,us,vs,I] = perform_warping(M,theta,P)

% sampling location
[Y,X] = meshgrid(1:P,1:P);
% projection on orthogonal direction
t = -sin(theta)*X(:) + cos(theta)*Y(:);
u = cos(theta)*X(:) + sin(theta)*Y(:);
% order points in increasing order
[t,I] = sort(t);
u = u(I);
% sorted signal
v = M(I);
% remove small entries
J = find( abs(v)>1e-2 );
vs = v(J);
ts = t(J);
us = u(J);