基于傅里叶变换中kreisDemod函数里面所需要的m文件

function sd=tSPHT（c）

%SPHT spiral phase transform
% sd=SPHT（c） computes the quadrture term of c still affected by the
% direction phase factor. Therefore for a real c=b*cos（phi）
% sd=SPHT（c）=i*exp（i*dir）*b*sin（phi）
% Ref: Kieran G. Larkin Donald J. Bone and Michael A. Oldfield “Natural
% demodulation of two-dimensional fringe patterns. I. General background of the spiral phase quadrature %transform“ J. Opt. Soc. Am. A 18 1862-1870 （2001） 

%   AQ 19/8/09
%   Copyright 2009 OM4M
%   $ Revision: 1.0.0.0 $
%   $ Date: 19-08-2009 $

try
    
    TH=max（abs（c（:）））;
    if mean（real（c（:）））>0.01*TH
        warning（‘OM4M:SPHT:OutOfRange‘ ...
            ‘input must be DC filtered‘）;
    end
    
    [NR NC]=size（c）;
    [uv]