Changeset 227 for trunk/src/set_object.m

Timestamp:

Mar 31, 2011, 1:42:51 PM (13 years ago)

Author:

sommeria

Message:

add function sub_field_series to apply the sub_field operation to a series of fileds (for instance subtracting a background to an image series)

File:

• trunk/src/set_object.m (modified) (4 diffs)

trunk/src/set_object.m

@@ -766 +766 @@
 
 %% plot the field projected on the object and store in the corresponding figue
-ProjData= proj_field(UvData.Field,ObjectData);%project the current interface field on ObjectData
+[ProjData,errormsg]= proj_field(UvData.Field,ObjectData);%project the current interface field on ObjectData
+if ~isempty(errormsg)
+    msgbox_uvmat('ERROR', errormsg)
+    return
+end
 PlotParam=read_plot_param(PlotHandles);
 [PlotType,Object_out{IndexObj}.PlotParam,plotaxes]=plot_field(ProjData,plotaxes,PlotParam);%update an existing field plot
@@ -846 +850 @@
 end
 msgbox_uvmat('CONFIRMATION',[answer{1}  ' saved'])
-%------------------------------------------------------------------------
+
 %------------------------------------------------------------------------
 % --- Executes on slider movement.
@@ -853 +857 @@
 Z_value=get(handles.z_slider,'Value');
 %rotation angles
-PlaneAngle(1)=str2num(get(handles.Phi,'String'));%first  angle in degrees
-PlaneAngle(2)=str2num(get(handles.Theta,'String'));%second  angle in degrees
-PlaneAngle(3)=str2num(get(handles.Psi,'String'));%second  angle in degrees
+PlaneAngle=[0 0 0]; 
+norm_plane=[0 0 1];
+cos_om=1;
+sin_om=0;
+
+PlaneAngle(1)=str2double(get(handles.Phi,'String'));%first  angle in degrees
+PlaneAngle(2)=str2double(get(handles.Theta,'String'));%second  angle in degrees
+PlaneAngle(3)=str2double(get(handles.Psi,'String'));%second  angle in degrees
+PlaneAngle=(pi/180)*PlaneAngle;
 om=norm(PlaneAngle);%norm of rotation angle in radians
-OmAxis=PlaneAngle/om; %unit vector marking the rotation axis
-cos_om=cos(pi*om/180);
-sin_om=sin(pi*om/180);
-coeff=OmAxis(3)*(1-cos_om);
-%components of the unity vector norm_plane normal to the projection plane
-norm_plane(1)=OmAxis(1)*coeff+OmAxis(2)*sin_om;
-norm_plane(2)=OmAxis(2)*coeff-OmAxis(1)*sin_om;
-norm_plane(3)=OmAxis(3)*coeff+cos_om;
+if isequal(om,0)
+    norm_plane=[0 0 1];
+else
+    OmAxis=PlaneAngle/om; %unit vector marking the rotation axis
+    cos_om=cos(om);
+    sin_om=sin(om);
+    coeff=OmAxis(3)*(1-cos_om);
+    %components of the unity vector norm_plane normal to the projection plane
+    norm_plane(1)=OmAxis(1)*coeff+OmAxis(2)*sin_om;
+    norm_plane(2)=OmAxis(2)*coeff-OmAxis(1)*sin_om;
+    norm_plane(3)=OmAxis(3)*coeff+cos_om;
+end
 
 %set new plane position and update graph
@@ -871 +885 @@
 set(handles.ZObject,'String',num2str(norm_plane(3)*Z_value,4))
 PLOT_Callback(hObject, eventdata, handles)
-%------------------------------------------------------------------------
+
 %------------------------------------------------------------------------
 % --- Executes on button press in HELP.
 function HELP_Callback(hObject, eventdata, handles)
+%------------------------------------------------------------------------
 path_to_uvmat=which ('uvmat');% check the path of uvmat
 pathelp=fileparts(path_to_uvmat);