This function implements a convolution for each pixel of datapack. 

  It convolves each pixel of signal with the function H of M elements, 
  so that if 
  H    = [h0 h1 h2 ... h(M-1)]
  then its Z transform is:
  M2=floor(M/2);
  $H[Z] = h0 Z^{+M2} + h1 Z^{+M2-1} + ...  h(M-1) Z^{+M2-(M-1)}$ 

  Thus H represents an impulse response with sample zero in h(M2).

  After starting the main routine just type the following command at the
  prompt:
  DATAOUT = datapack_conv(DATA,H)
  
  Input:
  DATA   is a speckle data pack. Where DATA is a 3D matrix created grouping NTIMES 
         intensity matrices with NLIN lines and NCOL columns. When N=size(DATA), then
         N(1,1) represents NLIN and
         N(1,2) represents NCOL and
         N(1,3) represents NTIMES.
  H      is a vector, where H represents a non-causal FIR filters centered in h(M2).

  
  Output:
  DATAOUT  is obtained from the output of H FIR filter, in non-causal form, 
           where its input is each pixel of datapack. 
  

  For help, bug reports and feature suggestions, please visit:
  http://nongnu.org/bsltl/

Package: bsltl