← Back to branch summary

~gerald-mwangi/+junk/ThesisMatlab

« back to all changes in this revision

Viewing changes to Matlab/StructureTensFeaturesAnal/chambollePockTVPrimalDual.m

  • Committer: gerald.mwangi at gmx
  • Date: 2018-03-21 13:31:26 UTC
  • Revision ID: gerald.mwangi@gmx.de-20180321133126-89jt0jr5vg7n8o3u
work on primal dual code

Show diffs side-by-side

added added

removed removed

Lines of Context:
Matlab/StructureTensFeaturesAnal/chambollePockTVPrimalDual.m
15
15
xgradf2=[gradf]%[-1 1];
16
16
ygradf2=[gradf ]'%[-1 1]';
17
17
theta=1.0;
18
 
taudual=0.75;%1.0/(4.0*lambda);
 
18
taudual=0.5;%1.0/(4.0*lambda);
19
19
tauprimal=0.1;
20
20
tauspace=0.0;%tau/5.0;
21
21
px=zeros(imsize);%imfilter(im,[gradf,0],'conv','replicate');
74
74
                    div_p0 = cxx+cyy;
75
75
            lagim=lambda_data*(ims-im0s)-lambda*div_p0;
76
76
            
 
77
%             u=ims-tauprimal*lambda*div_p0;
77
78
                 imsnew=ims-tauprimal*lagim;
78
79
        
79
 
        ims=imsnew;           
 
80
        ims=imsnew+theta*(imsnew-ims);           
80
81
            
81
82
            
82
83
%                         dx=imfilter(ims,[0.5 0 0.5],'conv','replicate');
137
138
     curv=(cxx-cyy).^2+(cxy+cyx).^2;
138
139
     intrinscurv=[intrinscurv,sum(sum(curv))];
139
140
 
140
 
     figure(2);
 
141
     figure(4);
141
142
 
142
143
    im=scale_image(ims,0,255);
143
144
     imshow([uint8(scale_image(im,0,255)),uint8(im0)]);