Changeset 883 for trunk/src/series

Timestamp:

Mar 13, 2015, 8:19:39 PM (10 years ago)

Author:

sommeria

Message:

bugs corrected in series, transform functions with input parameters introduced, civ/deformation corrected (sign)

Location:

trunk/src/series

Files:

• aver_stat.m (modified) (1 diff)
• civ_input.m (modified) (10 diffs)
• civ_series.m (modified) (2 diffs)
• merge_proj.m (modified) (1 diff)
• time_series.m (modified) (1 diff)

trunk/src/series/aver_stat.m

@@ -185 +185 @@
     transform_fct=str2func(Param.FieldTransform.TransformName);
     rmpath(Param.FieldTransform.TransformPath)
+    if isfield(Param,'TransformInput')
+            XmlData{1}.TransformInput=Param.TransformInput;
+    end     
 end
 

trunk/src/series/civ_input.m

@@ -688 +688 @@
 nbfield=siztime(1)-1;
 nbfield2=siztime(2)-1;
-indchosen=1;  %%first pair selected by default
+%indchosen=1;  %%first pair selected by default
 %displ_num used to define the indices of the civ_input pairs
 % in mode 'pair j1-j2', j1 and j2 are the file indices, else the indices
@@ -907 +907 @@
 nbpair=200;%default
 select=ones(size(1:nbpair));%flag for displayed pairs =1 for display
-nbpair=200; %default
+%nbpair=200; %default
 
 %% determine the menu display in .ListPairCiv1
@@ -915 +915 @@
             if select(ipair)
                 displ_pair{ipair}=['Di= ' num2str(-floor(ipair/2)) '|' num2str(ceil(ipair/2))];
+                displ_pair_dt{ipair}=displ_pair{ipair};
                 if ~checkframe
                     if size(Time,1)>=ref_i+1+ceil(ipair/2) && size(Time,2)>=ref_j+1&& ref_i-floor(ipair/2)>=0 && ref_j>=0
                         dt=Time(ref_i+1+ceil(ipair/2),ref_j+1)-Time(ref_i+1-floor(ipair/2),ref_j+1);%Time interval dtref_j+1
-                        displ_pair{ipair}=[displ_pair{ipair} ' :dt= ' num2str(dt*1000)];
+                        displ_pair_dt{ipair}=[displ_pair_dt{ipair} ' :dt= ' num2str(dt*1000)];
                     end
                 else
                     dt=ipair/1000;
-                    displ_pair{ipair}=[displ_pair{ipair} ' :dt= ' num2str(ipair)];
+                    displ_pair_dt{ipair}=[displ_pair_dt{ipair} ' :dt= ' num2str(ipair)];
                 end
             else
-                displ_pair{ipair}='...'; %pair not displayed in the menu
+                displ_pair{ipair}='...';
+                displ_pair_dt{ipair}='...'; %pair not displayed in the menu
             end
         end
@@ -932 +934 @@
             if select(ipair)
                 displ_pair{ipair}=['Dj= ' num2str(-floor(ipair/2)) '|' num2str(ceil(ipair/2))];
+                displ_pair_dt{ipair}=displ_pair{ipair};
                 if ~checkframe
                     if size(Time,2)>=ref_j+1+ceil(ipair/2) && size(Time,1)>=ref_i+1 && ref_j-floor(ipair/2)>=0 && ref_i>=0
                         dt=Time(ref_i+1,ref_j+1+ceil(ipair/2))-Time(ref_i+1,ref_j+1-floor(ipair/2));%Time interval dtref_j+1
-                        displ_pair{ipair}=[displ_pair{ipair} ' :dt= ' num2str(dt*1000)];
+                        displ_pair_dt{ipair}=[displ_pair_dt{ipair} ' :dt= ' num2str(dt*1000)];
                     end
                 else
                     dt=ipair/1000;
-                    displ_pair{ipair}=[displ_pair{ipair} ' :dt= ' num2str(dt*1000)];
+                    displ_pair_dt{ipair}=[displ_pair_dt{ipair} ' :dt= ' num2str(dt*1000)];
                 end
             else
                 displ_pair{ipair}='...'; %pair not displayed in the menu
+                displ_pair_dt{ipair}='...'; %pair not displayed in the menu
             end
         end
@@ -956 +960 @@
                 index_pair=index_pair+1;
                 displ_pair{index_pair}=['j= ' num2stra(numod_a,nom_type_ima) '-' num2stra(numod_b,nom_type_ima)];
+                displ_pair_dt{ipair}=displ_pair{ipair};
                 dt(index_pair)=numod_b-numod_a;%default dt
                 if size(Time,1)>ref_i && size(Time,2)>numod_b  % && ~checkframe
                     dt(index_pair)=Time(ref_i+1,numod_b+1)-Time(ref_i+1,numod_a+1);% Time interval dt
-                    displ_pair{index_pair}=[displ_pair{index_pair} ' :dt= ' num2str(dt(index_pair)*1000)];
+                    displ_pair_dt{index_pair}=[displ_pair_dt{index_pair} ' :dt= ' num2str(dt(index_pair)*1000)];
                 end
             end
@@ -966 +971 @@
         [tild,indsort]=sort(dt);
         displ_pair=displ_pair(indsort);
+        displ_pair_dt=displ_pair_dt(indsort);
     case 'displacement'
         displ_pair={'Di=Dj=0'};
