The probability mass function of absolute difference (PMFAD) represents the 
  probabilities of a random variable Z, being $Z=|i-j|$.
 
  $Pr(z+1)= P(|i-j|=z)$, $0 \leq z \leq 255$

  Where i is one intensity id line in COM matrix and j is one intensity id column 
  in COM. Thus $Pr(z+1)= P(|i-j|=z)$ is the probability of happen an absolute
  intensity jump of value z. This function calculates the difference probability [1].
  The COM matrix needs a matrix of 256x256 elements.

  References:
  [1]  R.M. Haralick, K. Shanmugam, and Its' Hak Dinstein. 'Textural Features 
       for Image Classification'. In: Systems, Man and Cybernetics, IEEE 
       Transactions on SMC-3.6 (Nov. 1973), pages 610-621. ISSN: 0018-9472. 

  After starting the main routine just type the following command at the
  prompt:
  Pr = pmfad(COM);
  
  Input:
  COM  is the Co-Occurrence matrix. A 2D matrix with 256 lines and 256 columns.

  Output:
  Pr   is the probability mass function, where  $Pr(z+1)= P(|i-j|=z)$ is the 
       probability of happen an absolute intensity jump of value z. 0<=z<=255
       Pr is a vector with 256 elements.
  Z    [OPTIONAL] is a vector with the absolute intensities jumps $|i-j|=z$.
       Z=[0:255];


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

Package: bsltl