Commit 655e54fd authored by Franco Buffa's avatar Franco Buffa
Browse files

source and data

parent 5e6955d0

dev/errmap.m

+1 −1
Original line number Original line Diff line number Diff line
@@ -142,7 +142,7 @@ function errmap(fconf) 
 axis('square');

 axis('square');

 colormap(cmap);

 colormap(cmap);

 colorbar;

 colorbar;

 title('error map (mm)');

 title('deflection map [mm]');

%print -dpng errmap.png

%print -dpng errmap.png

%   figure(9);

%   figure(9);

%   surf(u1);

%   surf(u1);

dev/map.m

+1 −0
Original line number Original line Diff line number Diff line 
function map(s0,ca,x1,lutname,fname,mp)

function map(s0,ca,x1,lutname,fname,mp)

%

%

%   map(s0,ca,x1,lutname,fname,mp)

%   lutname='axist_parab.txt';

%   lutname='axist_parab.txt';

%   fname='map';

%   fname='map';

%   s0=1;

%   s0=1;

dev/mresize.m

+18 −1
Original line number Original line Diff line number Diff line 
function m2=mresize(m1,mp,ca)

function m2=mresize(m1,r0,r1,mp,ca)

 [b blk]=mask(r0,r1,mp);

 n0=rows(m1);

 n0=rows(m1);

 xx=1:n0;

 xx=1:n0;

 x=linspace(1,n0,mp);

 x=linspace(1,n0,mp);

 y=x';

 y=x';

 m2=interp2(xx,xx,m1,x,y);

 m2=interp2(xx,xx,m1,x,y);

 m2=m2.*b;

 imagesc(m2);

 imagesc(m2);

 caxis([-ca ca]);

 caxis([-ca ca]);

 axis('square');

 axis('square');
@@ -11,3 +13,18 @@ function m2=mresize(m1,mp,ca) 
 colorbar;

 colorbar;

end

end

 

 

function [b blk]=mask(r0,r1,mp)

 b=ones(mp,mp);

 c=ones(mp,mp).*(1:mp);

 r=transpose(c);

 if(r0)

   i1=find(sqrt((r-mp/2).^2+(c-mp/2).^2)<=r0);

   b(i1)=0;

 end  

 if(r1)

  i1=find(sqrt((r-mp/2).^2+(c-mp/2).^2)>r1);

  b(i1)=0;

 end

 blk=find(b==0);

%

end
 
 No newline at end of file

dev/radrms.m

+10 −5
Original line number Original line Diff line number Diff line
@@ -31,6 +31,12 @@ function radrms(fconf) 
 end

 end





%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%



 if (~exist('dr','var'))

  aber=0;

  dr=5;

 end



 load zern.dat

 load zern.dat





%   i0=find(abs(u3)>tol4);

%   i0=find(abs(u3)>tol4);
@@ -39,11 +45,10 @@ function radrms(fconf) 
 x0=((1:mp)-mp/2);

 x0=((1:mp)-mp/2);

 [xm ym]=meshgrid(x0,x0);

 [xm ym]=meshgrid(x0,x0);

 R=sqrt(xm.^2+ym.^2); 

 R=sqrt(xm.^2+ym.^2); 

 d=5;

 x0=y=u=[];

 x0=y=u=[];

 for r=r0:r1-d

 for r=r0:r1-dr

  z=(R>=r);

  z=(R>=r);

  z=z.*(R<r+d);

  z=z.*(R<r+dr);

  i0=find(z>0);

  i0=find(z>0);

%    imagesc(z)

%    imagesc(z)

  y=[y ;std(u3(i0))];

  y=[y ;std(u3(i0))];
@@ -52,8 +57,8 @@ function radrms(fconf) 
 end  

 end  

 dx=D/(mp-1)/k0;

 dx=D/(mp-1)/k0;

 x=x0*dx;

 x=x0*dx;

 y(end-d:end)=y(end-d);

 y(end-dr:end)=y(end-dr);

 u(end-d:end)=u(end-d);

 u(end-dr:end)=u(end-dr);

 

 

 figure(1);

 figure(1);

 plot(x,y,'*-');

 plot(x,y,'*-');

dev/unwr2.m

0 → 100644
+47 −0
Original line number Original line Diff line number Diff line 
function unwr2(fconf)



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%

% Function:  unwr2.m

% Package:   HoP

% Author:    Franco Buffa

% Purpose:   Driver for: unwrap_phase

% Date:      

%	

% GNU Octave Version: 4.0.2

%

% revised by 

%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%

% unwr2.m  is a function that 

%

% Output:

% 

%

% Input:

%

%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%

 if nargin>0

  run(fconf); 

 else

  run('conf.m');

 end

 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% 

 clear all

 load phase.dat

 e2(blk)=NaN;

 u0=unwrap_phase(e2);

 u0(blk)=0; 

 figure(1);

 imagesc(u0);

 axis('square');

 colormap('gray');

 colorbar;

 title('unwrapped map');

 save unwr.dat u0 blk

end
 
 No newline at end of file