KSVD图像去噪

SAR图像去噪，希望对刚接触图像处理的同学又用吧，代码使用matlab写的

function [qXSNewInd]=AddOneElement（OmegaQXSYopt）
%AddOneElement adds one entry to the co-support updates the signal 
%estimate and the accumulated orthogonal set.
%  [qXSNewInd]=AddOneElement（OmegaQXSYopt）
%  ========================================================================
%  Input:
%  Omega - analysis dictionary.
%  Q     - matrix of the orthogonal set accumulated so far （in its rows）.
%  X     - the current estimate for the signal.
%  S     - the current estimate for the co-support.
%  Y     - the observed （possibly noisy） signal.
%  opt   - the type of optimality criterion used for adding elements to the 
%          co-support: 0 - THR 1 - BG 2 - OBG. 
%  ========================================================================
%  Output:
%  q      - the new row vector added to the orthogonal set.
%  X      - the new estimate for the signal.
%  S      - the new estimate for the co-support.
%  NewInd - the new element added to the co-support.
%  ========================================================================
%  Tomer Peleg
%  Department of Electrical Engineering
%  Technion Haifa 32000 Israel
%  tomerfa@tx.technion.ac.il
%
%  October 2012
%  ========================================================================
d=size（Omega2）;
inds=find（S==0）;
q=[];
NewInd=[];
if ~isempty（inds）
    Q_temp=Omega（inds:）;
    if opt==2
        if ~isempty（Q）
            % modified gram-schmidt for all possible directions
            for j=1:size（Q1）
                Q_temp=Q_temp-（Q_temp*Q（j:）‘）*Q（j:）;
            end
            Q_temp=Q_temp./repmat（sqrt（diag（Q_temp*Q_temp‘））1d）;
        end
    end
    % project on each of the possible directions
    if opt==0
        Proj=（Q_temp*Y）.^2;
    else
        Proj=（Q_temp*X）.^2;
    end
    % choose the direction with the smallest energy reduction
    [dummypos] = min（Proj）;
    q=Q_temp（pos:）;
    if opt<2
        if ~isempty（Q）        
            for j=1:size（Q1）
                q=q-（q*Q（j:）‘）*Q（j:）;
            end
            q=q/sqrt（q*q‘）;
        end
    end
    if nargout>1
        X=X-q‘*（q*X）;
    end
    if nargout>2
        S=abs（Omega*X）<1e-4;
    end
    if nargout>3
        NewInd=inds（pos）;
    end
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件     185727  2001-07-12 13:02  images\barbara.png

     文件     177762  2001-07-12 12:42  images\boats.png

     文件      34985  2005-09-11 14:44  images\house.png

     文件     151199  2009-10-04 09:08  images\lena.png

     文件      40181  2002-08-29 16:48  images\peppers256.png

     文件       8601  2011-11-01 15:21  images\PWCImage.png

     文件       6294  2012-10-28 16:19  AnalysisKSVD.m

     文件        703  2012-10-28 15:54  readme.txt

     文件       3115  2012-10-28 15:54  Contents.m

     文件       3611  2012-10-28 15:20  TestKSVDNaturalImages.m

     文件       3405  2012-10-28 15:15  TestKSVDPWCImage.m

     文件       1645  2012-10-28 15:12  TestKSVDRandomOmega.m

     文件       2283  2012-10-28 15:11  TestKSVDOmegaDIF.m

     文件       4049  2012-10-28 15:11  TestPursuitOmegaDIF.m

     文件       2243  2012-10-28 15:01  GenerateAnalysisSignals.m

     文件       1346  2012-10-28 15:00  GenerateOmegaDIF.m

     文件        966  2012-10-28 14:59  GenerateRandomRow.m

     文件       1063  2012-10-28 14:59  ComputeOrthoSet.m

     文件       2304  2012-10-28 14:59  AddOneElement.m

     文件       2064  2012-10-28 14:58  ErrorBGP.m

     文件       1898  2012-10-28 14:58  RankBGP.m

     文件       1268  2012-10-28 14:58  DisplayOmega.m

     文件       2117  2012-10-27 09:19  dictdist.m

     文件       7680  2011-05-16 09:08  col2imstep.mexw32

     文件       7168  2011-05-16 09:08  im2colstep.mexw32

     文件       3924  2011-05-16 08:48  col2imstep.c

     文件       1001  2011-05-16 08:48  col2imstep.m

     文件       1178  2011-05-16 08:48  countcover.m

     文件       3388  2011-05-16 08:48  im2colstep.c

     文件       1341  2011-05-16 08:48  im2colstep.m

............此处省略16个文件信息