Changeset 581 for trunk/src/plot_field.m

Timestamp:

Mar 12, 2013, 12:52:13 PM (11 years ago)

Author:

sommeria

Message:

clean the transform field functions

File:

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

trunk/src/plot_field.m

@@ -249 +249 @@
 function errormsg=plot_text(FieldData,CellInfo,htext)
 %-------------------------------------------------------------------
-errormsg=[];
+errormsg='';
 txt_cell={};
 Data={};
@@ -1045 +1045 @@
     
     %decimate by a factor 2 in vector mesh(4 in nbre of vectors)
+    check_decimate=0;
     if isfield(PlotParam.Vectors,'CheckDecimate4') && PlotParam.Vectors.CheckDecimate4
+        check_decimate=1;
         diffy=diff(vec_Y); %difference dy=vec_Y(i+1)-vec_Y(i)
         dy_thresh=max(abs(diffy))/2; 
@@ -1055 +1057 @@
         nb_sel=length(ind_sel);
         ind_sel=ind_sel(1:2:nb_sel);% take half the points on a line
+    elseif isfield(PlotParam.Vectors,'CheckDecimate16') && PlotParam.Vectors.CheckDecimate16
+        check_decimate=1;
+        diffy=diff(vec_Y); %difference dy=vec_Y(i+1)-vec_Y(i)
+        dy_thresh=max(abs(diffy))/2; 
+        ind_jump=find(abs(diffy) > dy_thresh); %indices with diff(vec_Y)> max/2, detect change of line
+        ind_sel=1:ind_jump(1);%select the first line
+        for i=2:4:length(ind_jump)-1
+            ind_sel=[ind_sel (ind_jump(i)+1:ind_jump(i+1))];% select the odd lines
+        end
+        nb_sel=length(ind_sel);
+        ind_sel=ind_sel(1:4:nb_sel);% take half the points on a line
+    end
+    if check_decimate
         vec_X=vec_X(ind_sel);
         vec_Y=vec_Y(ind_sel);