Changeset 1103 for trunk/src/plot_field.m

Timestamp:

May 8, 2021, 11:01:15 PM (3 years ago)

Author:

sommeria

Message:

colormap options introduced

File:

• trunk/src/plot_field.m (modified) (11 diffs)

trunk/src/plot_field.m

@@ -865 +865 @@
     
     %set for grey scale setting
+    ColorMap='default';
     if isfield(PlotParam.Scalar,'CheckBW') && ~isempty(PlotParam.Scalar.CheckBW)
-        BW=PlotParam.Scalar.CheckBW; %BW=0 color imposed, else gray scale imposed.
-    else % BW imposed automatically chosen
-        BW=(siz==2) && (isa(A,'uint8')|| isa(A,'uint16'));% non color images represented in gray scale by default
-        PlotParamOut.Scalar.CheckBW=BW;
-    end
-    
+        ColorMap=PlotParam.Scalar.CheckBW; %BW=0 color imposed, else gray scale imposed.
+    elseif ((siz==2) && (isa(A,'uint8')|| isa(A,'uint16')))% non color images represented in gray scale by default
+        ColorMap='grayscale';
+    end
+    PlotParamOut.Scalar.CheckBW=ColorMap;
     % determine the plot option 'image' or 'contours'
     CheckContour=0; %default
@@ -915 +915 @@
     PlotParamOut.Scalar.Npx=size(A,2);
     PlotParamOut.Scalar.Npy=size(A,1);
-    %     if siz==2
+    
     % case of contour plot
     if CheckContour
@@ -946 +946 @@
         y_cont=Coord_y(1):-sizpy:Coord_y(end); % pixel x coordinates for image display
         
-        %axes(haxes)% set the input axes handle as current axis
-        
-        % colormap(map);
         tag_axes=get(haxes,'Tag');% axes tag
         Opacity=1;
@@ -972 +969 @@
         %determine the color scale and map
         caxis([abscontmin abscontmax])
-        if BW
+        if strcmp(ColorMap,'grayscale')
             vec=linspace(0,1,(abscontmax-abscontmin)/interval);%define a greyscale colormap with steps interval
             map=[vec' vec' vec'];
             colormap(map);
-        else
-            colormap('jet'); % default matlab colormap ('jet')
+        elseif strcmp(ColorMap,'BuYlRd')
+            hh=load('BuYlRd.mat');
+            colormap(hh.BuYlRd);
+        else
+            colormap(ColorMap); 
         end
     else %usual images (no contour)
         % set  colormap for  image display
-        if BW
+        if strcmp(ColorMap,'grayscale')
             vec=linspace(0,1,255);%define a linear greyscale colormap
             map=[vec' vec' vec'];
@@ -988 +988 @@
                 A=uint16(sum(A,3));%sum the three color components for color images displayed with BW option
             end
+        elseif strcmp(ColorMap,'BuYlRd')
+            hh=load('BuYlRd.mat');
+            colormap(hh.BuYlRd);
         else
             if siz==3 && CheckFixScalar % true color images rescaled by MaxA
                   A=uint8(255*double(A)/double(MaxA));
             end
-            colormap('default'); % standard false colors for div, vort , scalar fields
+            colormap(ColorMap); % standard false colors for div, vort , scalar fields
         end
         
@@ -1123 +1126 @@
 if test_vec
    %vector scale representation
-    if size(vec_U,1)==numel(vec_Y) && size(vec_U,2)==numel(vec_X); % x, y  coordinate variables
+    if size(vec_U,1)==numel(vec_Y) && size(vec_U,2)==numel(vec_X) % x, y  coordinate variables
         [vec_X,vec_Y]=meshgrid(vec_X,vec_Y);
     end   
@@ -1201 +1204 @@
     
     % take flags into account: add flag colors to the list of colors
-    sizlist=size(colorlist);
-    nbcolor=sizlist(1);
+    nbcolor=size(colorlist,1);
     if test_black 
        nbcolor=nbcolor+1;
@@ -1223 +1225 @@
     %plot vectors:
     quiresetn(haxes,vec_X,vec_Y,vec_U,vec_V,scale,colorlist,col_vec);   
-
 else
     hvec=findobj(haxes,'Tag','vel');
@@ -1231 +1232 @@
     PlotParamOut=rmfield(PlotParamOut,'Vectors');
 end
-% nbvar=0;
 
 %store the coordinate extrema occupied by the field
@@ -1297 +1297 @@
 alpha=0.3 ;%length arrow
 rot=alpha*[cos(theta) -sin(theta); sin(theta) cos(theta)]';
+
 %find the existing lines
 h=findobj(haxes,'Tag','vel');% search existing lines in the current axes
 sizh=size(h);
-%set(h,'EraseMode','xor');
 set(haxes,'NextPlot','replacechildren');
 
-%drawnow
 %create lines (if no lines) or modify them
 if ~isequal(size(col_vec),size(x))
     col_vec=ones(size(x));% case of error in col_vec input
 end
-sizlist=size(colorlist);
-ncolor=sizlist(1);
-
-for icolor=1:ncolor
+nbcolor=size(colorlist,1);
+
+for icolor=1:nbcolor
     %determine the line positions for each color icolor
     ind=find(col_vec==icolor);
@@ -1363 +1361 @@
     end
 end
-if sizh(1) > 2*ncolor
-    for icolor=ncolor+1 : sizh(1)/2%delete additional objects
+if sizh(1) > 2*nbcolor
+    for icolor=nbcolor+1 : sizh(1)/2 %delete additional objects
         delete(h(2*icolor-1))
         delete(h(2*icolor))