Changeset 16053

Timestamp:

09/03/13 08:57:35 (12 years ago)

Author:

Mathieu Morlighem

Message:

NEW: automatically disseminate data, depending on cutoff value

File:

• issm/trunk-jpl/src/m/plot/plot_scatter.m (modified) (2 diffs)

issm/trunk-jpl/src/m/plot/plot_scatter.m

@@ -1 +1 @@
 function plot_scatter(x,y,level,varargin),
+%PLOT_SCATTER - scatter plot
+%
+%   Usage:
+%      plot_scatter(x,y,level,options);
+%
+%   Available options:
+%      'MarkerSize' : default is 3
+%      'Line'       : use line instead of circles
+%      'Cutoff'     : cut the line if the distance between 2 points is
+%                     greater than Cutoff (default is 1000)
 
-options=pairoptions(varargin{:}); 
+if nargin == 4,
+        options = varargin{1};
+else
+        options=pairoptions(varargin{:}); 
+end
 
 %check input
@@ -47 +61 @@
 hold on
 hp=[];
-for i=1:numcolors
-        pos=find(colorind==i);
-        hprime=plot3(x(pos),y(pos),ones(size(x(pos))),...
-                'o','MarkerSize',getfieldvalue(options,'MarkerSize',3),'MarkerEdgeColor',palette(i,:),...
-                'MarkerFaceColor',palette(i,:));
-        hp=[hp;hprime];
+if ~exist(options,'line'),
+        for i=1:numcolors
+                pos=find(colorind==i);
+                hprime=plot3(x(pos),y(pos),ones(size(x(pos))),...
+                        'o','MarkerSize',getfieldvalue(options,'MarkerSize',3),'MarkerEdgeColor',palette(i,:),...
+                        'MarkerFaceColor',palette(i,:));
+                hp=[hp;hprime];
+        end
+else
+        distances = sqrt( (x(1:end-1)-x(2:end)).^2 + (y(1:end-1)-y(2:end)).^2);
+        pos=find(distances>getfieldvalue(options,'Cutoff',1000));
+        x(pos,:)=NaN;
+        y(pos,:)=NaN;
+        for j=1:numcolors;
+                pos=find(colorind==j);
+                if(~isempty(pos) & pos(1)==1), pos(1)=[]; end
+                if ~isempty(pos),
+                        tempx = [x(pos-1) x(pos) NaN(size(pos))]';
+                        tempy = [y(pos-1) y(pos) NaN(size(pos))]';
+                        line(tempx(1:end-1),tempy(1:end-1),'color',palette(j,:),'linewidth',getfieldvalue(options,'LineWidth',2));
+                end
+        end
 end