Changeset 122 for trunk/src/uvmat.m
- Timestamp:
- Nov 14, 2010, 9:58:15 PM (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/uvmat.m
r121 r122 823 823 set(handles.transform_fct,'Value',1); % no transform by default 824 824 else 825 if (isfield(GeometryCalib,'R')&& ~isequal(GeometryCalib.R(2,1),0) && ~isequal(GeometryCalib.R(1,2),0)) | ...825 if (isfield(GeometryCalib,'R')&& ~isequal(GeometryCalib.R(2,1),0) && ~isequal(GeometryCalib.R(1,2),0)) ||... 826 826 (isfield(GeometryCalib,'kappa1')&& ~isequal(GeometryCalib.kappa1,0)) 827 827 set(handles.pxcm,'String','var') 828 828 set(handles.pycm,'String','var') 829 829 elseif isfield(GeometryCalib,'fx_fy') 830 pixcmx=GeometryCalib.fx_fy(1) %*GeometryCalib.R(1,1)*GeometryCalib.sx/(GeometryCalib.Tz*GeometryCalib.dpx);831 pixcmy=GeometryCalib.fx_fy(2) %*GeometryCalib.R(2,2)/(GeometryCalib.Tz*GeometryCalib.dpy);830 pixcmx=GeometryCalib.fx_fy(1);%*GeometryCalib.R(1,1)*GeometryCalib.sx/(GeometryCalib.Tz*GeometryCalib.dpx); 831 pixcmy=GeometryCalib.fx_fy(2);%*GeometryCalib.R(2,2)/(GeometryCalib.Tz*GeometryCalib.dpy); 832 832 set(handles.pxcm,'String',num2str(pixcmx)) 833 833 set(handles.pycm,'String',num2str(pixcmy)) … … 911 911 state_j='off'; %default 912 912 scan_option='i';%default 913 switch NomType 914 case {'_i_j','_i_j1-j2','_i1-i2_j','#_ab','#A','%01dA','%02dA','%03dA','%04dA'},% two navigation indices 913 NomTypeRaw=regexprep(NomType(2:end), '-', ''); 914 if numel(regexp(NomTypeRaw,'\D'))>=1 915 % case {'_i_j','_i_j1-j2','_i1-i2_j','#_ab','#A','%01dA','%02dA','%03dA','%04dA'},% two navigation indices 915 916 state_j='on'; 916 917 if isequal(nbfield,1) 917 scan_option='j'; 918 scan_option='j'; %scan j index by default if nbfield=1 918 919 end 919 920 end 920 921 if ~isempty(NomType_1) 921 switch NomType_1 922 case {'_i_j','_i_j1-j2','_i1-i2_j','#_ab'},% two navigation indices 922 NomTypeRaw=regexprep(NomType_1(2:end), '-', ''); 923 if numel(regexp(NomTypeRaw,'\D'))>=1 924 % switch NomType_1 925 % case {'_i_j','_i_j1-j2','_i1-i2_j','#_ab'},% two navigation indices 923 926 state_j='on'; 924 927 if isequal(nbfield,1) … … 1486 1489 ind_mask=findstr('mask',Name); 1487 1490 i=ind_mask-1; 1488 while val(i)==0 & i>01491 while val(i)==0 && i>0 1489 1492 i=i-1; 1490 1493 end … … 1668 1671 UvData=get(handles.uvmat,'UserData'); 1669 1672 1670 while get(handles.speed,'Value')~=0 & isequal(get(handles.Movie,'BusyAction'),'queue') % enable STOP command1673 while get(handles.speed,'Value')~=0 && isequal(get(handles.Movie,'BusyAction'),'queue') % enable STOP command 1671 1674 errormsg=runpm(hObject,eventdata,handles,increment); 1672 1675 if ~isempty(errormsg) … … 1902 1905 FileType='vol'; 1903 1906 else 1904 form=imformats(Ext( [2:end]));1907 form=imformats(Ext(2:end)); 1905 1908 if ~isempty(form)% if the extension corresponds to an image format recognized by Matlab 1906 1909 if isequal(NomType,'*'); … … 2059 2062 FieldName='image'; 2060 2063 else 2061 form=imformats(Ext( [2:end]));2064 form=imformats(Ext(2:end)); 2062 2065 if ~isempty(form)% if the extension corresponds to an image format recognized by Matlab 2063 2066 if isequal(NomType,'*'); … … 2106 2109 FieldName_1='image'; 2107 2110 else 2108 form=imformats(Ext_1( [2:end]));2111 form=imformats(Ext_1(2:end)); 2109 2112 if ~isempty(form)% if the extension corresponds to an image format recognized by Matlab 2110 2113 if isequal(NomType_1,'*'); … … 2623 2626 %use of mask 2624 2627 if isfield(ObjectData,'NbDim')&isequal(ObjectData.NbDim,2) 2625 if isfield(ObjectData,'Mask') & isfield(ObjectData,'A')2628 if isfield(ObjectData,'Mask') && isfield(ObjectData,'A') 2626 2629 flag_mask=double(ObjectData.Mask>200);%=0 for masked regions 2627 2630 AX=ObjectData.AX; … … 2629 2632 MaskX=ObjectData.MaskX; 2630 2633 MaskY=ObjectData.MaskY; 2631 if ~isequal(MaskX,AX)| ~isequal(MaskY,AY)2634 if ~isequal(MaskX,AX)||~isequal(MaskY,AY) 2632 2635 nxy=size(flag_mask); 2633 2636 sizpx=(ObjectData.MaskX(end)-ObjectData.MaskX(1))/(nxy(2)-1);%size of a mask pixel 2634 2637 sizpy=(ObjectData.MaskY(1)-ObjectData.MaskY(end))/(nxy(1)-1); 2635 x_mask= [ObjectData.MaskX(1):sizpx:ObjectData.MaskX(end)]; % pixel x coordinates for image display2636 y_mask= [ObjectData.MaskY(1):-sizpy:ObjectData.MaskY(end)];% pixel x coordinates for image display2638 x_mask=ObjectData.MaskX(1):sizpx:ObjectData.MaskX(end); % pixel x coordinates for image display 2639 y_mask=ObjectData.MaskY(1):-sizpy:ObjectData.MaskY(end);% pixel x coordinates for image display 2637 2640 %project on the positions of the scalar 2638 2641 npxy=size(ObjectData.A); 2639 2642 dxy(1)=(ObjectData.AY(end)-ObjectData.AY(1))/(npxy(1)-1);%grid mesh in y 2640 2643 dxy(2)=(ObjectData.AX(end)-ObjectData.AX(1))/(npxy(2)-1);%grid mesh in x 2641 xi= [ObjectData.AX(1):dxy(2):ObjectData.AX(end)];2642 yi= [ObjectData.AY(1):dxy(1):ObjectData.AY(end)];2644 xi=ObjectData.AX(1):dxy(2):ObjectData.AX(end); 2645 yi=ObjectData.AY(1):dxy(1):ObjectData.AY(end); 2643 2646 [XI,YI]=meshgrid(xi,yi);% creates the matrix of regular coordinates 2644 2647 flag_mask = interp2(x_mask,y_mask,flag_mask,XI,YI); … … 2650 2653 if isfield(ObjectData,'ListVarName') 2651 2654 for ilist=1:length(ObjectData.ListVarName) 2652 if isequal(ObjectData.ListVarName{ilist},'Mask')| isequal(ObjectData.ListVarName{ilist},'MaskX')|isequal(ObjectData.ListVarName{ilist},'MaskY')2655 if isequal(ObjectData.ListVarName{ilist},'Mask')||isequal(ObjectData.ListVarName{ilist},'MaskX')||isequal(ObjectData.ListVarName{ilist},'MaskY') 2653 2656 ind_off=[ind_off ilist]; 2654 2657 end … … 2841 2844 Tabchar={''};%default 2842 2845 Tabcell=[]; 2843 if isfield(Field{2},'ListGlobalAttribute')& ~isempty(Field{2}.ListGlobalAttribute)2846 if isfield(Field{2},'ListGlobalAttribute')&& ~isempty(Field{2}.ListGlobalAttribute) 2844 2847 for iline=1:length(Field{2}.ListGlobalAttribute) 2845 2848 Tabcell{iline,1}=Field{2}.ListGlobalAttribute{iline}; … … 3801 3804 index_code=get(handles.col_vec,'Value');% selected string index 3802 3805 col_code= list_code{index_code(1)}; % selected field 3803 if isequal(col_code,'black') | isequal(col_code,'white')3806 if isequal(col_code,'black') || isequal(col_code,'white') 3804 3807 set(handles.slider1,'Visible','off') 3805 3808 set(handles.slider2,'Visible','off') … … 4439 4442 Data_uvmat=get(handles.uvmat,'UserData'); 4440 4443 evalin('base','global Data_uvmat')%make CurData global in the workspace 4441 display( ['current field :'])4444 display('current field :') 4442 4445 evalin('base','Data_uvmat') %display CurData in the workspace 4443 4446 commandwindow; %brings the Matlab command window to the front … … 4670 4673 npx=size(UvData.Field.A,2); 4671 4674 npy=size(UvData.Field.A,1); 4672 xi= [0.5:npx-0.5];4673 yi= [0.5:npy-0.5];4675 xi=0.5:npx-0.5; 4676 yi=0.5:npy-0.5; 4674 4677 [Xi,Yi]=meshgrid(xi,yi); 4675 4678 if isfield(UvData,'Object')
Note: See TracChangeset
for help on using the changeset viewer.