+        displ_pair_dt={'Di=Dj=0'};
 end
 if index==1
-    set(handles.ListPairCiv1,'String',displ_pair');
+    set(handles.ListPairCiv1,'String',displ_pair_dt');
 end
 
@@ -976 +983 @@
 %ichoice=find(select,1);% index of first selected pair
 %if (isempty(ichoice) || ichoice < 1); ichoice=1; end;
+end_pair=regexp(PairCiv1Init,' :dt=');
+if ~isempty(end_pair)
+    PairCiv1Init=PairCiv1Init(1:end_pair-1);
+end
 ichoice=find(strcmp(PairCiv1Init,displ_pair'),1);
 if ~isempty(ichoice)
@@ -988 +999 @@
 %% determine the default selection in the pair menu for Civ2
 if strcmp(get(handles.ListPairCiv2,'Visible'),'on')
+    end_pair=regexp(PairCiv2Init,' :dt=');
+    if ~isempty(end_pair)
+        PairCiv2Init=PairCiv2Init(1:end_pair-1);
+    end
     ichoice=find(strcmp(PairCiv2Init,displ_pair'),1);
     if ~isempty(ichoice)
@@ -994 +1009 @@
         set(handles.ListPairCiv2,'Value',1)
     end
-    
-    %
-    %     initial=get(handles.ListPairCiv2,'Value');
-    %
-    %
-    %     if initial>length(displ_pair')%|~isequal(select(initial),1)
-    %         if ichoice <= length(displ_pair')
-    %             set(handles.ListPairCiv2,'Value',ichoice);% same pair proposed by default for civ2
-    %         else
-    %             set(handles.ListPairCiv2,'Value',1);% same pair proposed by default for civ2
-    %         end
-    %     end
 else
     set(handles.ListPairCiv2,'Value',get(handles.ListPairCiv1,'Value'))% initiate the choice of Civ2 as a reproduction of if civ1
 end
-set(handles.ListPairCiv2,'String',displ_pair');
+set(handles.ListPairCiv2,'String',displ_pair_dt');
 set(gcf,'Pointer','arrow')% Indicate that the process is finished
-
-
-    
-% %------------------------------------------------------------------------   
-% % call 'view_field.fig' to display the  field selected in the list of 'status'
-% function open_view_field(hObject, eventdata)
-% %------------------------------------------------------------------------
-% list=get(hObject,'String');
-% index=get(hObject,'Value');
-% rootroot=get(hObject,'UserData');
-% filename=list{index};
-% ind_dot=strfind(filename,'...');
-% filename=filename(1:ind_dot-1);
-% filename=fullfile(rootroot,filename);
-% delete(get(hObject,'parent'))%delete the display figure to stop the check process
-% if exist(filename,'file')%visualise the vel field if it exists
-%     uvmat(filename)
-%     set(gcbo,'Value',1)
-% end
 
 
@@ -1160 +1144 @@
 set(handles.ConfigSource,'String','NEW')
 set(handles.OK,'BackgroundColor',[1 0 1])
+if get(handles.CheckDeformation,'Value')
+    set(handles.num_CorrSmooth,'Visible','off')
+else
+    set(handles.num_CorrSmooth,'Visible','on')
+end
 
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

trunk/src/series/civ_series.m

@@ -1036 +1036 @@
                     yi=(1:mesh:size(image1_crop,1))';
                     [XI,YI]=meshgrid(xi-ceil(size(image1_crop,2)/2),yi-ceil(size(image1_crop,1)/2));
-                    XIant=XI-par_civ.DUDX(ivec)*XI-par_civ.DUDY(ivec)*YI+ceil(size(image1_crop,2)/2);
-                    YIant=YI-par_civ.DVDX(ivec)*XI-par_civ.DVDY(ivec)*YI+ceil(size(image1_crop,1)/2);
+                    XIant=XI-par_civ.DUDX(ivec)*XI+par_civ.DUDY(ivec)*YI+ceil(size(image1_crop,2)/2);
+                    YIant=YI+par_civ.DVDX(ivec)*XI-par_civ.DVDY(ivec)*YI+ceil(size(image1_crop,1)/2);
                     image1_crop=interp2(image1_crop,XIant,YIant);
                     image1_crop(isnan(image1_crop))=0;
@@ -1044 +1044 @@
                     image2_crop=interp2(image2_crop,xi,yi,'*spline');
                     image2_crop(isnan(image2_crop))=0;
-                    %par_civ.CorrSmooth=3;%%%%%%%%%%%%%%%%%%%
-                    %%
+                    par_civ.CorrSmooth=2;% use SUBPIX2DGAUSS (take into account more points near the max)
                 end
                 sum_square=sum(sum(image1_crop.*image1_crop));

trunk/src/series/merge_proj.m

@@ -179 +179 @@
         transform_fct=str2func(Param.FieldTransform.TransformName);
         cd (currentdir)
+        if isfield(Param,'TransformInput')
+            for iview=1:NbView
+            XmlData{iview}.TransformInput=Param.TransformInput;
+            end
+        end       
 end
 %%%%%%%%%%%% END STANDARD PART  %%%%%%%%%%%%

trunk/src/series/time_series.m

@@ -220 +220 @@
     transform_fct=str2func(Param.FieldTransform.TransformName);
     rmpath(Param.FieldTransform.TransformPath)
+    if isfield(Param,'TransformInput')
+        XmlData{1}.TransformInput=Param.TransformInput;
+    end
 end