Application: gvSIG desktop » gvSIG tools: ToolBox

package es.unex.sextante.gridStatistics.neighborhoodSkewness;

import es.unex.sextante.core.Sextante;

import es.unex.sextante.gridStatistics.base.NeighborhoodStatsBaseAlgorithm;

public class NeighborhoodSkewnessAlgorithm

         extends

            NeighborhoodStatsBaseAlgorithm {

   @Override

   public void defineCharacteristics() {

      setName(Sextante.getText("Skeweness__neighbourhood"));

      setGroup(Sextante.getText("Focal_statistics"));

      super.defineCharacteristics();

   }

   @Override

   protected double processValues() {

      int i;

      int iValidCells = 0;

      double dValue;

      double dMean = 0;

      double dDif;

      double dDif2 = 0;

      double dDif3 = 0;

      for (i = 0; i < m_dValues.length; i++) {

         dValue = m_dValues[i];

         if (dValue != NO_DATA) {

            dMean += dValue;

            iValidCells++;

         }

      }

      dMean /= iValidCells;

      if (iValidCells > 2) {

         for (i = 0; i < m_dValues.length; i++) {

            dDif = m_dValues[i] - dMean;

            dDif2 += (dDif * dDif);

            dDif3 += (dDif * dDif * dDif);

         }

         dDif2 /= iValidCells;

         dDif3 /= iValidCells;

         return (dDif3 / Math.pow(dDif2, 1.5));

      }

      else {

         return NO_DATA;

      }

   }

}