Changeset 911 for trunk/src/series
- Timestamp:
- Jun 15, 2015, 10:45:43 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/series/civ_series.m
r910 r911 1093 1093 xtable(ivec)=iref+utable(ivec)/2-0.5;% convec flow (velocity taken at the point middle from imgae 1 and 2) 1094 1094 ytable(ivec)=jref+vtable(ivec)/2-0.5;% and position of pixel 1=0.5 (convention for image coordinates=0 at the edge) 1095 iref=round(xtable(ivec) );% nearest image index for the middle of the vector1096 jref=round(ytable(ivec) );1095 iref=round(xtable(ivec)+0.5);% nearest image index for the middle of the vector 1096 jref=round(ytable(ivec)+0.5); 1097 1097 % eliminate vectors located in the mask 1098 if checkmask && par_civ.Mask(jref,iref)<200 && par_civ.Mask(jref,iref)>=1001098 if checkmask && (iref<1 || jref<1 ||iref>npx_ima || jref>npy_ima ||( par_civ.Mask(jref,iref)<200 && par_civ.Mask(jref,iref)>=100)) 1099 1099 utable(ivec)=0; 1100 1100 vtable(ivec)=0; … … 1129 1129 F=0; 1130 1130 [npy,npx]=size(result_conv); 1131 result_conv(result_conv<1)=1; %set to 1 correlation values smaller than 1 (to avoid divergence in the log)1131 result_conv(result_conv<1)=1; %set to 1 correlation values smaller than 1 (=0 by discretisation, to avoid divergence in the log) 1132 1132 %the following 8 lines are copyright (c) 1998, Uri Shavit, Roi Gurka, Alex Liberzon, Technion ï¿œ Israel Institute of Technology 1133 1133 %http://urapiv.wordpress.com 1134 1134 peaky = y; 1135 if y < = npy-1&& y > 11135 if y < npy && y > 1 1136 1136 f0 = log(result_conv(y,x)); 1137 1137 f1 = log(result_conv(y-1,x)); … … 1142 1142 end 1143 1143 peakx=x; 1144 if x < =npx-1 && x > 11144 if x < npx-1 && x > 1 1145 1145 f0 = log(result_conv(y,x)); 1146 1146 f1 = log(result_conv(y,x-1)); … … 1162 1162 result_conv(result_conv<1)=1; %set to 1 correlation values smaller than 1 (to avoid divergence in the log) 1163 1163 [npy,npx]=size(result_conv); 1164 if (x < = npx-1) && (y <= npy-1) && (x >= 1) && (y >=1)1164 if (x < npx) && (y < npy) && (x > 1) && (y > 1) 1165 1165 F=0; 1166 1166 for i=-1:1
Note: See TracChangeset
for help on using the changeset viewer.