Changeset 349 in ETALON

Timestamp:

Mar 2, 2016, 9:44:51 AM (9 years ago)

Author:

malovyts

Message:

Added colormap option in the scatter 2d plot, and removed repetition frequency

File:

• SPESO/ana2015/data_scatter2d.m (modified) (2 diffs)

SPESO/ana2015/data_scatter2d.m

@@ -1 @@
-function sch = data_scatter2d(var1,var2,repeat_frequency)
+function sch = data_scatter2d(var1,var2,magnitude_var)
 
         global allData;
+        global viridis;
+
+        colormap_viridis;
+        colormap(viridis);
 
         [valsused,vals]=data_get_vals(var1,var2);
 
-        if nargin>2
-                vals(1,:)=mod(vals(1,:),repeat_frequency);
-        end
-
         if (isempty(vals(1,:))||isempty(vals(2,:)))
-                  warning_sp([ ' Warning one of the two variable is empty: ' length(vals(1,:)) ' / '  length(vals(2,:))  ])
-                else
+                  warning_sp([ ' Warning one of the two variable is empty: '...
+                                 length(vals(1,:)) ' / '  length(vals(2,:))  ])
+        elseif nargin==2 % scatter
                 hold on
                 if (valsused(4,1))
@@ -21 +22 @@
                 end % there is a var in pos 4
 
-                title({ [ get_var_name(var1) ' vs ' get_var_name(var2) ];allData.selection})
+                title({ [ get_var_name(var1) ' vs ' get_var_name(var2) ];...
+                          allData.selection})
+                xlabel(get_var_name(var1))
+                ylabel(get_var_name(var2))
+        elseif nargin==3 % scatter colormap
+                [valsused, magnitude]=data_get_vals(magnitude_var);
+                magnitude=magnitude(1,:);
+
+                % assigning colors to the values
+                min_magnitude=min(magnitude);
+                color_coeff=(max(magnitude)-min(magnitude));
+                colors=round((magnitude-min_magnitude)./color_coeff.*255)+1;
+                color_to_magn=@(color_in)((color_in-1)/255*color_coeff+min_magnitude);
+
+                % rewriting colorbar ticks
+                tick5=sprintf('%4.2e',(color_to_magn(250)));
+                tick4=sprintf('%4.2e',(color_to_magn(200)));
+                tick3=sprintf('%4.2e',(color_to_magn(150)));
+                tick2=sprintf('%4.2e',(color_to_magn(100)));
+                tick1=sprintf('%4.2e',(color_to_magn(50)));
+                %tick6=sprintf('%4.2e',magnitude(find(colors==0)));
+
+                sch = scatter(vals(1,:),vals(2,:),[],colors, 's','filled');
+
+                cbh = colorbar('v');%, 'YTickLabel',{tick1, tick2, tick3, tick4, tick5, tick6})
+                %set(cbh,'YTick',[0:50:250])
+                set(cbh,'YTickLabel',{tick1, tick2, tick3, tick4, tick5})
+                title({ [ get_var_name(var1) ' vs ' get_var_name(var2) ];...
+                          get_var_name(magnitude_var);...
+                          allData.selection})
                 xlabel(get_var_name(var1))
                 ylabel(get_var_name(var2))