开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 图形图像 > 波变换 > 查看源码
demolift.m
资源名称：wavelet_toolbox.rar　[点击查看]
上传用户：haiyisale
上传日期：2013-01-09
资源大小：3246k
文件大小：45k
源码类别：

波变换

开发平台：

Matlab

demolift.m：源码内容
  1. function varargout = demolift(varargin)
  2. %DEMOLIFT Demonstrates Lifting functions in the Wavelet Toolbox.  
  3. %
  4. % This is a slideshow file for use with wshowdrv.m
  5. % To see it run, type 'wshowdrv demolift', 
  7. %   M. Misiti, Y. Misiti, G. Oppenheim, J.M. Poggi 10-Dec-2002.
  8. %   Last Revision: 04-Sep-2003.
  9. %   Copyright 1995-2004 The MathWorks, Inc.
  10. %   $Revision: 1.1.6.3 $  $Date: 2004/04/13 00:39:08 $
  12. % Initialization and Local functions if necessary.
  13. if nargin>0
  14. action = varargin{1};
  15.     switch action
  16.         case 'addHelp'
  17.             % Add Help Item.
  18.             %---------------
  19.             hdlFig = varargin{2};
  20.             wfighelp('addHelpItem',hdlFig,'Wavelet Decomposition','DW1D_DECOMPOS');
  21.             
  22.         case 'auto'    , wshowdrv('#autoMode',mfilename,'close');
  23.         case 'gr_auto' , wshowdrv('#gr_autoMode',mfilename,'close');
  24.         case 'getFigParam'
  25.             figName  = 'Lifting';
  26.             showType = 'mix10';
  27.             varargout = {figName,showType};
  28.             
  29.         case 'initShowViewer'
  30.             figHandle = varargin{2};
  31.             slideData = get(figHandle,'Userdata');
  32.             autoHndl = slideData.autoHndl;
  33.             propAuto = get(autoHndl,{'Units','Position','BackGroundColor'});
  34.             txtBkColor = get(slideData.slitxtHndl,'BackGroundColor');
  35.             pos  = propAuto{2};
  36.             hBtn = pos(4);
  37.             pos(2) = pos(2)-2*hBtn;
  38.             popstr = strvcat('db1','db2');
  39.             txtstr = 'Wavelet';
  40.             propUIC = {'Units',propAuto{1},'Position',pos,'BackGroundColor',txtBkColor};
  41.             slideData.txtLocHandle = uicontrol('Parent',figHandle,...
  42.                 'style','text','String',txtstr,propUIC{:});
  43.             pos(2) = pos(2)-hBtn+hBtn/2;
  44.             propUIC = {'Units',propAuto{1},'Position',pos,'BackGroundColor',propAuto{3}};
  45.             propAuto{2} = pos;
  46.             propUIC{4}  = propAuto{2};
  47.             slideData.popLocHandle = uicontrol('Parent',figHandle,...
  48.                 'style','popupmenu','String',popstr,propUIC{:}, ...
  49.                 'Tag','popWAV','Enable','Off');
  51.             pos(2) = pos(2)-hBtn;
  52.             popstr = strvcat('blocks','bumps','heavy sine','doppler','quadchirp','mishmash');
  53.             txtstr = 'Signal';
  54.             propUIC = {'Units',propAuto{1},'Position',pos,'BackGroundColor',txtBkColor};
  55.             slideData.txtLocHandle(2) = uicontrol('Parent',figHandle,...
  56.                 'style','text','String',txtstr,propUIC{:});
  57.             pos(2) = pos(2)-hBtn+hBtn/2;
  58.             propUIC = {'Units',propAuto{1},'Position',pos,'BackGroundColor',propAuto{3}};
  59.             propAuto{2} = pos;
  60.             propUIC{4}  = propAuto{2};
  61.             slideData.popLocHandle(2) = uicontrol('Parent',figHandle,...
  62.                 'style','popupmenu','String',popstr,propUIC{:}, ...
  63.                 'Tag','popSIG','Enable','Off');
  64.  
  65.             pos(2) = pos(2)-hBtn;
  66.             editstr = '';
  67.             txtstr = 'Nb. Kept Cfs.';
  68.             propUIC = {'Units',propAuto{1},'Position',pos,'BackGroundColor',txtBkColor};
  69.             slideData.txtLocHandle(3) = uicontrol('Parent',figHandle,...
  70.                 'style','text','String',txtstr,propUIC{:});
  71.             pos(2) = pos(2)-hBtn+hBtn/2;
  72.             propUIC = {'Units',propAuto{1},'Position',pos,'BackGroundColor',propAuto{3}};
  73.             pos(4) = pos(4)/2;
  74.             pos(2) = pos(2)+hBtn/2;
  75.             propUIC{4}  = pos;
  76.             edit_CB = ['demolift(''Local_FUN'',''edit_CB'');'];
  77.             slideData.editLocHandle = uicontrol('Parent',figHandle,...
  78.                 'style','edit','String',editstr,propUIC{:},...
  79.                 'Tag','ediCFS','Enable','Off','Callback',edit_CB);
  80.             
  81.             set(figHandle,'Userdata',slideData);
  82.  
  83.         case 'Local_FUN'
  84.             funNAME = varargin{2};
  85.             nbOUT = nargout;
  86.             varargout = cell(1,nbOUT);
  87.             switch nbOUT
  88.                 case 0 , feval(funNAME,varargin{3:end});
  89.                 case 1 , varargout{1} = feval(funNAME,varargin{3:end});
  90.                 case 2 , [varargout{1},varargout{2}] = feval(funNAME,varargin{3:end});    
  91.         end
  92. end
  93. return
  94. end
  96. if nargout<1,
  97.   wshowdrv(mfilename)
  98. else
  99.   idx = 0;
  100.   
  101.   %========== Slide 1 ==========
  102.   idx = idx+1;
  103.   slide(idx).code = {
  104. 'figHandle = gcf;',
  105.     'set(figHandle,''Position'',[240 100 752 612]);'
  106.     'demolift(''Local_FUN'',''set_UI_Graphics'',-1); ',
  107. '' };
  109.   slide(idx).text = {
  110. '',
  111. ' Press the "Start" button to see a demonstration of Lifting tools.',
  112. '',
  113. '',
  114. ' This demo uses Wavelet Toolbox functions.',
  115. ''};
  117.   %========== Slide 2 ==========
  118.   idx = idx+1;
  119.   slide(idx).code = {
  120. 'popWAV = findall(figHandle,''style'',''popupmenu'',''tag'',''popWAV'');',
  121.     'popSIG = findobj(figHandle,''style'',''popupmenu'',''Tag'',''popSIG'');',
  122. 'ediCFS = findobj(figHandle,''style'',''edit'',''Tag'',''ediCFS''); drawnow',
  123. 'set([popWAV,popSIG,ediCFS],''enable'',''Off''); ',
  124. 'set(popWAV,''Value'',1, ''String'',strvcat(''db1'',''db2'')); ',
  125. 'set(ediCFS,''String'',''''); ',
  126.     
  127.     'blanc = '' '';';  
  128.     'S1  = ''LIFTING (1)'';',
  129.     'S2  = ''-----------'';',
  130.     'S3  = blanc;',
  131.     'S4  = '' In the Wavelet Toolbox, the lifting functions are'';',
  132.     'S5  = '' designed to implement the discrete wavelet transform'';',
  133.     'S6  = '' and the inverse discrete wavelet transform using the'';',
  134.     'S7  = '' "LIFTING SCHEME" tool (See W. SWELDENS).'';',
  135.     'S8  = blanc;',
  136.     'S9  = '' The main feature of the Lifting Scheme Tool is that all'';',
  137.     'S10 = '' constructions are derived in the spatial domain. Staying in '';',
  138.     'S11 = '' spatial domain leads to two major advantages:'';',
  139.     'S12 = blanc;',    
  140.     'S13 = ''     1) the machinery of Fourier analysis is not required'';',
  141.     'S14 = ''     2) lifting leads to algorithms that can easily be  '';',
  142.     'S15 = ''        generalized to more complex geometric situations.'';',
  143.     'Infostr = strvcat(S1,S2,S3,S4,S5,S6,S7,S8,S9,S10,S11,S12,S13,S14,S15);',
  145.     'demolift(''Local_FUN'',''set_UI_Graphics'',0);',
  146.     'hdl_UI_Graphics = wtbxappdata(''get'',figHandle,''hdl_UI_Graphics'');',
  147.     'ax = hdl_UI_Graphics.ax ; t = hdl_UI_Graphics.t;',
  148.     'set(t(1),''String'',Infostr,''FontName'',''Courier New'',''FontSize'',10);',
  149. '' };
  151.   slide(idx).text = {
  152.     ' ',    
  153.     '   type: doc lifting',
  154.     ' ',    
  155.     '   or click on the "Info" button ',
  156.     ' ',
  157.     '   to see more information about lifting functions',
  158.   ''};
  160.   slide(idx).info = 'lifting';
  161.   
  162.   %========== Slide 3 ==========
  163.   idx = idx+1;
  164.   slide(idx).code = {
  165.     'S1  = ''LIFTING (2)'';',
  166.     'S2  = ''-----------'';',
  167.     'S3  = blanc;',
  168.     'S4  = '' To implement the discrete wavelet transform, the POLYPHASE'';',
  169.     'S5  = '' technique is used. So the Lifting Wavelet Transform (LWT) may be'';',
  170.     'S6  = '' expressed with ELEMENTARY LIFTING STEPS.'';',
  171.     'S7  = blanc;',
  172.     'S8  = '' There are two types of elementary lifting steps:'';',
  173.     'S9  = ''     "primal" lifting and "dual" lifting'';',
  174.     'S10 = blanc;',
  175.     'S11 = '' This implementation leads to three major advantages: '';',
  176.     'S12 = blanc;',
  177.     'S13 = ''     1) the traditional "well known" wavelets used in the DWT '';',
  178.     'S14 = ''        transform may have one or more "lifted" decompositions'';'
  179.     'S15 = ''     2) the inverse discrete wavelet transform is very easily'';',
  180.     'S16 = ''        computed'';',    
  181.     'S17 = ''     3) "generalized" wavelets can be easily constructed.'';',
  182.     
  183.     'Infostr = strvcat(S1,S2,S3,S4,S5,S6,S7,S8,S9,S10,S11,S12,S13,S14,S15,S16,S17);',
  184.     'set(t(1),''String'',Infostr,''FontName'',''Courier New'',''FontSize'',10);',
  185. '' };
  187.   slide(idx).text = {
  188.     ' ',    
  189.     '   type: doc lifting',
  190.     ' ',    
  191.     '   or click on the "Info" button ',
  192.     ' ',
  193.     '   to see more information about lifting functions',
  194.   ''};
  196.   slide(idx).info = 'lifting';
  197.   
  198.   %========== Slide 4 ==========
  199.   idx = idx+1;
  200.   slide(idx).code = {
  201.     'S1  = ''LIFTING (3)'';',
  202.     'S2  = ''-----------'';',
  203.     'S3  = blanc;',
  204.     'S4  = '' Each elementary lifting step is defined by:'';',
  205.     'S5  = ''     - a TYPE : "primal" (p) or "dual" (d)'';',
  206.     'S6  = ''     - a LAURENT POLYNOMIAL P(z) which defines the coefficients'';',
  207.     'S7  = ''       of the filter associated to the lifting step'';',
  208.     'S8  = blanc;',
  209.     'S9  = '' A LIFTING SCHEME (LS) is a set of elementary lifting steps'';',
  210.     'S10 = '' completed by an information about normalization coefficients.'';',
  211.     'S11 = '' (see W. SWELDENS papers for more information about theory)'';',
  212.     'S12  = blanc;',
  213.     'S13 = '' Example 1: lifting decomposition of DWT for the db1 wavelet '';',
  214.     'bcar = '' '';',
  215.     'LS = liftwave(''db1'');',
  216.     'TMP = displs(LS,''%8.4f''); nbROW = size(TMP,1);',
  217.     'S14 = [bcar(ones(nbROW,12)),TMP];',
  218.     'S15  = blanc;',
  219.     'S16 = '' See next slide for comments ...'';',    
  220.     'Infostr = strvcat(S1,S2,S3,S4,S5,S6,S7,S8,S9,S10,S11,S12,S13,blanc,S14,S15,S16);',
  221.     
  222.     'set(t(1),''String'',Infostr,''FontName'',''Courier New'',''FontSize'',10);',
  223. '' };
  225.   slide(idx).text = {
  226.     ' ',    
  227.     '   type: doc lifting',
  228.     ' ',    
  229.     '   or click on the "Info" button ',
  230.     ' ',
  231.     '   to see more information about lifting functions',
  232.   ''};
  234.   slide(idx).info = 'lifting';
  235.  
  236.   %========== Slide 5 ==========
  237.   idx = idx+1;
  238.   slide(idx).code = {
  239.     'S1  = ''LIFTING (3) - continu '';',
  240.     'S2  = ''----------------------'';',
  241.     'S3  = blanc;',
  242.     'S4  = '' Example 2: lifting decomposition of DWT for the wavelet db2'';',
  243.     'bcar = '' '';',
  244.     'LSdb2 = liftwave(''db2'');',
  245.     'PStr = disp(lp(LSdb2{2,2},LSdb2{2,3}));',
  246.     'PStr = demolift(''Local_FUN'',''wlift_str_util'',''replace'',PStr);',
  247.     'TMP = displs(LSdb2,''%8.4f''); nbROW = size(TMP,1);',
  248.     'S5  = [bcar(ones(nbROW,12)),TMP];',
  249.     'S6  = '' The last line contains the normalization coefficients.'';', 
  250.     'S7  = '' In the other lines:'';',
  251.     'S8  = ''    - The third column contains the Laurent Polynomial'';', 
  252.     'S9  = ''      maximal degree.'';', 
  253.     'S10 = ''    - The second column contains the Laurent Polynomial'';',
  254.     'S11 = ''      coefficients. For example in the second line: '';',
  255.     'S12 = [''         P'', PStr ''  ''];',
  256.     'Infostr = strvcat(S1,S2,S3,S4,blanc,S5,S6,S7,blanc,S8,S9,blanc,S10,S11,S12);',
  257.     
  258.     'set(t(1),''String'',Infostr,''FontName'',''Courier New'',''FontSize'',10);',
  259. '' };
  261.   slide(idx).text = {
  262.     ' ',    
  263.     '   type: doc lifting',
  264.     ' ',    
  265.     '   or click on the "Info" button ',
  266.     ' ',
  267.     '   to see more information about lifting functions',
  268.   ''};
  270.   slide(idx).info = 'lifting';
  271.  
  272.   %========== Slide 6 ==========
  273.   idx = idx+1;
  274.   slide(idx).code = {
  275.     'S1  = ''LIFTING (4)'';',
  276.     'S2  = ''------------'';',
  277.     'S4  = '' To implement LIFTING SCHEMES and POLYPHASE techniques'';',
  278.     'S5  = '' the Wavelet Toobox needs two types of object:'';',
  279.     'S7  = ''    - LP Object to manage Laurent Polynomials.'';',
  280.     'S8  = ''    - LM Object to manage Laurent Matrices.'';',
  281.     'S10 = '' The last one is used in particular for the Factorization'';',
  282.     'S11 = '' Algorithm of "usual" wavelets.'';',
  283.     'Infostr = strvcat(S1,S2,blanc,S4,S5,blanc,S7,blanc,S8,blanc,S10,S11);',
  284.     
  285.     'set(t(1),''String'',Infostr,''FontName'',''Courier New'',''FontSize'',10);',
  286. '' };
  288.   slide(idx).text = {
  289.     ' ',    
  290.     '   type: doc lifting',
  291.     ' ',    
  292.     '   or click on the "Info" button ',
  293.     ' ',
  294.     '   to see more information about lifting functions',
  295.   ''};
  297.   slide(idx).info = 'lifting';
  298.   
  299.   %========== Slide 7 ==========
  300.   idx = idx+1;
  301.   slide(idx).code = {
  302.     'S1  = ''LIFTING (5): lp OBJECT'';',
  303.     'S2  = ''-------------------------'';',
  304.     'HelpSTR = help(''lp'');',
  305.     'HelpSTR = demolift(''Local_FUN'',''wlift_str_util'',''replacePOW'',HelpSTR);',
  306.     'Infostr = strvcat(S1,S2,HelpSTR);',
  307.     'set(t(1),''String'',Infostr,''FontName'',''Courier New'',''FontSize'',9);',
  308. '' };
  310.   slide(idx).text = {
  311.     ' ',    
  312.     '   type: doc lp',
  313.     ' ',    
  314.     '   or click on the "Info" button ',
  315.     ' ',
  316.     '   to see more information about lp OBJECT',
  317.   ''};
  319.   slide(idx).info = 'lp'; 
  320.   
  321.   %========== Slide 8 ==========
  322.   idx = idx+1;
  323.   slide(idx).code = {
  324.     'S1  = ''LIFTING (6): lm OBJECT'';',
  325.     'S2  = ''-------------------------'';',
  326.     'HelpSTR = help(''lm'');',
  327.     'HelpSTR = demolift(''Local_FUN'',''wlift_str_util'',''replacePOW'',HelpSTR);',
  328.     'Infostr = strvcat(S1,S2,HelpSTR);',
  329.     'set(t(1),''String'',Infostr,''FontName'',''Courier New'',''FontSize'',10);',
  330. '' };
  332.   slide(idx).text = {
  333.     ' ',    
  334.     '   type: doc lm',
  335.     '   ',    
  336.     '   or click on the "Info" button ',
  337.     '   ',
  338.     '   to see more information about lm OBJECT',
  339.   ''};
  341.   slide(idx).info = 'lm';
  342.   
  343.   %========== Slide 9 ==========
  344.   idx = idx+1;
  345.   slide(idx).code = {
  346.   'set(popWAV,''enable'',''off'');',
  347.       
  348.     'S1  = ''LIFTING (7): MAIN LIFTING FUNCTIONS'';',
  349.     'S2  = ''-----------------------------------'';',
  350.     'S3  = blanc;',
  351.     
  352.     'S4  = ''addlift   - Adding lifting and dual-lifting steps.'';', 
  353.     'S5  = ''bswfun    - Biorthogonal scaling and wavelet functions.'';', 
  354.     'S6  = ''displs    - Display lifting scheme.'';',
  355.     'S7  = ''filt2ls   - Filters to lifting scheme.'';',
  356.     'S8  = ''ilwt      - Wavelet reconstruction 1-D using lifting.'';',
  357.     'S9  = ''ilwt2     - Wavelet reconstruction 2-D using lifting.'';',
  358.     'S10 = ''liftfilt  - Apply elementary lifting steps on filters.'';', 
  359.     'S11 = ''liftwave  - Lifting structure for wavelets.'';',
  360.     'S12 = ''lsinfo    - Information about lifting schemes.'';',
  361.     'S13 = ''lwt       - Wavelet decomposition 1-D using lifting.'';',
  362.     'S14 = ''lwt2      - Wavelet decomposition 2-D using lifting.'';',
  363.     'S15 = ''lwtcoef   - Extract or reconstruct 1-D wavelet coefficients.'';',
  364.     'S16 = ''lwtcoef2  - Extract or reconstruct 2-D wavelet coefficients.'';',
  365.     'S17 = ''wave2lp   - Laurent polynomial associated to a wavelet.'';',
  366.     'S18 = ''wavenames - Wavelet names information.'';',
  367.     'Infostr = strvcat(S1,S2,S3,S4,S5,S6,S7,S8,S9,S10,S11,S12,S13,S14,S15,S16,S17,S18);',
  368.     
  369.     'set(t(1),''String'',Infostr,''FontName'',''Courier New'',''FontSize'',10);',
  370. '' };
  371.     
  372.   slide(idx).text = {
  373.     ' ',    
  374.     '   type: doc lifting',
  375.     ' ',    
  376.     '   or click on the "Info" button ',
  377.     ' ',
  378.     '   to see more information about lifting functions',
  379.   ''};
  381.   slide(idx).info = 'lifting';
  382.   
  383.   %========== Slide 10 ==========
  384.   idx = idx+1;
  385.   slide(idx).code = {
  386. 'set(popWAV,''enable'',''On'');',
  387.     
  388.     'S1  = ''LIFTING (8): FILTERS and POLYNOMIALS ...'';',
  389.     'S2  = ''----------------------------------------'';',
  390.     'S3  = blanc;',
  391.     'S4  = '' Choose the Wavelet (db1 or db2) with the popupMenu'';',
  392.     'S5  = '' on the right part of the window.'';',
  393.     
  394.     'Infostr = strvcat(S1,S2,S3,S4,S5);',
  395.     'set(t(1),''String'',Infostr);',
  396.     'demolift(''Local_FUN'',''set_UI_Graphics'',0);', 
  397. '' };
  399.   slide(idx).text = {
  400.     ' ',    
  401.     '   type: doc lifting',
  402.     ' ',    
  403.     '   or click on the "Info" button ',
  404.     ' ',
  405.     '   to see more information about lifting functions',
  406.   ''};
  408.   slide(idx).info = 'lifting';
  410.   %========== Slide 11 ==========
  411.   idx = idx+1;
  412.   slide(idx).code = {
  413.   'set(popWAV,''enable'',''off'');',
  414. 'lstWAV = get(popWAV,''String''); idxWAV = get(popWAV,''Value'');',
  415. 'wname = deblankl(lstWAV(idxWAV,:));',
  416. 'comma = char(39);',
  417. 'newTXT = [''        wname = '',comma,wname,comma,'';''];';
  418. 'wshowdrv(''#modify_Comment'',figHandle,3,newTXT);',
  419.     'W = load(wname);',
  420.     'F = sqrt(2)*W.(wname);',
  421.     'FStr =  [ ''F = ['' num2str(F,4) '']'' ];',
  422.     
  423.     'S1  = ''LIFTING (9): FILTERS and POLYNOMIALS ...'';',
  424.     'S2  = ''----------------------------------------'';',
  425.     'S3  = blanc;',
  426.     'S4  = FStr;',
  427.     'Infostr = strvcat(S1,S2,S3,S4);', 
  428.     
  429.     'set(t(1),''String'',Infostr);',
  430.     'demolift(''Local_FUN'',''set_UI_Graphics'',0);', 
  431. '' };
  433.   slide(idx).text = {
  434. ' % Load wavelet filter.',
  435. '',
  436. '        wname = ''db1'';',
  437. '        load(wname);',
  438. '',
  439. ' % Compute the associated normalized filter F.',
  440.     '',
  441.     '        F = sqrt(2)*wname;',
  442. ''};
  444.   slide(idx).info = 'lifting';
  446.   %========== Slide 12 ==========
  447.   idx = idx+1;
  448.   slide(idx).code = {
  449.     'H = lp(F);',
  450.     'HStr = disp(H);',
  451.     'HStr = demolift(''Local_FUN'',''wlift_str_util'',''replace'',HStr);',
  453.     'S1  = ''LIFTING (10): FILTERS and POLYNOMIALS ...'';',
  454.     'S2  = ''-----------------------------------------'';',    
  455.     'S3  = blanc;',
  456.     'S4  = FStr;',
  457.     'S5  = blanc;',
  458.     'S6  = HStr;',
  459.     'Infostr = strvcat(S1,S2,S3,S4,S5,S6);', 
  460.     
  461. 'set(t(1),''String'',Infostr);',
  462. '' };
  464.   slide(idx).text = {
  465. ' % Compute the associated reconstruction low-pass',
  466.     ' % Laurent polynomial H.',
  467. '',
  468.   '        H = lp(F);',
  469. ''};
  471.   slide(idx).info = 'lp';
  472.   
  473. %   %========== Slide 13 ==========
  474.   idx = idx+1;
  475.   slide(idx).code = {
  476.     'G = lp(1,-1)*modulate(reflect(H));',
  477.     'GStr = disp(G);',
  478.     'GStr = demolift(''Local_FUN'',''wlift_str_util'',''replace'',GStr);',
  479.     
  480.     'S1  = ''LIFTING (11): FILTERS and POLYNOMIALS ...'';',
  481.     'S2  = ''-----------------------------------------'';',
  482.     'S3  = blanc;',
  483.     'S4  = FStr;',
  484.     'S5  = HStr;',
  485.     'S6  = blanc;',
  486.     'S7  = GStr;',
  487.     'Infostr = strvcat(S1,S2,S3,S4,S5,S6,S7);', 
  488.     
  489. 'set(t(1),''String'',Infostr);',
  490. '' };
  492.   slide(idx).text = {
  493. ' % Compute the associated reconstruction high-pass',
  494.     ' % Laurent polynomial G.',
  495. '',
  496.   '        G = lp(1,-1)*modulate(reflect(H));',
  497. ''};
  499.   slide(idx).info = 'lp/modulate';
  500.   
  501.   %========== Slide 14 ==========
  502.   idx = idx+1;
  503.   slide(idx).code = {
  504.     'E_H = even(H); H_EvenPart = disp(E_H);',
  505.     'O_H = odd(H);  H_OddPart  = disp(O_H);',
  506.     'H_EvenPart = demolift(''Local_FUN'',''wlift_str_util'',''replace'',H_EvenPart);',
  507.     'H_OddPart  = demolift(''Local_FUN'',''wlift_str_util'',''replace'',H_OddPart);',
  508.   
  509.     'S1  = ''LIFTING (12): FILTERS and POLYNOMIALS ...'';',
  510.     'S2  = ''-----------------------------------------'';',
  511.     'S3  = blanc;',
  512.     'S4  = FStr;',
  513.     'S5  = HStr;',
  514.     'S6  = GStr;',
  515.     'S7  = blanc;',
  516.     'S8  = H_EvenPart;',
  517.     'S9  = H_OddPart;',
  518.     'Infostr = strvcat(S1,S2,S3,S4,S5,S6,S7,S8,S9);', 
  519.     
  520. 'set(t(1),''String'',Infostr);',
  521. '' };
  523.   slide(idx).text = {
  524. ' % Compute the even and odd parts of H ...',
  525. '',
  526.   '        E_H = even(H);',
  527.     '        O_H = odd(H); ',
  528. ''};
  530.   slide(idx).info = 'lp/even';
  532.   %========== Slide 15 ==========
  533.   idx = idx+1;
  534.   slide(idx).code = {
  535.     'E_G = even(G); G_EvenPart = disp(E_G);',
  536.     'O_G = odd(G);  G_OddPart  = disp(O_G);',
  537.     'G_EvenPart = demolift(''Local_FUN'',''wlift_str_util'',''replace'',G_EvenPart);',
  538.     'G_OddPart  = demolift(''Local_FUN'',''wlift_str_util'',''replace'',G_OddPart);',
  539.     
  540.     'S1  = ''LIFTING (13): FILTERS and POLYNOMIALS ...'';',
  541.     'S2  = ''-----------------------------------------'';',
  542.     'S3  = blanc;',
  543.     'S4  = FStr;',
  544.     'S5  = HStr;',
  545.     'S6  = GStr;',
  546.     'S7  = blanc;',
  547.     'tmpSTR_E = strvcat(H_EvenPart,G_EvenPart);',
  548.     'tmpSTR_O = strvcat(H_OddPart,G_OddPart);',
  549.     'tmpSTR_EO = [tmpSTR_E , bcar(ones(2,8)) , tmpSTR_O];';
  550.     'Infostr = strvcat(S1,S2,S3,S4,S5,S6,S7,tmpSTR_EO);', 
  551.     
  552. 'set(t(1),''String'',Infostr);',
  553. '' };
  555.   slide(idx).text = {
  556. ' % Compute the even and odd parts of G ...',
  557. '',
  558.   '        E_G = even(G);',
  559.     '        O_G = odd(G); ',
  560. ''};
  562.   slide(idx).info = 'lp/odd';
  564.   %========== Slide 16 ==========
  565.   idx = idx+1;
  566.   slide(idx).code = {
  567.     'S1  = ''LIFTING (14): FILTERS and POLYNOMIALS ...'';',
  568.     'S2  = ''-----------------------------------------'';',
  569.     'S3  = blanc;',
  570.     'S4  = FStr;',
  571.     'S5  = HStr;',
  572.     'S6  = GStr;',
  573.     'S7  = blanc;',
  574.     'Infostr = strvcat(S1,S2,S3,S4,S5,S6,S7,tmpSTR_EO);', 
  575.     
  576. 'set(t(1),''String'',Infostr);',
  577. '' };
  579.   slide(idx).text = {
  580. ' % Now we have the even and odd parts of H and G.',
  581. ' % So we can built the Polyphase Matrix.',
  582.   '',
  583.   '        |  E_H    E_G  | ',
  584.     '        |  O_H    O_G  | ',
  585. ''};
  587.   slide(idx).info = 'lp/ppm';
  589.   %========== Slide 17 ==========
  590.   idx = idx+1;
  591.   slide(idx).code = {
  592.     'PolyMAT = ppm(H,G);',
  593.     'M_dispSTR = disp(PolyMAT);',
  594.     
  595.     'S1  = ''LIFTING (15): POLYPHASE MATRIX ...'';',
  596.     'S2  = ''----------------------------------'';',
  597.     'S3  = blanc;',
  598.     'S4  = M_dispSTR;',
  599.     'if isequal(wname,''db2''),',
  600.     '   S4  = demolift(''Local_FUN'',''wlift_str_util'',''replaceMAT'',S4);',
  601.     'end',
  602.     
  603.     'Infostr = strvcat(S1,S2,S3,S4);', 
  604.     
  605. 'set(t(1),''String'',Infostr);',
  606. '' };
  608.   slide(idx).text = {
  609. ' % Compute the associated polyphase matrix: PolyMAT',
  610. '',
  611.   '        PolyMAT = ppm(H,G);',
  612. ''};
  614.   slide(idx).info = 'lp/ppm';
  616.   %========== Slide 18 ==========
  617.   idx = idx+1;
  618.   slide(idx).code = {
  619.     'MatFACT = ppmfact(H,G);',
  620.     'DEC = MatFACT{1}; tmpSTR = [];',
  621.     'for j = 1:length(DEC), ',
  622.     '    tmp{j} = disp(DEC{j},''M''); tmp{j}(:,[1:4]) = '' '';',
  623.     '    if isequal(wname,''db2'')',
  624.     '        tmp{j}(2:3,:) = [];',
  625.     '        tmp{j} = demolift(''Local_FUN'',''wlift_str_util'',''replaceMAT_2'',tmp{j});',    
  626.     '    end',
  627.     'end',
  628.     'if isequal(wname,''db2'')',
  629.     '   tmpSTR = strvcat([tmp{1},tmp{2}],'' '','' '',[tmp{3},tmp{4}]);',
  630.     '   interSTR = ''none'';',    
  631.     'else',
  632.     '   for j = 1:length(DEC), tmpSTR = [tmpSTR , tmp{j}]; end',
  633.     '   interSTR = ''tex'';',
  634.     'end',
  635.     
  636.     'S1  = ''LIFTING (16): POLYPHASE MATRIX DECOMPOSITION'';',
  637.     'S2  = ''--------------------------------------------'';',
  638.     'S3  = blanc;',
  639.     'S4  = '' The polyphase matrix may be decomposed in '';',
  640.     'S5  = '' elementary lifting steps matrices:'';',
  641.     'S6  = blanc;',
  642.     'S7  = tmpSTR;',
  643.     'S8  = blanc;',
  644.     'S9  = '' The last matrix is devoted to a normalization operation.'';',
  645.     
  646.     'Infostr = strvcat(S1,S2,S3,S4,S5,S6,S7,S8,S9);', 
  647.     
  648. 'set(t(1),''Interpreter'',interSTR,''String'',Infostr);',
  649. '' };
  651.   slide(idx).text = {
  652. ' % Compute Matrix decomposition(s) ...',
  653. '',
  654.   '        MatFACT = ppmfact(H,G);',
  655. ''};
  657.   slide(idx).info = 'lp/ppmfact';
  659.   %========== Slide 19 ==========
  660.   idx = idx+1;
  661.   slide(idx).code = {
  662.     '[LS_d,LS_p] = pmf2ls(MatFACT,''t''); LSd = LS_d{1}; LSp = LS_p{1};',
  663.     'LS_d_STR  = displs(LSd,''%8.4f'');',
  664.     'bcar = '' '';',
  665.     'nbROW = size(LS_d_STR,1);',
  666.     'LS_d_STR = [bcar(ones(nbROW,4)),LS_d_STR];',
  667.     'LS_p_STR  = displs(LSp,''%8.4f''); ',
  668.     'bcar = '' '';',
  669.     'nbROW = size(LS_p_STR,1);',
  670.     'LS_p_STR = [bcar(ones(nbROW,4)),LS_p_STR];',
  671.     
  672.     'S1  = ''LIFTING (17): POLYPHASE MATRIX DECOMPOSITION'';',
  673.     'S2  = ''--------------------------------------------'';',
  674.     'S3  = blanc;',
  675.     'S4  = ''Now we can build the Lifting Scheme(s) associated'';',
  676.     'S5  = ''with the previous decomposition:'';',
  677.     'S6  = blanc;',
  678.     'S7  = LS_d_STR;',
  679.     'S8  = ''  or '';',
  680.     'S9  = blanc;',
  681.     'S10  = LS_p_STR;',
  682.     'Infostr = strvcat(S1,S2,S3,S4,S5,S6,S7,S8,S9,S10);', 
  683.     
  684. 'set(t(1),''Interpreter'',''tex'',''FontSize'',10,''String'',Infostr);',
  685. '' };
  687.   slide(idx).text = {
  688. ' % Computing and Viewing lifting schemes table ...',
  689. '',
  690.   '        [LS_d,LS_p] = pmf2ls(MatFACT,''t'');',
  691. ''};
  693.   slide(idx).info = 'pmf2ls';
  695.   %========== Slide 20 ==========
  697.   idx = idx+1;
  698.   slide(idx).code = {
  699.     'LS = liftwave(wname);',
  700.     'LS_STR  = displs(LS,''%8.4f'');',
  701.     'bcar = '' '';',
  702.     'nbROW = size(LS_STR,1);',
  703.     'LS_STR = [bcar(ones(nbROW,4)),LS_STR];',
  705.     'S1  = ''LIFTING (18): LIFTING SCHEME'';',
  706.     'S2  = ''----------------------------'';',
  707.     'S3  = blanc;',
  708.     'S4  = ''For "classical" wavelets we can obtain directly '';',
  709.     'S5  = ''the associated Lifting Scheme using the LIFTWAVE function:'';',
  710.     'S6  = blanc;',
  711.     'S7  = ''   LS = liftwave(wname);'';',
  712.     'S8  = blanc;',    
  713.     'S9  = ''So we get:'';',
  714.     'S10  = blanc;',
  715.     'S11 = LS_STR;',
  716.     
  717.     'Infostr = strvcat(S1,S2,S3,S4,S5,S6,S7,S8,S9,S10,S11);',
  718.     
  719. 'set(t(1),''FontSize'',10,''String'',Infostr);',
  720.     'demolift(''Local_FUN'',''set_UI_Graphics'',0);',
  721.     
  722. '' };
  724.   slide(idx).text = {
  725. ' % Compute the lifting scheme  ...',
  726. '',
  727.   '        ',
  728. ''};
  730.   slide(idx).info = 'liftwave';
  732.   %========== Slide 21 ==========
  733.   idx = idx+1;
  734.   slide(idx).code = {
  735.     'popstr = strvcat(''db1'',''db2'');',
  736.     'set(popWAV,''String'',popstr,''Value'',idxWAV)',      
  737.     'LS = liftwave(wname);',
  738.     
  739.     'demolift(''Local_FUN'',''set_UI_Graphics'',1);',
  740.     'demolift(''Local_FUN'',''plot_phi_psi'',LS,ax);', 
  741. '' };
  743.   slide(idx).text = {
  744. ' % Viewing Scaling function and Wavelet ... ',
  745. '',
  746.   '        ',
  747. ''};
  749.   slide(idx).info = 'liftwave';
  750.   
  751.   %========== Slide 22 ==========
  752.   idx = idx+1;
  753.   slide(idx).code = {
  754.     'liftSTEP = {''p'' , [-1 2 -1]/4 , [1] };',
  755.     'LSnew_1 = addlift(LS,liftSTEP);',
  756.     'LS_STR = displs(LSnew_1,''%8.4f'');',
  757.     'bcar = '' '';',
  758.     'nbROW = size(LS_STR,1);',
  759.     'LS_STR = [bcar(ones(nbROW,4)),LS_STR];',
  760.     
  761.     'newWAV = [wname ''_LS1''];',
  762.     'popstr = strvcat(''db1'',''db2'',newWAV);',
  763.     'set(popWAV,''String'',popstr,''Value'',3)',
  764.     
  765.     'S1  = ''LIFTING (19): LIFTING SCHEME'';',
  766.     'S2  = ''----------------------------'';',
  767.     'S3  = blanc;',
  768.     'S4  = ''Compute a new Lifting Scheme adding a new '';',
  769.     'S5  = ''elementary lifting step:'';',
  770.     'S6  = blanc;',
  771.     'S7  = LS_STR(end-2,:);',
  772.     'S8  = blanc;',    
  773.     'S9  = ''So we get:'';',
  774.     'S10 = blanc;',
  775.     'S11 = LS_STR;',
  776.     'S12 = blanc;',
  777.     'S13 = ''Plot the new "Scaling function" and "Wavelet".'';', 
  778.     'Infostr = strvcat(S1,S2,S3,S4,S5,S6,S7,S8,S9,S10,S11,S12,S13);',
  779.     
  780. 'set(t(1),''FontSize'',10,''String'',Infostr);',
  781.     'demolift(''Local_FUN'',''set_UI_Graphics'',0);', 
  782. '' };
  784.   slide(idx).text = {
  785. ' % Viewing a new Lifting Scheme ... ',
  786.   ' ',
  787. ' '};
  789.   slide(idx).info = 'addlift';
  790.   
  791.   %========== Slide 23 ==========
  792.   idx = idx+1;
  793.   slide(idx).code = {
  794.     'demolift(''Local_FUN'',''set_UI_Graphics'',1);', 
  795.     'demolift(''Local_FUN'',''plot_phi_psi'',LSnew_1,ax);', 
  796. '' };
  798.   slide(idx).text = {
  799. ' % Viewing new Scaling function and new Wavelet ... ',
  800.   ' ',
  801. ' '};
  803.   slide(idx).info = 'addlift';
  805.   %========== Slide 24 ==========
  806.   idx = idx+1;
  807.   slide(idx).code = {
  808.     'liftSTEP = { ''p'' , [3 2 1 7 -1 2 -1]/4  ,  [1] };',   
  809.     'LSnew_2 = addlift(LS,liftSTEP);',
  810.     'LS_STR = displs(LSnew_2,''%6.2f'');',
  811.     'bcar = '' '';',
  812.     'nbROW = size(LS_STR,1);',
  813.     'LS_STR = [bcar(ones(nbROW,4)),LS_STR];',
  814.     
  815.     'Lift_1 = [wname ''_LS1'']; Lift_2 = [wname ''_LS2''];',
  816.     'popstr = strvcat(''db1'',''db2'',Lift_1,Lift_2);',
  817.     'set(popWAV,''String'',popstr,''Value'',4)',    
  818.     
  819.     'S1  = ''LIFTING (20): LIFTING SCHEME'';',
  820.     'S2  = ''----------------------------'';',
  821.     'S3  = blanc;',
  822.     'S4  = ''Compute a new Lifting Scheme adding a new '';',
  823.     'S5  = ''elementary lifting step:'';',
  824.     'S6  = blanc;',
  825.     'S7  = LS_STR(end-2,:);',
  826.     'S8  = blanc;',    
  827.     'S9  = ''So we get:'';',
  828.     'S10 = blanc;',
  829.     'S11 = LS_STR;',
  830.     'S12 = blanc;',
  831.     'S13 = ''Plot the new "Scaling function" and "Wavelet".'';', 
  832.     'Infostr = strvcat(S1,S2,S3,S4,S5,S6,S7,S8,S9,S10,S11,S12,S13);',
  833.     
  834. 'set(t(1),''FontSize'',10,''String'',Infostr);',
  835.     'demolift(''Local_FUN'',''set_UI_Graphics'',0);', 
  836. '' };
  838.   slide(idx).text = {
  839. ' % Viewing a new Lifting Scheme ... ',
  840.   ' ',
  841. ' '};
  843.   slide(idx).info = 'addlift';
  844.   
  845.   %========== Slide 25 ==========
  846.   idx = idx+1;
  847.   slide(idx).code = {
  848.     'set([popWAV,popSIG],''Enable'',''Off'')',  
  849.     
  850.     'demolift(''Local_FUN'',''set_UI_Graphics'',1);', 
  851.     'demolift(''Local_FUN'',''plot_phi_psi'',LSnew_2,ax);', 
  852. '' };
  854.   slide(idx).text = {
  855. ' % Viewing new "Scaling function" and new "Wavelet" ... ',
  856.     ' ',    
  857.   '   This is not a wavelet in the usal sense since:'
  858.     ' ',    
  859.     '       mean(PSI) is not zero ! ',
  860.     ' ',
  861. ' '};
  863.   slide(idx).info = 'addlift';
  864.     
  865.   %========== Slide 26 ==========
  866.   idx = idx+1;
  867.   slide(idx).code = {
  868.     'set([popWAV,popSIG],''Enable'',''On'')',
  869.     
  870.     'S1  = ''LIFTING (21): ANALYSIS'';',
  871.     'S2  = ''----------------------'';',
  872.     'S3  = blanc;',
  873.     'S4  = ''Choose the SIGNAL to analyze and the WAVELET using '';',
  874.     'S5  = ''the popupmenus in the right part of the window.'';',
  875.     'S6  = blanc;',
  876.     'S7  = ''The LEVEL of decomposition is 5.'';',    
  877.     'S8  = blanc;',
  878.     'S9  = ''Then Press the "Next" button to perform the analysis.'';',
  879.     'Infostr = strvcat(S1,S2,S3,S4,S5,S6,S7,S8,S9);',
  880.     
  881. 'set(t(1),''FontSize'',10,''String'',Infostr);',
  882.     'demolift(''Local_FUN'',''set_UI_Graphics'',0);', 
  883. '' };
  885.   slide(idx).text = {
  886. ' % Analysis with a Lifting Scheme ... ',
  887.   ' ',
  888. ' '};
  890.   slide(idx).info = 'addlift';
  891.   
  892.   %========== Slide 27 ==========
  893.   idx = idx+1;
  894.   slide(idx).code = {
  895.     'set([popWAV,popSIG],''Enable'',''Off'')',
  896. 'lstSIG = get(popSIG,''String''); idxSIG = get(popSIG,''Value'');',
  897. 'sigNAM = deblankl(lstSIG(idxSIG,:));',
  898. 'idxWAV = get(popWAV,''Value'');',
  899.     'switch idxWAV ',
  900.     '    case {1,2} , LSnew = LS;',
  901.     '    case 3     , LSnew = LSnew_1;',
  902.     '    case 4     , LSnew = LSnew_2;',
  903.     'end', 
  904.     'x  = wnoise(sigNAM,10); level = 5;',
  905. 'lx = length(x);',
  906. 'xDEC = lwt(x,LSnew,level);',
  907. 'a = zeros(level,lx);',
  908. 'd = zeros(level,lx);',
  909.     'for k = 1:level',
  910.        'a(k,:) = lwtcoef(''a'',xDEC,LSnew,level,k);',
  911.        'd(k,:) = lwtcoef(''d'',xDEC,LSnew,level,k);' ,
  912.     'end',
  913.     'demolift(''Local_FUN'',''set_UI_Graphics'',2);',
  914.     'demolift(''Local_FUN'',''plot_ANALYSIS'',x,a,d,1,ax(4:6));',
  915.     'demolift(''Local_FUN'',''store_ANALYSIS'',x,xDEC,LSnew,level,ax(4:7),lx,''1D'');',
  916.     
  917. '' };
  919.   slide(idx).text = {
  920. ' % Compute the decomposition at level 5.',
  921.   ' ',
  922. '       xDEC = lwt(x,LSnew,level);',
  923. '       a = zeros(level,lx);',
  924. '       d = zeros(level,lx);',
  925.     '       for k = 1:level',
  926.     '           a(k,:) = lwtcoef(''a'',xDEC,LSnew,level,k);',
  927.     '           d(k,:) = lwtcoef(''d'',xDEC,LSnew,level,k);' ,
  928.     '       end',
  929.    ' ',   
  930. ' % And view the decomposition at level 1.',
  931.   ' ',
  932. ' '};
  934.   slide(idx).info = 'lwtcoef';
  935.   
  936.   %========== Slide 28 ==========
  937.   idx = idx+1;
  938.   slide(idx).code = {
  939.     'errorSTR = [];',
  940.     'for k = 1:level,',
  941.       'er = x-a(k,:);',
  942.       'for j=1:k , er = er-d(j,:); end',
  943.       'err(k) = max(abs(er));',
  944.       'errLEV = [''    level:'' int2str(k) '' - error = '' num2str(err(k))];',
  945.       'errorSTR = strvcat(errorSTR,errLEV);',
  946.     'end',
  947.     
  948. 'wshowdrv(''#modify_Comment'',figHandle,9,errorSTR);',
  949. 'set(popSIG,''enable'',''off'');',      
  950.     'demolift(''Local_FUN'',''set_UI_Graphics'',3);',
  951.     'demolift(''Local_FUN'',''plot_ANALYSIS'',x,a,d,2,ax(4:7));', 
  952. '' };
  954.   slide(idx).text = {
  955. ' % View the decomposition at level 2.',
  956.   ' ',
  957.   ' % S = A1 + D1 = A2 + D2 + D1 = A3 + D3 + D2 + D1 ...',
  958.   ' ',
  959.   ' % err_Level_1 = max( | S - A1 - D1 | )',
  960.   ' % err_Level_2 = max( | S - A1 - D1 - D2 | )',
  961.   ' % ...',
  962.   ' ',
  963.   ' ',    
  964. ' '};
  966.   slide(idx).info = 'lwtcoef';
  967.   
  968.   %========== Slide 29 ==========
  969.   idx = idx+1;
  970.   slide(idx).code = {
  971.     'set(ediCFS,''String'','' '',''Enable'',''Off'');',      
  972. 'wshowdrv(''#modify_Comment'',figHandle,3,errorSTR);',      
  973.     'demolift(''Local_FUN'',''set_UI_Graphics'',4);', 
  974.     'demolift(''Local_FUN'',''plot_ANALYSIS'',x,a,d,3,ax(4:8));',        
  975. '' };
  977.   slide(idx).text = {
  978. ' % View the decomposition at level 3.',
  979.   ' ',
  980.   ' ',    
  981. ' '};
  983.   slide(idx).info = 'addlift';
  984.   
  985.   %========== Slide 30 ==========
  986.   idx = idx+1;
  987.   slide(idx).code = {
  988. 'keptCFSMAX = length(x);',
  989. 'set(ediCFS,''String'',int2str(keptCFSMAX),''Enable'',''On'');',
  990.     'demolift(''Local_FUN'',''set_UI_Graphics'',5);', 
  991.     'demolift(''Local_FUN'',''plot_DECOMP'',x,xDEC,LSnew,level,ax(4:7),keptCFSMAX,''1D'');',
  992. '' };
  994.   slide(idx).text = {
  995.   ' ',      
  996. '   Choose the number of kept coefficients in the edit.',
  997.     ' ',
  998.   '   Then:' , 
  999.   '              Press the "Enter" key to execute the edit button callback.' , 
  1000.   '   Or:',    
  1001.   '              Press the "Next" button.',
  1002. ' '};
  1004.   slide(idx).info = 'addlift';
  1005.   
  1006.   %========== Slide 31 ==========
  1007.   idx = idx+1;
  1008.   slide(idx).code = {
  1009.     'img = findall(ax(4:7),''type'',''image''); delete(img);',
  1010.     'keptDEF = 10;',
  1011.     'keptCFS = str2num(get(ediCFS,''String''));',
  1012.     ['if isempty(keptCFS) | keptCFS<0 | keptCFSMAX<keptCFS | keptCFS~=fix(keptCFS)' ...
  1013.            ' , keptCFS = keptDEF; end'],
  1014. 'set(ediCFS,''String'',int2str(keptCFS));',      
  1015.     'demolift(''Local_FUN'',''set_UI_Graphics'',5);', 
  1016.     'demolift(''Local_FUN'',''plot_DECOMP'',x,xDEC,LSnew,level,ax(4:7),keptCFS,''1D'');',
  1017. '' };
  1019.   slide(idx).text = {
  1020.   ' ',
  1021. '   Choose the number of kept coefficients in the edit then ',
  1022.   '   Press the "Enter" key to execute the edit button callback.' , 
  1023.   ' ',    
  1024.   '   Or',
  1025.   ' ',    
  1026.   '   Press the "Prev" button.',
  1027. '   Choose the number of kept coefficients in the edit.',
  1028.   '   Then press the "Next" button.',
  1029. ' '};
  1031.   slide(idx).info = 'addlift';
  1032.   
  1033.   %========== Slide 32 ==========
  1034.   idx = idx+1;
  1035.   slide(idx).code = {
  1036.     'S1  = ''LIFTING (22): 2D ANALYSIS'';',
  1037.     'S2  = ''-------------------------'';',
  1038.     'S3  = blanc;',
  1039.     'S4  = ''  Like for one dimensional signal, you can perform '';',
  1040.     'S5  = ''  2D analysis using Lifting Tools.'';',
  1041.     'S6  = ''  Look at LWT2, ILWT2, LWTCOEF2 functions to get more information.'';',
  1042.     'S7  = blanc;',
  1043.     'S8  = ''  As an example:'';',  
  1044.     'S9  = ''      Load an image (woman)'';', 
  1045.     'S10 = ''      And perform a 2D wavelet analysis at level 8'';',  
  1046.     'S11 = ''      using the Haar wavelet.'';',  
  1047.     'S12 = blanc;',
  1048.     'S13 = '' '';',
  1049.     'Infostr = strvcat(S1,S2,S3,S4,S5,S6,S7,S8,S9,S10,S11,S12,S13);',
  1050.     
  1051. 'set(t(1),''FontSize'',10,''String'',Infostr);',
  1052.     'demolift(''Local_FUN'',''set_UI_Graphics'',0);', 
  1053. '' };
  1055.   slide(idx).text = {
  1056. ' % Compute the decomposition at level 8.',
  1057.   ' ',
  1058.     '       LSdb1 = liftwave(''db1''); level2D = 8;',
  1059.     '       St = load(''woman.mat'');',
  1060. '       xDEC = lwt2(St.X,LSdb1,level2D);',
  1061.   ' ',    
  1062. ' % And then view the decomposition.',
  1063.   ' ',
  1064. '      Choose the number of kept coefficients in the edit then ',
  1065.   '      Press the "Enter" key to execute the edit button callback.' , 
  1066.   ' ',    
  1067.   ' ',
  1068. ' '};
  1070.   slide(idx).info = 'lwt2';
  1072.   %========== Slide 32 ==========
  1073.   idx = idx+1;
  1074.   slide(idx).code = {
  1075.     'lin = findall(ax(4:7),''type'',''line''); delete(lin);',
  1076.     'LSdb1 = liftwave(''db1''); level2D = 8;',
  1077.     'St = load(''woman.mat'');',
  1078. 'xDEC = lwt2(St.X,LS,level2D);',
  1079. 'keptCFSMAX = prod(size(St.X));',
  1080. 'set(ediCFS,''String'',int2str(keptCFSMAX),''Enable'',''On'');',
  1081.     'demolift(''Local_FUN'',''set_UI_Graphics'',6);', 
  1082.     'demolift(''Local_FUN'',''plot_DECOMP'',St.X,xDEC,LS,level2D,ax(4:7),keptCFSMAX,''2D'');',
  1083. '' };
  1085.   slide(idx).text = {
  1086.   ' ',
  1087. '      Choose the number of kept coefficients in the edit then ',
  1088.   '      Press the "Enter" key to execute the edit button callback.' , 
  1089.   ' ',    
  1090.   ' ',
  1091. ' '};
  1093.   slide(idx).info = 'lwt2';
  1094.   
  1095.   %==============================
  1097.   
  1098.   varargout{1} = slide;
  1099.   
  1100. end
  1102. %--------------------------------------------------------------
  1103. function set_UI_Graphics(axMode)
  1105. hdl_UI_Graphics = wtbxappdata('get',gcf,'hdl_UI_Graphics');
  1106. if isempty(hdl_UI_Graphics)
  1107.     pos = [0.0200  0.3700  0.7700  0.6000 ;  ...
  1108.            0.0750  0.4327  0.2848  0.4973 ;  ...
  1109.            0.4652  0.4327  0.2848  0.4973 ;  ...
  1110.            0.1300  0.4527  0.7750  0.1297 ;  ...
  1111.            0.1300  0.2813  0.7750  0.1297 ;  ...
  1112.            0.1300  0.1100  0.7750  0.1297 ;  ...
  1113.            0.1300  0.6240  0.7750  0.1297 ;  ...
  1114.            0.1300  0.7953  0.7750  0.1297    ...
  1115.            ];
  1116.     ax(1) = axes('Position',pos(1,:), ...
  1117.         'Xtick',[],'Ytick',[],'Color',[1 1 1],'Box','On');
  1118.     rectangle('LineWidth',10,'EdgeColor',[1 0.75 0.75]);
  1119.     t(1) = text(0.04,0.95,' ', ...
  1120.         'Interpreter','tex','VerticalAlignment','top',...
  1121.         'Color',[0 0 0],'tag','LiftDEMO');
  1122.     for k = 2:size(pos,1)
  1123.         ax(k) = axes('Position',pos(k,:),'Visible','Off');
  1124.     end
  1125.     %-------------------------------------------------------
  1126.     hdl_UI_Graphics.mode = axMode; 
  1127.     hdl_UI_Graphics.ax = ax; 
  1128.     hdl_UI_Graphics.t  = t;
  1129.     wtbxappdata('set',gcf,'hdl_UI_Graphics',hdl_UI_Graphics);
  1130. else
  1131.     axModeOLD = hdl_UI_Graphics.mode;
  1132.     hdl_UI_Graphics.mode = axMode; 
  1133.     ax = hdl_UI_Graphics.ax;
  1134.     t  = hdl_UI_Graphics.t;
  1135.     pos = get(ax,'Position'); pos = cat(1,pos{:});
  1136. end
  1138. switch axMode
  1139.     case -1
  1140.         setPOS = [];
  1141.         vis ={'Off','Off','Off','Off','Off','Off','Off','Off'};
  1142.         
  1143.     case 0
  1144.         setPOS = [];
  1145.         vis ={'On','Off','Off','Off','Off','Off','Off','Off'};
  1146.         
  1147.     case 1
  1148.         setPOS = [];
  1149.         vis ={'Off','On','On','Off','Off','Off','Off','Off'};
  1150.         
  1151.     case 2
  1152.         setPOS = [4:6]; 
  1153.         pos(setPOS,:) = [...
  1154.                 0.07  0.80  0.70  0.15
  1155.                 0.07  0.60  0.70  0.15
  1156.                 0.07  0.40  0.70  0.15
  1157.                 ];
  1158.         vis ={'Off','Off','Off','On','On','On','Off','Off'}; 
  1159.         
  1160.     case 3
  1161.         setPOS = [4:7];         
  1162.         pos(setPOS,:) = [  ...
  1163.                 0.07  0.85  0.70  0.11
  1164.                 0.07  0.70  0.70  0.11
  1165.                 0.07  0.55  0.70  0.11
  1166.                 0.07  0.40  0.70  0.11
  1167.                 ];
  1168.         vis ={'Off','Off','Off','On','On','On','On','Off'}; 
  1169.         
  1170.     case 4
  1171.         setPOS = [4:8];         
  1172.         pos(setPOS,:) = [  ...
  1173.                 0.07  0.88  0.70  0.09
  1174.                 0.07  0.76  0.70  0.09
  1175.                 0.07  0.64  0.70  0.09
  1176.                 0.07  0.52  0.70  0.09
  1177.                 0.07  0.40  0.70  0.09
  1178.                 ];
  1179.         vis ={'Off','Off','Off','On','On','On','On','On'}; 
  1180.         
  1181.     case {5,6}
  1182.         setPOS = [4:7];                 
  1183.         pos(setPOS,:) = [  ...
  1184.                 0.06  0.72  0.30  0.23                 
  1185.                 0.45  0.72  0.30  0.23
  1186.                 0.06  0.41  0.30  0.23
  1187.                 0.45  0.41  0.30  0.23
  1188.                 ];
  1189.         vis ={'Off','Off','Off','On','On','On','On','Off'}; 
  1190. end
  1192. switch axMode
  1193.     case {2,3,4}
  1194.         hdl_DEN = findobj(ax(4),'type','line','Tag','denoSIG');
  1195.         if ~isempty(hdl_DEN), delete(hdl_DEN); end
  1196.     case 5 , set(ax(4),'Ydir','Normal');
  1197.     case 6 , set(ax(4),'Ydir','Reverse');
  1198. end
  1200. for j=1:length(setPOS)
  1201.     k = setPOS(j);
  1202.     set(ax(k),'Position',pos(k,:));
  1203. end
  1204. for k = 1:size(pos,1)
  1205.     set(findall(ax(k)),'Visible',vis{k});
  1206. end
  1207. %--------------------------------------------------------------
  1208. function plot_phi_psi(LS,ax)
  1210. N = 10;
  1211. phi = zeros(1,16); phi(9) = 1;
  1212. for k = 1:N , phi = dyadup(phi,0); phi = ilwt(phi,LS,1); end;
  1213. phi = (sqrt(2)^N)*phi;
  1215. psi = zeros(1,16); psi(8) = 1;
  1216. psi = ilwt(psi,LS,1);
  1217. for k = 1:N, psi = dyadup(psi,0); psi = ilwt(psi,LS,1); end;
  1218. psi = -(sqrt(2)^(N+1))*psi;
  1220. yLimPHI = [min(phi) max(phi)];
  1221. yLimPSI = [min(psi) max(psi)];
  1223. set(findall(ax),'Visible','Off')
  1224. set(ax(2:3),'Visible','On')
  1225. axes(ax(2)); plot(phi,'r','LineWidth',2); grid
  1226. set(ax(2),'Ylim',yLimPHI);
  1227. xlabel('Scaling Function');
  1228. axes(ax(3)); plot(psi,'g','LineWidth',2); grid
  1229. set(ax(3),'Ylim',yLimPSI);
  1231. xlabel('Wavelet');
  1232. set(ax(2:3),'Xlim',[1 length(psi)],'XtickLabel','')
  1233. %--------------------------------------------------------------
  1234. function plot_ANALYSIS(x,a,d,level,ax)
  1236. lx = length(x);
  1237. for k = 1:level
  1238.    axes(ax(1)) ; plot(x,'r') ; 
  1239.    title(''); yl = ylabel('S');
  1240.    p  = get(yl,'Position'); p(1) = -55;
  1241.    set(yl,'Rotation',0,'Position',p)
  1242.    set(ax(1),'Ylim',[min(x) max(x)]);
  1244.    axes(ax(2)) ; plot(a(k,:),'b') ; 
  1245.    yl = ylabel(['A' int2str(k)]);
  1246.    p  = get(yl,'Position'); p(1) = -55;
  1247.    set(yl,'Rotation',0,'Position',p)
  1248.    
  1249.    set(ax(2),'Ylim',[min(a(k,:)) max(a(k,:))]);
  1250.    for j =1:k
  1251.        idxAX = 3+k-j;
  1252.        axes(ax(idxAX)) ; plot(d(j,:),'g') ; 
  1253.        yl = ylabel(['D' int2str(j)]);
  1254.        p  = get(yl,'Position'); p(1) = -55;
  1255.        set(yl,'Rotation',0,'Position',p)
  1256.        set(ax(idxAX),'Ylim',[min(d(j,:)) max(d(j,:))]);
  1257.    end
  1258.    set(ax,'Xlim',[1 lx])
  1260. end
  1261. %--------------------------------------------------------------
  1262. function plot_DECOMP(x,xDEC,LS,level,ax,keptCFS,DIM)
  1264. switch DIM
  1265.     case '1D'
  1266.         NbMaxCFS = length(x);
  1267.         idxMIN = NbMaxCFS-keptCFS;
  1268.         T = sort(abs(xDEC(:)));
  1269.         
  1270.         xDEC_Thr = xDEC;
  1271.         if idxMIN>0 , thr = T(idxMIN); else , thr = 0; end 
  1272.         xDEC_Thr(abs(xDEC_Thr)<=thr) = 0;
  1273.         xREC = ilwt(xDEC_Thr,LS,level);
  1274.         
  1275.         yLX    = [min(x) max(x)];
  1276.         yLXREC = [min(xREC) max(xREC)];
  1277.         yLX_XREC = [min([yLX(1),yLXREC(1)]) , max([yLX(2),yLXREC(2)])];
  1278.         
  1279.         axes(ax(1)) ; plot(x,'r') ; title('Signal')
  1280.         hdl_DEN = findobj(ax(1),'type','line','Tag','denoSIG');
  1281.         if (idxMIN>0)
  1282.             if isempty(hdl_DEN)
  1283.                 axes(ax(1)) ; hold on; hdl_DEN = plot(xREC,'y','Tag','denoSIG');
  1284.             else
  1285.                 set(hdl_DEN,'Ydata',xREC);
  1286.             end
  1287.             set(ax(1),'Ylim',yLX_XREC);
  1288.         else
  1289.             delete(hdl_DEN)
  1290.             set(ax(1),'Ylim',yLX);
  1291.         end
  1292.         
  1293.         axes(ax(3)) ; plot(xDEC,'b')
  1294.         xlabel('Original Lifted Decomposition (in place)','Color','k')
  1295.         
  1296.         axes(ax(4)) ; plot(xDEC_Thr,'b')
  1297.         xlabel('Modified Lifted Decomposition (in place)','Color','k')
  1298.         
  1299.         axes(ax(2)) ;  plot(xREC,'y') ; set(ax(2),'Ylim',yLXREC);
  1300.         strTITLE = ['Nb. Cfs ==> total: ' int2str(NbMaxCFS) ...
  1301.                 ' - kept: ' int2str(NbMaxCFS-idxMIN)];
  1302.         title(strTITLE,'Color','c');
  1303.         xlabel(['Threshold = ' , num2str(thr)],'Color','c')
  1304.         
  1305.         set(ax,'Xlim',[1,NbMaxCFS]);
  1306.         wtbxappdata('set',gcf,'analysis_PARAM_1D',{x,xDEC,LS,level,ax,keptCFS,DIM});        
  1307.         
  1308.     case '2D'
  1309. LSType  = [];
  1310. NBC = min(max(max(x)),256);
  1311. NbMaxCFS = prod(size(x));
  1312. idxMIN = NbMaxCFS-keptCFS;
  1313. T = sort(abs(xDEC(:)));
  1314. xDEC_Thr = xDEC;
  1315. if idxMIN>0 , thr = T(idxMIN); else , thr = 0; end 
  1316. xDEC_Thr(abs(xDEC_Thr)<=thr) = 0;
  1317. xREC = ilwt2(xDEC_Thr,LS,level);
  1319. set(gcf,'Colormap',pink(NBC));
  1320. axes(ax(1)) ; image(x); title('Original Image')
  1321.         set(ax(1),'Xlim',[1,size(x,2)],'Ylim',[1,size(x,1)])
  1322. axes(ax(3)); 
  1323. if isempty(LSType) , image(abs(xDEC)); else , imagesc(abs(xDEC)); end
  1324. xlabel('Original Lifted Decomposition (in place)','Color','k')
  1325. axes(ax(4));
  1326. if isempty(LSType) , image(abs(xDEC_Thr)); else , imagesc(abs(xDEC_Thr)); end
  1327. xlabel('Modified Lifted Decomposition (in place)','Color','k')
  1329. axes(ax(2)); image(xREC);
  1330. strTITLE = ['Nb. Cfs ==> total: ' int2str(NbMaxCFS) ...
  1331.                     ' - kept: ' int2str(NbMaxCFS-idxMIN)];
  1332. title(strTITLE,'Color','c');
  1333. xlabel(['Threshold = ' , num2str(thr)],'Color','c')
  1334.         wtbxappdata('set',gcf,'analysis_PARAM_2D',{x,xDEC,LS,level,ax,keptCFS,DIM});        
  1335. end
  1336. wtbxappdata('set',gcf,'analysis_DIM',DIM);
  1337. %--------------------------------------------------------------
  1338. function store_ANALYSIS(x,xDEC,LSnew,level,ax,keptCFS,DIM)
  1340. switch DIM
  1341.     case '1D' , wtbxappdata('set',gcf,'analysis_PARAM_1D',{x,xDEC,LSnew,level,ax,keptCFS,DIM});
  1342.     case '2D' , wtbxappdata('set',gcf,'analysis_PARAM_1D',{x,xDEC,LSnew,level,ax,keptCFS,DIM});
  1343. end
  1344. wtbxappdata('set',gcf,'analysis_DIM',DIM);
  1345. %--------------------------------------------------------------
  1346. function edit_CB
  1348. DIM = wtbxappdata('get',gcf,'analysis_DIM');
  1349. switch DIM
  1350.     case '1D' , analysis_PARAM = wtbxappdata('get',gcf,'analysis_PARAM_1D');
  1351.     case '2D' , analysis_PARAM = wtbxappdata('get',gcf,'analysis_PARAM_2D');
  1352. end
  1353. keptCFS = str2num(get(gcbo,'String'));
  1354. if isempty(keptCFS)
  1355.     keptCFS = analysis_PARAM{end-1};
  1356.     set(gcbo,'String',int2str(keptCFS));
  1357. end
  1358. analysis_PARAM{end-1} = keptCFS;
  1359. demolift('Local_FUN','plot_DECOMP',analysis_PARAM{:});  
  1360. %--------------------------------------------------------------
  1361. function S = wlift_str_util(option,S)
  1362. %WLIFT_STR_UTIL String utilities for lifting demos. 
  1363. switch option
  1364.     case 'replace'
  1365.         idx = find(S=='('); S(idx(2:end)) = '{'; 
  1366.         idx = find(S==')'); S(idx(2:end)) = '}';
  1367.         idx = find(S=='*'); S(idx) = ' ';
  1368.     case 'replaceMAT'
  1369.         idx = find(S=='('); S(idx) = '{';
  1370.         idx = find(S==')'); S(idx) = '}';
  1371.         idx = find(S=='*'); S(idx) = ' ';
  1372.         idx = find(S=='|'); S(idx) = ' ';
  1373.     case 'replaceMAT_1'
  1374.         idx = find(S=='('); S(idx) = '{';
  1375.         idx = find(S==')'); S(idx) = '}';
  1376.         idx = find(S=='*'); S(idx) = ' ';
  1377.     case 'replaceMAT_2'
  1378.         idx = find(S=='('); S(idx) = '{';
  1379.         idx = find(S==')'); S(idx) = '}';
  1380.     case 'replacePOW'
  1381.         idx_1 = find(S=='^'); 
  1382.         idx_2 = find(S=='('); idx_3 = find(S==')'); 
  1383.         for k = 1:length(idx_1)
  1384.             j = idx_1(k);
  1385.             if isequal(S(j+1),'(')
  1386.                 S(j+1) = '{';
  1387.                 tmp = find(idx_3>j);
  1388.                 jj = idx_3(tmp(1));
  1389.                 S(jj) = '}';
  1390.             end
  1391.         end
  1392.         idx = find(S=='*'); S(idx) = ' ';
  1393. end
  1394. %--------------------------------------------------------------