invfrt_zp2.m

来自「beamlet变化的工具箱」· M 代码 · 共 39 行

function x = InvFRT_Zp2(WC)
% InvFRT_Zp2: Fast Inverse Ridgelet Transform for Zp2
% Usage
%	WC = InvFRT_Zp2(x)
% Inputs
%   WC  p by (p+1) array
% Outputs
%	x	p by p array, p prime
%
p = size(WC,1)-1;
S6 = MakeONFilter('Symmlet',6);
Rad = zeros(p,p+1);
Tmp = zeros(p+1,1);
for i=1:(p+1),
  Tmp = IWT_PO(WC(:,i),3,S6); % Emmanuel's Favorite
  Rad(:,i) = Tmp(1:p,1);
end
Ave = mean(Rad(:,1));
Rad(:,1) = Rad(:,1) - Ave;
x = InvRadon_Zp2_Fast(Rad);
x = x + Ave;

% Copyright (c) 2000 David L. Donoho
  

%
% Part of BeamLab Version:200
% Built:Friday,23-Aug-2002 00:00:00
% This is Copyrighted Material
% For Copying permissions see COPYING.m
% Comments? e-mail [e-mail]
%
%
% Part of BeamLab Version:200
% Built:Saturday,14-Sep-2002 00:00:00
% This is Copyrighted Material
% For Copying permissions see COPYING.m
% Comments? e-mail [e-mail]
%