Changeset 153 for trunk/src/civ.m
- Timestamp:
- Dec 9, 2010, 11:31:04 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/civ.m
r150 r153 68 68 testall=0; 69 69 %default input parameters: 70 num 1=1; % set of field i numbers71 num 2=1; % set of field i numbers72 num_ a=1; % set of field j numbers (fields a)73 num_ b=1; % second set of field j numbers (fields b)70 num_i1=1; % set of field i numbers 71 num_i2=1; % set of field i numbers 72 num_j1=1; % set of field j numbers (fields a) 73 num_j2=1; % second set of field j numbers (fields b) 74 74 subdir='A'; % subdir for the netcdf result files 75 75 ind_opening=1; % proposed operation number (1=civ1,2=fix1,3=patch1,4=civ2,5=fix2,6=patch2) … … 78 78 filebase=param.RootName; 79 79 nom_type_read=param.NomType; 80 num1=param.num1; 81 num2=param.num2; 82 num_a=param.num_a; 83 num_b=param.num_b; 80 num_i1=param.num1; 81 if isnan(num_i1),num_i1=1;end 82 num_i2=param.num2; 83 if isnan(num_i2),num_i2=num_i1;end 84 num_j1=param.num_a; 85 if isnan(num_j1),num_j1=1;end 86 num_j2=param.num_b; 87 if isnan(num_j2),num_j2=num_j1;end 84 88 subdir=param.SubDir; 85 89 ind_opening=param.IndOpening; 86 90 ext=param.ImaExt; 87 91 end 88 browse.num_i1=num 1;89 browse.num_i2=num 2;90 browse.num_j1=num_ a;91 browse.num_j2=num_ b;92 browse.num_i1=num_i1; 93 browse.num_i2=num_i2; 94 browse.num_j1=num_j1; 95 browse.num_j2=num_j2; 92 96 if ~isempty(ext) && (~isempty(imformats(ext(2:end)))||strcmpi(ext,'.avi'));%if an image file has been opened by uvmat 93 97 set(handles.ImaExt,'String',ext) … … 195 199 196 200 % set the range of fields (1:1 by default) and selected pair 197 if is nan(num2)||isequal(num2,num1)198 num_ref_i=num 1;199 else 200 num_ref_i=floor((num 1+num2)/2);201 browse.incr_pair(1)=num 2-num1;201 if isequal(num_i2,num_i1) 202 num_ref_i=num_i1; 203 else 204 num_ref_i=floor((num_i1+num_i2)/2); 205 browse.incr_pair(1)=num_i2-num_i1; 202 206 browse.incr_pair(2)=0; 203 207 end 204 if is nan(num_b)||isequal(num_a,num_b)205 if isnan(num_ a)208 if isequal(num_j1,num_j2) 209 if isnan(num_j1) 206 210 num_ref_j=1; 207 211 else 208 num_ref_j=num_ a;209 end 210 else 211 num_ref_j=floor((num_ a+num_b)/2);212 browse.incr_pair(2)=num_ b-num_a;212 num_ref_j=num_j1; 213 end 214 else 215 num_ref_j=floor((num_j1+num_j2)/2); 216 browse.incr_pair(2)=num_j2-num_j1; 213 217 end 214 218 set(handles.first_i,'String',num2str(num_ref_i)); … … 658 662 % nbdetect=0;%test of detected images 659 663 field_i=browse.num_i2; 664 imagename=name_generator(filebase,field_i,1,ext_search,nom_type_search); 665 imagename_plus=''; 660 666 idetect=1; 661 while idetect ==1%look for the maximum file number in the series662 imagename =name_generator(filebase,field_i+1,1,ext_search,nom_type_search);663 idetect=(exist(imagename ,'file')==2);667 while idetect %look for the maximum file number in the series 668 imagename_plus=name_generator(filebase,field_i+1,1,ext_search,nom_type_search); 669 idetect=(exist(imagename_plus,'file')==2)&& ~strcmp(imagename,imagename_plus); 664 670 if idetect 665 671 field_i=field_i+1; … … 670 676 nbfield=field_i;% last detected field number 671 677 field_i=browse.num_i1;%look for the minimum file number in the series 678 imagename_min=''; 672 679 idetect=1; 673 680 while idetect==1 674 imagename =name_generator(filebase,field_i-1,1,ext_search,nom_type_search);675 idetect=(exist(imagename ,'file')==2);681 imagename_min=name_generator(filebase,field_i-1,1,ext_search,nom_type_search); 682 idetect=(exist(imagename_min,'file')==2)&& ~strcmp(imagename,imagename_min); 676 683 if idetect 677 684 field_i=field_i-1; … … 682 689 field_i=browse.num_i1; 683 690 field_j=browse.num_j2; 691 imagename_plus=''; 684 692 jdetect=1; 685 693 while jdetect==1 %look for the maximum file number in the series 686 imagename =name_generator(filebase,field_i,field_j,ext_search,nom_type_search);687 jdetect=(exist(imagename ,'file')==2);694 imagename_plus=name_generator(filebase,field_i,field_j+1,ext_search,nom_type_search); 695 jdetect=(exist(imagename_plus,'file')==2)&& ~strcmp(imagename,imagename_plus); 688 696 if jdetect 689 697 field_j=field_j+1; … … 2479 2487 vers=0; 2480 2488 subdir_civ1_new=subdir_civ1; 2481 ind_test=0; 2482 while detect==1 && ind_test<10%create a new subdir if the netcdf files already exist 2489 while detect==1 %create a new subdir if the netcdf files already exist 2483 2490 for ifile=1:nbfield 2484 2491 for j=1:nbslice
Note: See TracChangeset
for help on using the changeset viewer.