Changeset 1129 for trunk/src/uvmat.m

Timestamp:

Mar 1, 2024, 8:09:20 PM (5 months ago)

Author:

sommeria

Message:

LIF calibration update and extrct_multitif

File:

• trunk/src/uvmat.m (modified) (6 diffs)

trunk/src/uvmat.m

@@ -1657 +1657 @@
 
 %% rescale the image
-[nby,nbx]=size(UvData.Field.A);
-x=linspace(UvData.Field.Coord_x(1),UvData.Field.Coord_x(2),nbx)-nbx/2;
-y=linspace(UvData.Field.Coord_y(1),UvData.Field.Coord_y(2),nby)-nby/2;
-[X,Y]=meshgrid(x,y);
+% [nby,nbx]=size(UvData.Field.A);
+% x=linspace(UvData.Field.Coord_x(1),UvData.Field.Coord_x(2),nbx)-nbx/2;
+% y=linspace(UvData.Field.Coord_y(1),UvData.Field.Coord_y(2),nby)-nby/2;
+% [X,Y]=meshgrid(x,y);
 %coeff_quad=0.15*4/(nbx*nbx);% image luminosity reduced by 10% at the edge
 %UvData.Field.A=double(UvData.Field.A).*(1+coeff_quad*(X.*X+Y.*Y));
@@ -1686 +1686 @@
     ylabel('radius from light source')
     title('ref line in polar coordinates')
-    azimuth_ima=linspace(DataPol.Coord_y(1),DataPol.Coord_y(2),size(DataPol.A,1))-360;%array of angular indices on the transformed image
-    dist_source = interp1(theta_ref,r_ref,azimuth_ima);% get the polar position of the reference line
-    dist_source_pixel=round(size(DataPol.A,2)*(dist_source-DataPol.Coord_x(1))/(DataPol.Coord_x(2)-DataPol.Coord_x(1)));
+    %azimuth_ima=linspace(DataPol.Coord_y(1),DataPol.Coord_y(2),size(DataPol.A,1));%array of angular indices on the transformed image
+   % dist_source = interp1(theta_ref,r_ref,azimuth_ima);% get the polar position of the reference line
+    dist_source = interp1(theta_ref,r_ref,DataPol.theta);% get the polar position of the reference line
+    dist_source_pixel=floor(size(DataPol.A,2)*(dist_source-DataPol.radius(1))/(DataPol.radius(end)-DataPol.radius(1)))+1;
     line_nan= isnan(dist_source);
     dist_source_pixel(line_nan)=1;
@@ -1710 +1711 @@
 %% get the image A*r
 [npy,npx]=size(Anorm);
-AX=DataPol.Coord_x;
-AY=DataPol.Coord_y;
-dX=(AX(2)-AX(1))/(npx-1);
-dY=(AY(1)-AY(2))/(npy-1);%mesh of new pixels
-[R,Theta]=meshgrid(linspace(AX(1),AX(2),npx),linspace(AY(1),AY(2),npy));
+AX=DataPol.radius;
+AY=DataPol.theta;
+% dX=(AX(2)-AX(1))/(npx-1);
+% dY=(AY(1)-AY(2))/(npy-1);%mesh of new pixels
+[R,Theta]=meshgrid(linspace(AX(1),AX(end),npx),linspace(AY(1),AY(end),npy));%matrix of radius and angles with the same size as DataPol
 A=R.*Anorm;
 A=(A<=0).*ones(npy,npx)+(A>0).*A; %replaces zeros by ones
@@ -1734 +1735 @@
 gamma_coeff=NaN(1,npy);
 fitlength=NaN(1,npy);
+%[ThetaMask,RMask] = cart2pol(MaskData.Coord(:,1)-x0,MaskData.Coord(:,2)-y0);
+%ThetaMask=ThetaMask*180/pi
 for iY=1:npy% loop on the y index of the image in polar coordinate
     ALine=A(iY,:);%profile of image luminosity log (vs radial index)
@@ -1752 +1755 @@
     p = polyfit(RLine,ALine,1);
     gamma_coeff(iY)=-p(1);
-    fitlength(iY)=numel(find(check_good));
-end
+   fitlength(iY)=numel(find(check_good));
+end
+
+%% plot the decay coef versus theta
+figure(13)
+plot(Theta(:,1),100*gamma_coeff)
+xlabel('angle(degree)')
+ylabel('decay coeff(m-1)')
+
 
 %% show and store the rescaled image
@@ -1760 +1770 @@
 view_field(DataPol);
 imwrite(DataPol.A,NewImageName,'BitDepth',16)
-figure(13)
-plot(Theta(:,1),100*gamma_coeff)
-xlabel('angle(degree)')
-ylabel('decay coeff(m-1)')
+
 
 %% keep the average of gamma_coeff