# VMT_FilterBS_IntensityRS

## PURPOSE

Not implemented in current version

## SYNOPSIS

function A = VMT_FilterBS_IntensityRS(z,A)

## DESCRIPTION

     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

## CROSS-REFERENCE INFORMATION

This function calls:

This function is called by:

## SOURCE CODE

    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