Not implemented in current version This routine filters the backscatter data. (adapted from code by J. Czuba) Version for Ricardo for use with intensity data rather than backscatter P.R. Jackson, USGS, 2-15-11
0001 function A = VMT_FilterBS_IntensityRS(z,A) 0002 % Not implemented in current version 0003 % This routine filters the backscatter data. 0004 % 0005 % (adapted from code by J. Czuba) 0006 % 0007 % Version for Ricardo for use with intensity data rather than backscatter 0008 % 0009 % P.R. Jackson, USGS, 2-15-11 0010 0011 0012 %% Filter 0013 % filter backscatter when climbing walls 0014 0015 for zi = 1 : z 0016 0017 A(zi).Clean.backstandf=nan(double(A(zi).Sup.nBins),A(zi).Sup.noe); 0018 A(zi).Clean.bsf=nan(double(A(zi).Sup.nBins),A(zi).Sup.noe); 0019 0020 % Determine the standard deviation of the backscatter for the four beams in 0021 % one bin in one ensemble 0022 back=A(zi).Wat.backscatter; 0023 back(back>=255) = NaN; 0024 A(zi).Clean.backstandf = std(back,0,3); 0025 0026 % Remove backscatter intensities >= 255 (bad value) and backscatter 0027 % intensities with a standard deviation greater than 10 (mainly climbing 0028 % walls) 0029 A(zi).Clean.bsf = nanmean(A(zi).Wat.backscatter,3); 0030 %A(zi).Clean.bsf(A(zi).Clean.bsf>=255) = NaN; 0031 %A(zi).Clean.bsf(A(zi).Clean.backstandf>10)=NaN; 0032 0033 A(zi).Nav.depth(A(zi).Nav.depth==-32768)=NaN; 0034 0035 end