MATLAB /matlab_tools/Converted/kkwt.m

Language MATLAB Lines 237
MD5 Hash b93af65935117e3f4ee6f4278bc28b05
Repository git://github.com/aludnam/MATLAB.git View Raw File View Project SPDX
  1
  2
  3
  4
  5
  6
  7
  8
  9
 10
 11
 12
 13
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26
 27
 28
 29
 30
 31
 32
 33
 34
 35
 36
 37
 38
 39
 40
 41
 42
 43
 44
 45
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
%kkwt '5D wavelet transform'
% This MatLab function was automatically generated by a converter (KhorosToMatLab) from the Khoros kwt.pane file
%
% Parameters: 
% InputFile: i 'Input Object ', required: 'Input object filename '
% InputFile: fh 'Highpass', required: 'Input highpass filter coefficients'
% InputFile: fl 'Lowpass ', required: 'Input lowpass filter coefficients'
% OutputFile: o 'Output Object ', required: 'Output object filename'
% Toggle: w 'Width', default: 0: 'filter in the WIDTH direction'
% Toggle: h 'Height', default: 0: 'filter in the HEIGHT direction'
% Toggle: d 'Depth', default: 0: 'filter in the DEPTH direction'
% Toggle: t 'Time', default: 0: 'filter in the TIME direction'
% Toggle: e 'Elements', default: 0: 'filter in the ELEMENTS direction'
%
% Example: o = kkwt({i, fh, fl}, {'i','';'fh','';'fl','';'o','';'w',0;'h',0;'d',0;'t',0;'e',0})
%
% Khoros helpfile follows below:
%
%  PROGRAM
% kwt - Discrete Wavelet Transform (Forward and Inverse)
%
%  DESCRIPTION
% "kwt" performs the wavelet decomposition of the
% source object defined by \fBInput" in any dimension or a combination of the 5
% dimensions.  The filters for both transform directions (forward and
% inverse) are defined by \fBHighpass" for the highpass filter and
% \fBLowpass" for the lowpass filter.  These filter must be defined
% acordlingly to the transform direction.  The output object of "kwt"
% corresponds to the coefficients of the wavelet expansion.  Default
% values are: width and height dimensions and transform direction forward.
% NOTE: Filters must be 1D.
% 
% "kwt" uses tensorial product to perform the multidimensional
% wavelet transform.  That is, the wavelet transform routine
% assumes separability (decoupling) in the multidimensional case.
% 
%  "Data Type" 5
% The "kwt" routine accepts any kind of data type for the source 
% objects, except complex.  Internally, "kwt" casts the data of
% every source object to double before processing and changes the
% destination object data type (output) accordingly.  The data is
% processed using double despite the data type of any source object.  
% Note that data is not casted to a lower type for processing.
% 
%  "Map Data" 5
% If any source object contains map data as well as value data, then
% "kwt" removes the map.  The destination object is always unmapped.
% 
%  "Mask, Time and Location Segments" 5
% If any source object contains a mask, time or location segment, then an
% error message is printed out.  Mask, time and location segments have to
% be removed before processing.
%
%  
%
%  EXAMPLES
% 2D forward wavelet transform of a 2D object using daub_2_f_h (highpass 
% Daubechies filter, order 2, forward direction) as a highpass filter and
% daub_2_f_l (lowpass Daubechies filter, order 2, forward direction) as a
% lowpass filter.  
% kwt -i image_input -w -h -fh daub_2_f_h -fl daub_2_f_l -dir 0 -o image_output
% 2D inverse wavelet transform of a 2D object using daub_2_i_h (highpass 
% Daubechies filter, order 2, inverse direction) as a highpass filter and
% daub_2_i_l (lowpass Daubechies filter, order 2, inverse direction) as a
% lowpass filter.  
% kwt -i image_input -w -h -fh daub_2_i_h -fl daub_2_i_l -dir 1 -o image_output
%
%  "SEE ALSO"
% kmwt
%
%  RESTRICTIONS 
% Dimension of data record must be a power of 2.
% Filters must be 1D.
%
%  REFERENCES 
% 1) WAVELETS ANALYSIS AND ITS APPLICATIONS, Vol. I, II and III, by
% Charles K. Chui, Academis Press, 1992.
% 2) ADAPTED WAVELET ANALYSIS FROM THEORY TO SOFTWARE,
% Mladen Victor Wickerhauser, A K Peters, 1994.
% 3) TEN LECTURES ON WAVELETS, I. Daubechies, Siam, 1992.
% 4) DIFFERENT PERSPECTIVAS ON WAVELETS, Proceedings of Symposia in
% Applied Mathematics, 1993.
%
%  COPYRIGHT
% Copyright (C) 1993 - 1997, Khoral Research, Inc.  All rights reserved.
% 


function varargout = kkwt(varargin)
if nargin ==0
  Inputs={};arglist={'',''};
elseif nargin ==1
  Inputs=varargin{1};arglist={'',''};
elseif nargin ==2
  Inputs=varargin{1}; arglist=varargin{2};
else error('Usage: [out1,..] = kkwt(Inputs,arglist).');
end
if size(arglist,2)~=2
  error('arglist must be of form {''ParameterTag1'',value1;''ParameterTag2'',value2}')
 end
