Changeset 89 for trunk/src/uvmat.m

Timestamp:

May 20, 2010, 1:48:25 AM (14 years ago)

Author:

sommeria

Message:

many bug corrections and cleaning. Activation of the BW option in uvmat. Improvement of the interaction of get_field with uvmat.

File:

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

trunk/src/uvmat.m

@@ -165 +165 @@
 % Object{iobj}.plotaxes and nbobj representation graphs  Object{iobj}.plothandles(:) (where nbobj is the
 % nbre of current objects opened in uvmat. Note that Object{iobj}.plothandles(iobj)=[] : an object is not represented in its own projection field;
-%-------------------------------------------------------------------
-%-------------------------------------------------------------------
+
+%------------------------------------------------------------------------
+%------------------------------------------------------------------------
 %  I - MAIN FUNCTION UVMAT (DO NOT MODIFY)
-%-------------------------------------------------------------------
-%-------------------------------------------------------------------
+%------------------------------------------------------------------------
+%------------------------------------------------------------------------
 function varargout = uvmat(varargin)
 
@@ -191 +192 @@
 % End initialization code - DO NOT EDIT
 
-%-------------------------------------------------------------------
+%------------------------------------------------------------------------
 % --- Executes just before uvmat is made visible.
 function uvmat_OpeningFcn(hObject, eventdata, handles, input )
-%-------------------------------------------------------------------
+%------------------------------------------------------------------------
 global nb_builtin
 
@@ -221 +222 @@
 set(hObject,'DeleteFcn',{@closefcn})%
 
-%TRANSFORM menu: loads the information stored in prefdir to initiate the browser and the list of functions
+%refresh projection plane
+UvData.Object{1}.Style='plane';%main plotting plane
+UvData.Object{1}.ProjMode='projection';%main plotting plane
+if ~isfield(UvData.Object{1},'plotaxes')
+    UvData.Object{1}.plotaxes=handles.axes3;%default plotting axis
+    set(handles.list_object_1,'Value',1);
+    set(handles.list_object_1,'String',{'1-PLANE'});
+end
+
+%TRANSFORM menu: builtin fcts
 menu_str={'';'phys';'px';'phys_polar'};
 nb_builtin=numel(menu_str); %number of functions
@@ -238 +248 @@
 end
 rmpath(fullfile(path_uvmat,'transform_field'))
-
-%refresh projection plane
-UvData.Object{1}.Style='plane';%main plotting plane
-UvData.Object{1}.ProjMode='projection';%main plotting plane
-if ~isfield(UvData.Object{1},'plotaxes')
-    UvData.Object{1}.plotaxes=handles.axes3;%default plotting axis
-    set(handles.list_object_1,'Value',1);
-    set(handles.list_object_1,'String',{'1-PLANE'});
-end
 
 %load the list of previously browsed files in menus Open and Open_1
@@ -279 +280 @@
              if exist(file,'file')
                 h_func=str2func(file);
-                testexist=[testexist 1]; %#ok<AGROW>
+                testexist=[testexist 1]; 
              else
                 h_func=[];
-                testexist=[testexist 0];  %#ok<AGROW>
+                testexist=[testexist 0]; 
              end
-             fct_handle=[fct_handle; {h_func}];%#ok<AGROW> %concatene the list of paths
+             fct_handle=[fct_handle; {h_func}]; %concatene the list of paths
              rmpath(path)
-             menu_str=[menu_str; {file}]; %#ok<AGROW>
+             menu_str=[menu_str; {file}]; 
          end
       end
@@ -344 +345 @@
         set(handles.Fields,'Value',1)
         set(handles.Fields,'String',{'get_field...'})    
-       % set(handles.Fields,'UserData',Field)
         testinputfield=1;
         
@@ -577 +577 @@
             [ff,rr,FileBase_1,ii,FileExt_1,SubDir_1]=read_file_boxes_1(handles);
             NomType_1=get(handles.FileIndex_1,'UserData');     
-            FileName_1=name_generator(FileBase_1,stra2num(i1),stra2num(i2),FileExt_1,NomType_1,1,stra2num(str_a),stra2num(str_b),SubDir_1);
+            FileName_1=name_generator(FileBase_1,str2double(i1),str2double(i2),FileExt_1,NomType_1,1,stra2num(str_a),stra2num(str_b),SubDir_1);
             if exist(FileName_1,'file')
-                FileIndex_1=name_generator('',stra2num(i1),stra2num(i2),'',NomType_1,1,stra2num(str_a),stra2num(str_b),'');
+                FileIndex_1=name_generator('',str2double(i1),str2double(i2),'',NomType_1,1,stra2num(str_a),stra2num(str_b),'');
                 set(handles.FileIndex_1,'String',FileIndex_1)
             else
@@ -608 +608 @@
 end
 
-%-------------------------------------------------------------------
+%------------------------------------------------------------------------
 function RootPath_Callback(hObject,eventdata,handles)
+%------------------------------------------------------------------------
 update_rootinfo(hObject,eventdata,handles);
 
-%-------------------------------------------------------------------
+%------------------------------------------------------------------------
 %-- called by action in RootFile edit box
-%-------------------------------------------------------------------
 function SubDir_Callback(hObject, eventdata, handles)
+%------------------------------------------------------------------------
 %refresh the menu of input fields
 Fields_Callback(hObject, eventdata, handles);
@@ -621 +622 @@
 run0_Callback(hObject, eventdata, handles); 
 
-
-%-------------------------------------------------------------------
-%-- called by action in RootFile edit box
-%-------------------------------------------------------------------
+%------------------------------------------------------------------------
+% --- called by action in RootFile edit box
 function RootFile_Callback(hObject, eventdata, handles)
+%------------------------------------------------------------------------
 update_rootinfo(hObject,eventdata,handles)
 
-%-------------------------------------------------------------------
-%-- called by action in FileIndex edit box
-%-------------------------------------------------------------------
+%------------------------------------------------------------------------
+% --- called by action in FileIndex edit box
 function FileIndex_Callback(hObject, eventdata, handles)
+%------------------------------------------------------------------------
 FileIndices=get(handles.FileIndex,'String');
 if isempty(str2num(FileIndices))
@@ -647 +647 @@
 run0_Callback(hObject, eventdata, handles)
 
-%-------------------------------------------------------------------
-%-- called by action in FileIndex_1 edit box
-%-------------------------------------------------------------------
+%------------------------------------------------------------------------
+% --- called by action in FileIndex_1 edit box
 function FileIndex_1_Callback(hObject, eventdata, handles)
-% FileIndices=get(handles.FileIndex_1,'String');
-% if isempty(str2num(FileIndices))
-%     [pp,ff,str1,str2,str_a,str_b]=name2display(FileIndices)
-% else
-%     str1=FileIndices;
-%     str2='';
-%     str_a='';
-%     str_b='';
-% end
-% set(handles.i1,'String',str1);
-% set(handles.i2,'String',str2);
-% set(handles.j1,'String',str_a);
-% set(handles.j2,'String',str_b);
+%------------------------------------------------------------------------
 run0_Callback(hObject, eventdata, handles)
 
-%-------------------------------------------------------------------
-% -- update information about a new field series (indices to scan, timing, calibration from an xml file, then refresh current plots
-%-------------------------------------------------------------------
+%------------------------------------------------------------------------
+% --- update information about a new field series (indices to scan, timing, calibration from an xml file, then refresh current plots
 function update_rootinfo(hObject,eventdata,handles)
-
+%------------------------------------------------------------------------
 set(handles.RootPath,'BackgroundColor',[1 1 0])
 drawnow
@@ -678 +664 @@
 UvData.TestInputFile=1;
 set(handles.fix_pair,'Value',1) % activate by default the comp_input '-'input window
-%FileIndex_Callback(hObject, eventdata, handles)% update field counters
 
 [FileName,RootPath,FileBase,FileIndices,FileExt,SubDir]=read_file_boxes(handles);
@@ -693 +678 @@
 TimeUnit=[];%default
 testima=0; %test for image input
+imainfo=[];
+ColorType='falsecolor'; %default
+hhh='';
 if isequal(lower(FileExt),'.avi') %.avi file
     testima=1;
-    info=aviinfo([FileBase FileIndices FileExt]);
-    nbfield=info.NumFrames;
+    imainfo=aviinfo([FileBase FileIndices FileExt]);
+    nbfield=imainfo.NumFrames;
     nburst=1;
-    set(handles.Dt_txt,'String',['Dt=' num2str(1000/info.FramesPerSecond) 'ms']);%display the elementary time interval in millisec
-    XmlData.Time=(0:1/info.FramesPerSecond:(info.NumFrames-1)/info.FramesPerSecond)';
+    set(handles.Dt_txt,'String',['Dt=' num2str(1000/imainfo.FramesPerSecond) 'ms']);%display the elementary time interval in millisec
+    XmlData.Time=(0:1/imainfo.FramesPerSecond:(imainfo.NumFrames-1)/imainfo.FramesPerSecond)';
     TimeUnit='s';
-    set(handles.npx,'String',num2str(info.Width));%fills nbre of pixels x box
-    set(handles.npy,'String',num2str(info.Height));%fills nbre of pixels y box
     hhh=which('mmreader');
-    if ~isequal(hhh,'')&& mmreader.isPlatformSupported()% if the function is found (recent version of matlab)
-        UvData.MovieObject=mmreader([FileBase FileIndices FileExt]);
-    end
+    ColorType=imainfo.ImageType;%='truecolor' for color images
 elseif ~isempty(imformats(FileExt(2:end))) || isequal(FileExt,'.vol')%&& isequal(NomType,'*')% multi-frame image
     testima=1;
@@ -715 +699 @@
         imainfo=imfinfo([FileBase FileIndices FileExt]);
     end
+    ColorType=imainfo.ColorType;%='truecolor' for color images
     if length(imainfo) >1 %case of image with multiple frames
         nbfield=length(imainfo);
         nburst=1;
     end 
-    if isfield(UvData,'MovieObject')
-        UvData=rmfield(UvData,'MovieObject');
-    end
-else 
-    %set(handles.last_i,'String','');%fills last_field box
+end
+if ~strcmp(hhh,'')&& mmreader.isPlatformSupported()% if the function is found (recent version of matlab)
+    UvData.MovieObject=mmreader([FileBase FileIndices FileExt]);
+elseif isfield(UvData,'MovieObject')
+    UvData=rmfield(UvData,'MovieObject');
+end
+if isfield(imainfo,'Width') && isfield(imainfo,'Height')
+    set(handles.npx,'String',num2str(imainfo.Width));%fills nbre of pixels x box
+    set(handles.npy,'String',num2str(imainfo.Height));%fills nbre of pixels x box
+else
     set(handles.npx,'String','');%fills nbre of pixels x box
-    set(handles.npy,'String','');%fills nbre of pixels y box
-    if isfield(UvData,'MovieObject')
-        UvData=rmfield(UvData,'MovieObject');
-    end
-end
+    set(handles.npy,'String','');%fills nbre of pixels x box
+end
+set(handles.BW,'Value',strcmp(ColorType,'grayscale'))% select handles.BW if grayscale image
 
 % read parameters (time, geometric calibration..) from a documentation file (.xml advised)
@@ -821 +809 @@
         set(handles.transform_fct,'Value',1); %  no transform by default
     else
-        if (isfield(GeometryCalib,'R')& ~isequal(GeometryCalib.R(2,1),0) & ~isequal(GeometryCalib.R(1,2),0)) |...
-            (isfield(GeometryCalib,'kappa1')& ~isequal(GeometryCalib.kappa1,0))
+        if (isfield(GeometryCalib,'R')&& ~isequal(GeometryCalib.R(2,1),0) && ~isequal(GeometryCalib.R(1,2),0)) |...
+            (isfield(GeometryCalib,'kappa1')&& ~isequal(GeometryCalib.kappa1,0))
             set(handles.pxcm,'String','var')
             set(handles.pycm,'String','var')
@@ -835 +823 @@
         end
         if isfield(GeometryCalib,'SliceCoord')
+            
            siz=size(GeometryCalib.SliceCoord);
            if siz(1)>1
@@ -847 +836 @@
            end
            slices_Callback(hObject, eventdata, handles)
-    %        Coord=UvData.XmlData.GeometryCalib.SliceCoord;
-    %        ZIndex=num_i1-NbSlice*(floor((num_i1-1)/NbSlice));
-        end
-
-            
+        end           
     end
 end
 
 %update the data attached to the uvmat interface
-% set(handles.nb_slice,'String',num2str(NbSlice)) 
 if ~isempty(TimeUnit)
     set(handles.time_txt,'String',['time (' TimeUnit ')'])
 end
-%set(handles.DtUnit,'String',['10^(-3)' TimeUnit])
 UvData.TimeUnit=TimeUnit;
 UvData.XmlData=XmlData;
@@ -986 +969 @@
 end
 sizf=size(fileinput_1);
-if (~ischar(fileinput_1)|~isequal(sizf(1),1)),return;end
+if (~ischar(fileinput_1)||~isequal(sizf(1),1)),return;end
 
 % refresh the current displayed field
@@ -1197 +1180 @@
 end
 TimeUnit_1=[];
-% testima=0; %test for image input
+hhh='';%default, test for  movie reading with mmreader
+imainfo=[];
 if isequal(lower(FileExt),'.avi') %.avi file
-%     testima=1;
-    info=aviinfo([FileBase FileIndices FileExt]);
-    nbfield_1=info.NumFrames;
+    imainfo=aviinfo([FileBase FileIndices FileExt]);
+    nbfield_1=imainfo.NumFrames;
     nburst_1=1;
-    %set(handles.last_i,'String',num2str(info.NumFrames));%fills last_field box
     set(handles.Dt_txt,'String',['Dt=' num2str(1000/info.FramesPerSecond) 'ms']);%display the elementary time interval in millisec
-%     set(handles.dt,'UserData',info.FramesPerSecond);
-    time=[0:1/info.FramesPerSecond:(info.NumFrames-1)/info.FramesPerSecond]';
-    npx=get(handles.npx,'String');
-    npy=get(handles.npy,'String');
-    if isequal(npx,'') || isequal(npy,'')
-        set(handles.npx,'String',num2str(info.Width));%fills nbre of pixels x box
-        set(handles.npy,'String',num2str(info.Height));%fills nbre of pixels y box
-    end
-elseif ~isempty(imformats(FileExt(2:end))) %&& isequal(NomType,'*')% multi-frame image
+    time=[0:1/imainfo.FramesPerSecond:(imainfo.NumFrames-1)/imainfo.FramesPerSecond]';
+    ColorType=imainfo.ImageType;%='truecolor' for color images
+    hhh=which('mmreader');
+elseif ~isempty(imformats(FileExt(2:end)))|| isequal(FileExt,'.vol')
     testima=1;
-    imainfo=imfinfo([FileBase FileIndices FileExt]);  
+    if ~isequal(SubDir,'')
+        RootFile=get(handles.RootFile,'String');
+        imainfo=imfinfo([fullfile(RootPath,SubDir,RootFile) FileIndices FileExt]);
+    else
+        imainfo=imfinfo([FileBase FileIndices FileExt]);
+    end
+    ColorType=imainfo.ColorType;%='truecolor' for color images
     if length(imainfo) >1 %case of image with multiple frames
         nbfield_1=length(imainfo);
         nburst_1=1;
-        %set(handles.last_i,'String',num2str(length(imainfo)));%fills last_field box
-    end       
-else 
-    %set(handles.last_i,'String','');%fills last_field box
-    set(handles.npx,'String','');%fills nbre of pixels x box
-    set(handles.npy,'String','');%fills nbre of pixels y box 
+    end 
+end
+if ~strcmp(hhh,'')&& mmreader.isPlatformSupported()% if the function is found (recent version of matlab)
+    UvData.MovieObject_1=mmreader([FileBase FileIndices FileExt]);
+elseif isfield(UvData,'MovieObject_1')
+    UvData=rmfield(UvData,'MovieObject_1');
+end
+if strcmp(get(handles.npx,'String'),'') || strcmp(get(handles.npy,'String'),'')%update npx and npy if it is not already filled by the first input field
+    if  isfield(imainfo,'Width') && isfield(imainfo,'Height')
+        set(handles.npx,'String',num2str(imainfo.Width));%fills nbre of pixels x box
+        set(handles.npy,'String',num2str(imainfo.Height));%fills nbre of pixels x box
+    else
+        set(handles.npx,'String','');%fills nbre of pixels x box
+        set(handles.npy,'String','');%fills nbre of pixels x box
+    end
+    set(handles.BW,'Value',strcmp(ColorType,'grayscale'))% select handles.BW if grayscale image
 end
 
@@ -1849 +1842 @@
 end
 
-num_i1=stra2num(get(handles.i1,'String'));
+num_i1=str2double(get(handles.i1,'String'));
 num_j1=stra2num(get(handles.j1,'String'));
-num_i2=stra2num(get(handles.i2,'String'));
+num_i2=str2double(get(handles.i2,'String'));
 num_j2=stra2num(get(handles.j2,'String'));
-if isempty(num_j2)
+if isnan(num_j2)
     if isempty(num_i2)   
         msgbox_uvmat('ERROR', 'a second image index i2 or j2 is needed to show the pair as a movie')
@@ -1862 +1855 @@
     end
 end
-if isempty(num_i2)
+if isnan(num_i2)
     num_i2=num_i1;%repeat the index i1 by default
 end
@@ -2238 +2231 @@
     if isequal(FileType,'netcdf')  %read the first nc field
         if isequal(FieldName,'get_field...')% read the field names on the interface get_field.
-            %VelType=get(handles.Fields,'UserData');
-            hget_field=findobj(allchild(0),'Name','get_field')%find the get_field... GUI
+            hget_field=findobj(allchild(0),'Name','uvmat_field');%find the get_field... GUI
             if isempty(hget_field)
-                hget_field= get_field(filename);%open the get_field GUI    
-            end
-            if numel(hget_field)>1%case of several get_field GUI
-                hget_field=hget_field(1);
+                hget_field= get_field(filename);%open the get_field GUI   
+                set(hget_field,'name','uvmat_field')
+            elseif UvData.NewSeries
+                delete(hget_field)
+                hget_field= get_field(filename);%open the get_field GUI 
+                set(hget_field,'name','uvmat_field')
             end
             hhget_field=guidata(hget_field);
-            set(hhget_field.inputfile,'String',filename)% update the list of input fields in get_field
-            set(hhget_field.ACTION,'Value',1)% PLOT option selected
-            set(hhget_field.list_fig,'Value',2)% plotting axes =uvmat selected
-            [Field{1},errormsg]=read_get_field(hget_field); %read the names of the variables to plot in the get_field GUI
-            if ~isempty(errormsg)
-                errormsg=['error in uvmat/run0_Callback/read_get_field: ' errormsg];
-                return
-            end
+            funct_list=get(hhget_field.ACTION,'UserData');
+            funct_index=get(hhget_field.ACTION,'Value');
+            funct=funct_list{funct_index};
+            Field{1}=funct(hget_field); %read the names of the variables to plot in the get_field GUI
             CivStage=0;
-            VelType_out=[];         
+            VelType_out=[];    
         else
             [Field{1},VelType_out]=read_civxdata(filename,InputField,VelType);
@@ -2269 +2259 @@
     if ~isempty(filename_1) && ~test_keepdata_1 && isequal(FileType_1,'netcdf') %read the second file
         if isequal(FieldName_1,'get_field...')% read the field names on the interface get_field.
-            hget_field=findobj(allchild(0),'Name','get_field_1');%find the get_field... GUI
-             if isempty(hget_field)
-                 hget_field= get_field(filename_1);%open the get_field GUI
-                 set(hget_field,'name','get_field_1')
-%                 enable_transform(handles,'off')% no field transform (possible transform in the GUI get_field)
-             end
-            hhget_field=guidata(hget_field);%handles of GUI elements in get_field
-            SubField=get_field('read_var_names',hObject,eventdata,hhget_field); %read the names of the variables to plot in the get_field GUI 
-            [Field{2},var_detect]=nc2struct(filename_1,SubField.ListVarName); %read the corresponding input data                
-            Field{2}.VarAttribute=SubField.VarAttribute;
-            %update the display on get_field
-            set(hhget_field.inputfile,'String',filename_1)
-            set(hhget_field.variables,'Value',1)
-            Tabchar={''};%default
-            Tabcell=[];
-            if isfield(Field{2},'ListGlobalAttribute')& ~isempty(Field{2}.ListGlobalAttribute)
-                for iline=1:length(Field{2}.ListGlobalAttribute)
-                    Tabcell{iline,1}=Field{2}.ListGlobalAttribute{iline};
-                    if isfield(Field{2}, Field{2}.ListGlobalAttribute{iline})
-                        eval(['val=Field{2}.' Field{2}.ListGlobalAttribute{iline} ';'])
-                        if ischar(val);
-                            Tabcell{iline,2}=val;
-                        else
-                            Tabcell{iline,2}=num2str(val);
-                        end
-                    end
-                end
-                if ~isempty(Tabcell)
-                    Tabchar=cell2tab(Tabcell,'=');
-                    Tabchar=[{''};Tabchar];
-                end
+            hget_field_1=findobj(allchild(0),'Name','uvmat_field_1');%find the get_field... GUI
+            if isempty(hget_field_1)
+                 hget_field_1= get_field(filename_1);%open the get_field GUI
+                 set(hget_field_1,'name','uvmat_field_1')
             end
-            set(hhget_field.attributes,'String',Tabchar);%update list of global attributes in get_field 
+            hhget_field_1=guidata(hget_field_1);%handles of GUI elements in get_field
+            funct_list=get(hhget_field_1.ACTION,'UserData');
+            funct_index=get(hhget_field_1.ACTION,'Value');
+            funct=funct_list{funct_index};
+            Field{2}=funct(hget_field_1); %read the names of the variables to plot in the get_field GUI
         else
             [Field{2},VelType_out_1]=read_civxdata(filename_1,InputField_1,VelType_1);
@@ -2379 +2346 @@
        if siz(1)>1
            NbSlice=siz(1);
-%            set(handles.slices,'Visible','on')
-%            set(handles.slices,'Value',1)
        else
            NbSlice=1;
        end
-       %set(handles.nb_slice,'String',num2str(NbSlice))
-%        slices_Callback(handles.uvmat, [], handles)
 end
 
@@ -2510 +2473 @@
 end
 %case of structured coordinates
-if isfield(UvData.Field,'AX') & isfield(UvData.Field,'AY')& isfield(UvData.Field,'A')
+if isfield(UvData.Field,'AX') && isfield(UvData.Field,'AY')&& isfield(UvData.Field,'A')
     UvData.XMax=max(UvData.Field.AX);
     UvData.XMin=min(UvData.Field.AX);
@@ -2518 +2481 @@
     UvData.Mesh=sqrt((UvData.XMax-UvData.XMin)*(UvData.YMax-UvData.YMin)/((np_A(1)-1) * (np_A(2)-1))) ; 
 end
-if  isempty(coord_x)&~isempty(CellVarIndex)
+if  isempty(coord_x) && ~isempty(CellVarIndex)
     VarIndex=CellVarIndex{imax}; % list of variable indices
     DimIndex=UvData.Field.VarDimIndex{VarIndex(1)}; %list of dim indices for the variable
@@ -2640 +2603 @@
         %Projeter les champs sur l'objet:*
         ObjectData=proj_field(UvData.Field,UvData.Object{iobj},iobj);
-%         if imap==2
-%             UvData.ProjField_2=ObjectData;%store the projection field on uvmat: ***** WILL REPLACE THE FIELD SORED ON THE AXES: AxeData *****
-%         else
-%             UvData.ProjField_1=ObjectData;%store the projection field on view_field
-%         end
         %use of mask
         if isfield(ObjectData,'NbDim')&isequal(ObjectData.NbDim,2)
@@ -2778 +2736 @@
 testimedoc=0;
 if isfield(UvData,'XmlData') && isfield(UvData.XmlData,'Time')
-    if isempty(num_i2)
+    if isempty(num_i2)||isnan(num_i2)
         num_i2=num_i1;
     end
-    if isempty(num_j1)
+    if isempty(num_j1)||isnan(num_j1)
         num_j1=1;
     end
-    if isempty(num_j2)
+    if isempty(num_j2)||isnan(num_j2)
         num_j2=num_j1;
     end
@@ -2828 +2786 @@
     end
 end
-
+%update the input file name in the get_field GUI
+if isequal(FieldName,'get_field...')
+    set(hhget_field.inputfile,'String',filename)
+    Tabchar={''};%default
+    Tabcell=[];
+    if isfield(Field{1},'ListGlobalAttribute')& ~isempty(Field{1}.ListGlobalAttribute)
+        for iline=1:length(Field{1}.ListGlobalAttribute)
+            Tabcell{iline,1}=Field{1}.ListGlobalAttribute{iline};
+            if isfield(Field{1}, Field{1}.ListGlobalAttribute{iline})
+                eval(['val=Field{1}.' Field{1}.ListGlobalAttribute{iline} ';'])
+                if ischar(val);
+                    Tabcell{iline,2}=val;
+                else
+                    Tabcell{iline,2}=num2str(val);
+                end
+            end
+        end
+        if ~isempty(Tabcell)
+            Tabchar=cell2tab(Tabcell,'=');
+            Tabchar=[{''};Tabchar];
+        end
+    end
+    set(hhget_field.attributes,'String',Tabchar);%update list of global attributes in get_field 
+end
+if isequal(FieldName_1,'get_field...')
+    set(hhget_field_1.inputfile,'String',filename_1)
+    Tabchar={''};%default
+    Tabcell=[];
+    if isfield(Field{2},'ListGlobalAttribute')& ~isempty(Field{2}.ListGlobalAttribute)
+        for iline=1:length(Field{2}.ListGlobalAttribute)
+            Tabcell{iline,1}=Field{2}.ListGlobalAttribute{iline};
+            if isfield(Field{2}, Field{2}.ListGlobalAttribute{iline})
+                eval(['val=Field{2}.' Field{2}.ListGlobalAttribute{iline} ';'])
+                if ischar(val);
+                    Tabcell{iline,2}=val;
+                else
+                    Tabcell{iline,2}=num2str(val);
+                end
+            end
+        end
+        if ~isempty(Tabcell)
+            Tabchar=cell2tab(Tabcell,'=');
+            Tabchar=[{''};Tabchar];
+        end
+    end
+    set(hhget_field_1.attributes,'String',Tabchar);%update list of global attributes in get_field 
+end
 
 
@@ -2858 +2862 @@
     set(handles.auto_xy,'BackgroundColor',[1 1 0])
     cla(handles.axes3)
-    update_plot(handles)
+    update_plot(handles);
 else
     set(handles.auto_xy,'BackgroundColor',[0.7 0.7 0.7])
-    update_plot(handles)
+    update_plot(handles);
 %     axis(handles.axes3,'image')
 end
@@ -3039 +3043 @@
 end
 
-% %----------------------------------------------
-% %read the data displayed for the input rootfile windows (new)
-% %-------------------------------------------------
+%------------------------------------------------------------------------
+% --- read the data displayed for the input rootfile windows (new)
 function [FileName,RootPath,FileBase,FileIndices,FileExt,SubDir]=read_file_boxes(handles)
+%------------------------------------------------------------------------
 RootPath=get(handles.RootPath,'String');
 FileName=RootPath; %default
@@ -3114 +3118 @@
      set_veltype_display(veltype_handles,0) % unvisible civ buttons
      filename=read_file_boxes(handles);
-     hget_field=findobj(allchild(0),'name','get_field');
-     if ~isempty(hget_field)
-         delete(hget_field)
+     hget_field=findobj(allchild(0),'name','uvmat_field');
+     if isempty(hget_field)
+         hget_field=get_field(filename);
+         set(hget_field,'Name','uvmat_field')        
      end
-     get_field(filename)
     return %no action
 end
@@ -3849 +3853 @@
 % set(handles.slider1,'Value',col) 
 set_vec_col_bar(handles)
-update_plot(handles)
+update_plot(handles);
 
 %----------------------------------------------------------------
@@ -3859 +3863 @@
 % slider2_Callback(hObject, eventdata, handles)
 set_vec_col_bar(handles)
-update_plot(handles)
+update_plot(handles);
 %------------------------------------------------------------
 %update the slider values after displaying vectors
@@ -3885 +3889 @@
 set_vec_col_bar(handles)
 
-% %--------------------------------------------
-% %update the display of color code for vectors
-% %--------------------------------------------
-% function set_vec_col_bar(handles)
-% %get the image of the color display button 'vec_col_bar' in pixels
-% uni=get(handles.vec_col_bar,'Unit');
-% set(handles.vec_col_bar,'Unit','pixel')
-% pos_vert=get(handles.vec_col_bar,'Position');
-% set(handles.vec_col_bar,'Unit','Normalized')
-% width=ceil(pos_vert(3));
-% height=ceil(pos_vert(4));
-% %get slider indications
-% colcode.min=get(handles.slider1,'Min');
-% colcode.max=get(handles.slider1,'Max');
-% colcode.colcode1=get(handles.slider1,'Value');
-% colcode.colcode2=get(handles.slider2,'Value');
-% colcode.option=get(handles.vec_col_bar,'Value');
-% colcode.auto=1;
-% list_code=get(handles.col_vec,'String');% list menu fields
-% index_code=get(handles.col_vec,'Value');% selected string index
-% colcode.CName= list_code{index_code(1)}; % selected field used for vector color
-% vec_C=colcode.min+(colcode.max-colcode.min)*[0.5:width-0.5]/width;%sample of vec_C values from min to max
-% [colorlist,col_vec]=set_col_vec(colcode,vec_C);
-% oneheight=ones(1,height);
-% A1=colorlist(col_vec,1)*oneheight;
-% A2=colorlist(col_vec,2)*oneheight;
-% A3=colorlist(col_vec,3)*oneheight;
-% A(:,:,1)=A1';
-% A(:,:,2)=A2';
-% A(:,:,3)=A3';
-% set(handles.vec_col_bar,'Cdata',A)
-
-%--------------------------------------------------------
-% --- Executes on button press in cal.
-function cal_Callback(hObject, eventdata, handles)
-
-
-
 %-------------------------------------------------------------
 % --- Executes on selection change in transform_fct.
@@ -3940 +3906 @@
     dlg_title = 'user defined transform';
     num_lines= 1;
-    [FileName, PathName, filterindex] = uigetfile( ...
+    [FileName, PathName] = uigetfile( ...
        {'*.m', ' (*.m)';
         '*.m',  '.m files '; ...
@@ -4015 +3981 @@
 list_object_2_Callback(hObject, eventdata, handles)
 
-
-
 %delete mask if it is displayed 
 if isequal(get(handles.mask_test,'Value'),1)%if the mask option is on
@@ -4115 +4079 @@
 set(handles.AutoScal,'Value',1) %suppress auto mode
 set(handles.AutoScal,'BackgroundColor',[1 1 0])
-update_plot(handles)
+update_plot(handles);
 
 %-----------------------------------------------------------------
@@ -4122 +4086 @@
 set(handles.AutoScal,'Value',1) %suppress auto mode
 set(handles.AutoScal,'BackgroundColor',[1 1 0])
-update_plot(handles)
+update_plot(handles);
 
 %-----------------------------------------------
@@ -4133 +4097 @@
     set(handles.AutoScal,'BackgroundColor',[0.7 0.7 0.7])
     update_plot(handles);
-%     set(handles.MinA,'String',num2str(ScalOut.MinA,3))
-%     set(handles.MaxA,'String',num2str(ScalOut.MaxA,3))
 end
 
@@ -4140 +4102 @@
 function BW_Callback(hObject, eventdata, handles)
 %-------------------------------------------------------------------
-update_plot(handles)
+update_plot(handles);
 
 %-------------------------------------------------------------------
@@ -4153 +4115 @@
     set(handles.IncrA,'Visible','off')
 end
-update_plot(handles)
+update_plot(handles);
 
 %-------------------------------------------------------------------
 function IncrA_Callback(hObject, eventdata, handles)
 %-------------------------------------------------------------------
-update_plot(handles)
+update_plot(handles);
 
 %-------------------------------------------------------------------
 function HideWarning_Callback(hObject, eventdata, handles)
 %-------------------------------------------------------------------
-update_plot(handles)
+update_plot(handles);
 
 %-------------------------------------------------------------------
 function HideFalse_Callback(hObject, eventdata, handles)
 %-------------------------------------------------------------------
-update_plot(handles)
+update_plot(handles);
 
 %-------------------------------------------------------------------
@@ -4175 +4137 @@
 set(handles.AutoVec,'Value',1);
 set(handles.AutoVec,'BackgroundColor',[1 1 0])
-update_plot(handles)
+update_plot(handles);
 
 %-------------------------------------------------------------------
@@ -4189 +4151 @@
 end
 
-%-------------------------------------------------------
+%------------------------------------------------------------------------
 % --- Executes on selection change in decimate4 (nb_vec/4).
-%-------------------------------------------------------
 function decimate4_Callback(hObject, eventdata, handles)
-update_plot(handles)
-
-
-%-------------------------------------------------------
+%------------------------------------------------------------------------
+update_plot(handles);
+
+%------------------------------------------------------------------------
 % --- Executes on selection change in color_code menu
-%-------------------------------------------------------
 function color_code_Callback(hObject, eventdata, handles)
+%------------------------------------------------------------------------
 set_vec_col_bar(handles)
 update_plot(handles);
 
-%-------------------------------------------------------
+%------------------------------------------------------------------------
 % --- Executes on button press in AutoVecColor.
-%-------------------------------------------------------
 function AutoVecColor_Callback(hObject, eventdata, handles)
+%------------------------------------------------------------------------
 test=get(handles.AutoVecColor,'Value');
 if test
@@ -4215 +4176 @@
     set(handles.AutoVecColor,'BackgroundColor',[0.7 0.7 0.7])
 end
-%set_vec_col_bar(handles)
-
-%-------------------------------------------------------
+
+%------------------------------------------------------------------------
 % --- Executes on selection change in max_vec.
-%-------------------------------------------------------
 function min_vec_Callback(hObject, eventdata, handles)
+%------------------------------------------------------------------------
 max_vec_Callback(hObject, eventdata, handles)
 
+%------------------------------------------------------------------------
 % --- Executes on selection change in max_vec.
 function max_vec_Callback(hObject, eventdata, handles)
+%------------------------------------------------------------------------
 set(handles.AutoVecColor,'Value',1)
 AutoVecColor_Callback(hObject, eventdata, handles)
@@ -4237 +4199 @@
 update_plot(handles);
 
-%-------------------------------------------------------------------
-%update the display of color code for vectors
+%------------------------------------------------------------------------
+% --- update the display of color code for vectors
 function set_vec_col_bar(handles)
-%-------------------------------------------------------------------
+%------------------------------------------------------------------------
 %get the image of the color display button 'vec_col_bar' in pixels
 set(handles.vec_col_bar,'Unit','pixel');
@@ -4283 +4245 @@
 AxeData=get(haxes,'UserData');
 PlotParam=read_plot_param(handles);
+PlotParam.Scalar
 [PlotType,PlotParamOut]= plot_field(AxeData,haxes,PlotParam,1);
 write_plot_param(handles,PlotParamOut); %update the auto plot parameters
@@ -4446 +4409 @@
 %------------------------------------------------------
 function MenuExportField_Callback(hObject, eventdata, handles)
-global CurData
-CurData=get(handles.uvmat,'UserData');
-evalin('base','global CurData')%make CurData global in the workspace
+global Data_uvmat
+Data_uvmat=get(handles.uvmat,'UserData');
+evalin('base','global Data_uvmat')%make CurData global in the workspace
 display(['current field :'])
-evalin('base','CurData') %display CurData in the workspace
+evalin('base','Data_uvmat') %display CurData in the workspace
 commandwindow; %brings the Matlab command window to the front
 
@@ -4976 +4939 @@
        '*.xml',  '.xml files '; ...
         '*.mat',  '.mat matlab files '}, ...
-        'Pick a file',get(handles.RootPath,'String'));
+        'Pick an xml Object file',get(handles.RootPath,'String'));
 fileinput=[PathName FileName];%complete file name 
 testblank=findstr(fileinput,' ');%look for blanks