← Back to branch summary

~gerald-mwangi/+junk/ThesisMatlab

« back to all changes in this revision

Viewing changes to Matlab/ThesisCode/determinantNoiseAnalysis/noiseDeterminant.m

  • Committer: gerald.mwangi at gmx
  • Date: 2016-05-25 18:10:22 UTC
  • Revision ID: gerald.mwangi@gmx.de-20160525181022-6v3d34cdbr1438it
started noise analysis2

Show diffs side-by-side

added added

removed removed

Lines of Context:
Matlab/ThesisCode/determinantNoiseAnalysis/noiseDeterminant.m
18
18
 
19
19
dx=imfilter(image,xgradf,'conv','replicate');
20
20
dy=imfilter(image,ygradf,'conv','replicate');
21
 
%[dx,dy]=gradient(image);
 
21
[dx,dy]=gradient(image);
22
22
nx=imfilter(noise,xgradf,'conv','replicate');
23
23
ny=imfilter(noise,ygradf,'conv','replicate');
 
24
[nx,ny]=gradient(noise);
24
25
%         figure(2);clf;
25
26
% %         imshow([uint8(im),uint8(im0),uint8(im-im0)]);
26
27
%         subplot(2,2,1); imagesc(dx); axis image; colormap(gray);
84
85
% dy2g=imfilter(dx2,GVyVy,'conv','replicate')-imfilter(dxdy,GVxVy,'conv','replicate').*2+imfilter(dy2,GVxVx,'conv','replicate');% Rotation component
85
86
% dxdyg=-imfilter(dx2,GVxVy,'conv','replicate')+imfilter(dxdy,GVxVx,'conv','replicate')-imfilter(dxdy,GVyVy,'conv','replicate')+imfilter(dy2,GVxVy,'conv','replicate');% mixture
86
87
% dx2g=imfilter(dx2,GVxVx,'conv','replicate')+2.*imfilter(dxdy,GVxVy,'conv','replicate')+imfilter(dy2,GVyVy,'conv','replicate');% Dilation component
87
 
 
 
88
dxg=imfilter(dx,G,'conv','replicate');
 
89
dyg=imfilter(dy,G,'conv','replicate');
88
90
dy2g=imfilter(dy2,G,'conv','replicate');
89
91
dxdyg=imfilter(dxdy,G,'conv','replicate');
90
92
dx2g=imfilter(dx2,G,'conv','replicate');
111
113
det3=2*dxnxg.*ny2g+2*nx2g.*dynyg-2*(nxnyg.*(dxnyg+dynxg));
112
114
det4=nx2g.*ny2g-nxnyg.^2;
113
115
 
 
116
%det3=(2*dxg.*dy2g-2*dxdyg.*dyg);%.*(dx.*dx+dy.*dy)+det.*2.*dx;
 
117
%det4=(2*dx2g.*dyg-2*dxdyg.*dxg);%.*(dx.*dx+dy.*dy)+det.*2.*dy;
114
118
% det1=dxnxg;
115
119
% det2=dxnyg;
116
120
% det3=dynxg;