Commit 55c91e0a authored by Franco Buffa's avatar Franco Buffa
Browse files

source and data

parent 6c71b951

dev/phase.m

+2 −1
Original line number Diff line number Diff line
@@ -56,7 +56,8 @@ function phase(fconf) 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 



 r=abs(y); %sqrt(y.*conj(y));

 f=[r(m/2+1:m,m/2+1:m) r(m/2+1:m,1:m/2); r(1:m/2,m/2+1:m) r(1:m/2,1:m/2)];

%f=[r(m/2+1:m,m/2+1:m) r(m/2+1:m,1:m/2); r(1:m/2,m/2+1:m) r(1:m/2,1:m/2)];

 f=fftshift(r);

 

 sc=sr=0;

 gf=~(pc&pr);

dev/uv.m

+15 −5
Original line number Diff line number Diff line
@@ -92,15 +92,25 @@ function uv(fconf) 
 xlabel('u');

 ylabel('v');

%print -dpng preproc.png



 pr0=mp/2-fix(size(er)(1)/2);

 pc0=mp/2-fix(size(er)(2)/2);

%

% 

 [pr0 pc0]=find(emod==0);

 [nr0 nc0]=size(emod);

 z=zeros(mp,mp);

 z(pr0:pr0+size(er)(1)-1,pc0:pc0+size(er)(2)-1)=er;

 z(mp/2-pr0:mp/2-pr0+nr0-1,mp/2-pc0:mp/2-pc0+nc0-1)=er;

 er=z;

 z=zeros(mp,mp);

 z(pr0:pr0+size(ei)(1)-1,pc0:pc0+size(ei)(2)-1)=ei;

 z(mp/2-pr0:mp/2-pr0+nr0-1,mp/2-pc0:mp/2-pc0+nc0-1)=ei;

 ei=z;



%   pr0=mp/2-fix(size(er)(1)/2);

%   pc0=mp/2-fix(size(er)(2)/2);

%   z=zeros(mp,mp);

%   z(pr0:pr0+size(er)(1)-1,pc0:pc0+size(er)(2)-1)=er;

%   er=z;

%   z=zeros(mp,mp);

%   z(pr0:pr0+size(ei)(1)-1,pc0:pc0+size(ei)(2)-1)=ei;

%   ei=z;

%

 save uuvv.dat er ei

%