narglist={'i', '__input';'fh', '__input';'fl', '__input';'o', '__output';'w', 0;'h', 0;'d', 0;'t', 0;'e', 0};
maxval={0,0,0,0,0,0,0,0,0};
minval={0,0,0,0,0,0,0,0,0};
istoggle=[0,0,0,0,1,1,1,1,1];
was_set=istoggle * 0;
paramtype={'InputFile','InputFile','InputFile','OutputFile','Toggle','Toggle','Toggle','Toggle','Toggle'};
% identify the input arrays and assign them to the arguments as stated by the user
if ~iscell(Inputs)
Inputs = {Inputs};
end
NumReqOutputs=1; nextinput=1; nextoutput=1;
  for ii=1:size(arglist,1)
  wasmatched=0;
  for jj=1:size(narglist,1)
   if strcmp(arglist{ii,1},narglist{jj,1})  % a given argument was matched to the possible arguments
     wasmatched = 1;
     was_set(jj) = 1;
     if strcmp(narglist{jj,2}, '__input')
      if (nextinput > length(Inputs)) 
        error(['Input ' narglist{jj,1} ' has no corresponding input!']); 
      end
      narglist{jj,2} = 'OK_in';
      nextinput = nextinput + 1;
     elseif strcmp(narglist{jj,2}, '__output')
      if (nextoutput > nargout) 
        error(['Output nr. ' narglist{jj,1} ' is not present in the assignment list of outputs !']); 
      end
      if (isempty(arglist{ii,2}))
        narglist{jj,2} = 'OK_out';
      else
        narglist{jj,2} = arglist{ii,2};
      end

      nextoutput = nextoutput + 1;
      if (minval{jj} == 0)  
         NumReqOutputs = NumReqOutputs - 1;
      end
     elseif isstr(arglist{ii,2})
      narglist{jj,2} = arglist{ii,2};
     else
        if strcmp(paramtype{jj}, 'Integer') & (round(arglist{ii,2}) ~= arglist{ii,2})
            error(['Argument ' arglist{ii,1} ' is of integer type but non-integer number ' arglist{ii,2} ' was supplied']);
        end
        if (minval{jj} ~= 0 | maxval{jj} ~= 0)
          if (minval{jj} == 1 & maxval{jj} == 1 & arglist{ii,2} < 0)
            error(['Argument ' arglist{ii,1} ' must be bigger or equal to zero!']);
          elseif (minval{jj} == -1 & maxval{jj} == -1 & arglist{ii,2} > 0)
            error(['Argument ' arglist{ii,1} ' must be smaller or equal to zero!']);
          elseif (minval{jj} == 2 & maxval{jj} == 2 & arglist{ii,2} <= 0)
            error(['Argument ' arglist{ii,1} ' must be bigger than zero!']);
          elseif (minval{jj} == -2 & maxval{jj} == -2 & arglist{ii,2} >= 0)
            error(['Argument ' arglist{ii,1} ' must be smaller than zero!']);
          elseif (minval{jj} ~= maxval{jj} & arglist{ii,2} < minval{jj})
            error(['Argument ' arglist{ii,1} ' must be bigger than ' num2str(minval{jj})]);
          elseif (minval{jj} ~= maxval{jj} & arglist{ii,2} > maxval{jj})
            error(['Argument ' arglist{ii,1} ' must be smaller than ' num2str(maxval{jj})]);
          end
        end
     end
     if ~strcmp(narglist{jj,2},'OK_out') &  ~strcmp(narglist{jj,2},'OK_in') 
       narglist{jj,2} = arglist{ii,2};
     end
   end
   end
   if (wasmatched == 0 & ~strcmp(arglist{ii,1},''))
        error(['Argument ' arglist{ii,1} ' is not a valid argument for this function']);
   end
end
% match the remaining inputs/outputs to the unused arguments and test for missing required inputs
 for jj=1:size(narglist,1)
     if  strcmp(paramtype{jj}, 'Toggle')
        if (narglist{jj,2} ==0)
          narglist{jj,1} = ''; 
        end;
        narglist{jj,2} = ''; 
     end;
     if  ~strcmp(narglist{jj,2},'__input') && ~strcmp(narglist{jj,2},'__output') && istoggle(jj) && ~ was_set(jj)
          narglist{jj,1} = ''; 
          narglist{jj,2} = ''; 
     end;
     if strcmp(narglist{jj,2}, '__input')
      if (minval{jj} == 0)  % meaning this input is required
        if (nextinput > size(Inputs)) 
           error(['Required input ' narglist{jj,1} ' has no corresponding input in the list!']); 
        else
          narglist{jj,2} = 'OK_in';
          nextinput = nextinput + 1;
        end
      else  % this is an optional input
        if (nextinput <= length(Inputs)) 
          narglist{jj,2} = 'OK_in';
          nextinput = nextinput + 1;
        else 
          narglist{jj,1} = '';
          narglist{jj,2} = '';
        end;
      end;
     else 
     if strcmp(narglist{jj,2}, '__output')
      if (minval{jj} == 0) % this is a required output
        if (nextoutput > nargout & nargout > 1) 
           error(['Required output ' narglist{jj,1} ' is not stated in the assignment list!']); 
        else
          narglist{jj,2} = 'OK_out';
          nextoutput = nextoutput + 1;
          NumReqOutputs = NumReqOutputs-1;
        end
      else % this is an optional output
        if (nargout - nextoutput >= NumReqOutputs) 
          narglist{jj,2} = 'OK_out';
          nextoutput = nextoutput + 1;
        else 
          narglist{jj,1} = '';
          narglist{jj,2} = '';
        end;
      end
     end
  end
end
if nargout
   varargout = cell(1,nargout);
else
  varargout = cell(1,1);
end
global KhorosRoot
if exist('KhorosRoot') && ~isempty(KhorosRoot)
w=['"' KhorosRoot];
else
if ispc
  w='"C:\Program Files\dip\khorosBin\';
else
[s,w] = system('which cantata');
w=['"' w(1:end-8)];
end
end
[varargout{:}]=callKhoros([w 'kwt"  '],Inputs,narglist);
Back to Top