/trunk/libraries/libRaster/src/org/gvsig/raster/dataset/properties/DatasetStatistics.java - Annotate - | gvSIG desktop 1 - gvSIG

root / trunk / libraries / libRaster / src / org / gvsig / raster / dataset / properties / DatasetStatistics.java @ 11196

History | View | Annotate | Download (8.34 KB)

-            nacho
+/* gvSIG. Sistema de Informaci?n Geogr?fica de la Generalitat Valenciana
+ *
  * Copyright (C) 2006 IVER T.I. and Generalitat Valenciana.
+ *
  * This program is free software; you can redistribute it and/or
  * modify it under the terms of the GNU General Public License
  * as published by the Free Software Foundation; either version 2
  * of the License, or (at your option) any later version.
+ *
  * This program is distributed in the hope that it will be useful,
  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  * GNU General Public License for more details.
+ *
  * You should have received a copy of the GNU General Public License
  * along with this program; if not, write to the Free Software
  * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307,USA.
  */
 package org.gvsig.raster.dataset.properties;
 import java.io.File;
 import java.util.Date;
 import java.util.Hashtable;
 import org.gvsig.raster.dataset.FileNotOpenException;
 import org.gvsig.raster.dataset.IBuffer;
 import org.gvsig.raster.dataset.IStatistics;
 import org.gvsig.raster.dataset.InvalidSetViewException;
 import org.gvsig.raster.dataset.RasterDataset;
 import org.gvsig.raster.dataset.RasterDriverException;
 /**
  * Estadisticas asociadas a un fichero raster.
+ *
  * @author Nacho Brodin (nachobrodin@gmail.com)
  */
 public class DatasetStatistics implements IStatistics{
         /*
          * Esta a false si las estadisticas no son del fichero completo. Esto es posible porque podemos
          * tener unas estad?sticas calculadas a partir de una petici?n con subsampleo. Hay que tener en
          * cuenta que el raster puede ser muy grande y este calculo muy costoso.
          */
         protected boolean                                complete = false;
         protected double[]                         max = null;
         protected double[]                         min = null;
         protected double[]                         secondMax = null;
         protected double[]                         secondMin = null;
         protected double[]                         mean = null;
         protected double[]                         variance = null;
         protected String                                 fName = null;
         protected RasterDataset                grf = null;
         protected boolean                                calculated = false;
         protected Hashtable                        tailTrim = new Hashtable();
-            bsanchez
+        private boolean                                        canceled = false;
             nacho
         /**
          * Constructor. Asigna el fichero asociado.
          */
         public DatasetStatistics(RasterDataset grf){
                 this.grf = grf;
+        }
         /**
          * Asigna el valor m?ximo del grid
          * @return Valor m?ximo
          */
         public void setMax(double[] max) {
                 this.max = max;
+        }
         /**
          * Asigna el valor del segundo m?ximo
          * @return Valor del segundo m?ximo
          */
         public void setSecondMax(double[] smax) {
                 this.secondMax = smax;
+        }
         /**
          * Asigna el valor m?dio del grid
          * @return Valor medio
          */
         public void setMean(double[] mean) {
                 this.mean = mean;
+        }
         /**
          * Asigna el valor m?ximo del grid
          * @return Valor m?nimo
          */
         public void setMin(double[] min) {
                 this.min = min;
+        }
         /**
          * Asigna el valor del segundo m?nimo
          * @return Valor del segundo m?nimo
          */
         public void setSecondMin(double[] smin) {
                 this.secondMin = smin;
+        }
         /**
          * Asigna la varianza
          * @return Varianza
          */
         public void setVariance(double[] variance) {
                 this.variance = variance;
+        }
         /*
          *  (non-Javadoc)
          * @see org.gvsig.fmap.driver.IStatistics#getMax()
          */
         public double[] getMax() {
                 return max;
+        }
         /*
          *  (non-Javadoc)
          * @see org.gvsig.fmap.driver.IStatistics#getSecondMax()
          */
         public double[] getSecondMax() {
                 return secondMax;
+        }
         /*
          *  (non-Javadoc)
          * @see org.gvsig.fmap.driver.IStatistics#getSecondMin()
          */
         public double[] getSecondMin() {
                 return secondMin;
+        }
         /*
          *  (non-Javadoc)
          * @see org.gvsig.fmap.driver.IStatistics#getMaximun()
          */
         public double getMaximun(){
                 double m = Double.MIN_VALUE;
                 for(int i = 0; i < max.length; i++)
                         m = Math.max(m, max[i]);
                 return m;
+        }
         /*
          *  (non-Javadoc)
          * @see org.gvsig.fmap.driver.IStatistics#getMinimun()
          */
         public double getMinimun(){
                 double m = Double.MAX_VALUE;
                 for(int i = 0; i < min.length; i++)
                         m = Math.min(m, min[i]);
                 return m;
+        }
         /*
          *  (non-Javadoc)
          * @see org.gvsig.fmap.driver.IStatistics#getMean()
          */
         public double[] getMean() {
                 return mean;
+        }
         /*
          *  (non-Javadoc)
          * @see org.gvsig.fmap.driver.IStatistics#getMin()
          */
         public double[] getMin() {
                 return min;
+        }
         /*
          *  (non-Javadoc)
          * @see org.gvsig.fmap.driver.IStatistics#getVariance()
          */
         public double[] getVariance() {
                 return variance;
+        }
         /*
          *  (non-Javadoc)
          * @see org.gvsig.fmap.driver.IStatistics#getBandCount()
          */
         public int getBandCount(){
                 return grf.getBandCount();
+        }
         /*
          *  (non-Javadoc)
          * @see org.gvsig.fmap.driver.IStatistics#calcFullStatistics()
          */
         public void calcFullStatistics() throws FileNotOpenException, RasterDriverException{
                 long t2;
-            bsanchez
+                long t1 = new Date().getTime();
-            nacho
+                int type = grf.getDataType();
                 max = new double[grf.getBandCount()];
                 min = new double[grf.getBandCount()];
                 secondMax = new double[grf.getBandCount()];
                 secondMin = new double[grf.getBandCount()];
                 mean = new double[grf.getBandCount()];
                 variance = new double[grf.getBandCount()];
                 int iValues = 0;
                 byte[] b = null;short[] s = null;int[] i = null;float[] f = null;double[] d = null;
-            bsanchez
+                for (int iBand = 0; iBand < grf.getBandCount(); iBand ++){
-            nacho
+                        iValues = 0;
                         max[iBand] = Double.MIN_VALUE;
                         min[iBand] = Double.MAX_VALUE;
                         secondMax[iBand] = Double.MIN_VALUE;
                         secondMin[iBand] = Double.MAX_VALUE;
                         mean[iBand] = variance[iBand] = 0;
-            bsanchez
+                        for (int line = 0; line < grf.getHeight(); line ++) {
-            nacho
+                                Object buf = null;
                                 try {
                                         buf = grf.readCompleteLine(line, iBand);
                                 } catch (InvalidSetViewException e) {
                                         //La vista se asigna autom?ticamente
+                                }
                                 switch(type){
                                 case IBuffer.TYPE_BYTE:                b = (byte[])buf;break;
                                 case IBuffer.TYPE_SHORT:         s = (short[])buf;break;
                                 case IBuffer.TYPE_INT:                 i = (int[])buf;break;
                                 case IBuffer.TYPE_FLOAT:         f = (float[])buf;break;
                                 case IBuffer.TYPE_DOUBLE:         d = (double[])buf;break;
+                                }
-            bsanchez
+                                for (int col = 0; col < grf.getWidth(); col ++){
-            nacho
+                                        double z = (b != null) ? (b[col] & 0xff) : ((s != null) ? (s[col] & 0xffff) : ((i != null) ? (i[col] & 0xffffff)  : (f != null) ? f[col] : d[col]));
-            bsanchez
+                                        if (line == 0 && col == 0){
-            nacho
+                                                min[iBand] = z;
                                                 secondMin[iBand] = z;
                                                 max[iBand] = z;
                                                 secondMax[iBand] = z;
                                                 continue;
+                                        }
-            bsanchez
+                                        if ( min[iBand] > z ){
-            nacho
+                                                secondMin[iBand] = min[iBand];
                                                 min[iBand] = z;
+                                        }
-            bsanchez
+                                        if ( secondMin[iBand] == min[iBand] && min[iBand] < z)
-            nacho
+                                                secondMin[iBand] = z;
-            bsanchez
+                                        if ( max[iBand] < z ){
-            nacho
+                                                secondMax[iBand] = max[iBand];
                                                 max[iBand] = z;
+                                        }
-            bsanchez
+                                        if ( secondMax[iBand] == max[iBand] && max[iBand] > z)
-            nacho
+                                                secondMax[iBand] = z;
                                         mean[iBand] += z;
                                         variance[iBand] += z * z;
                                         iValues++;
             bsanchez
                                 if (isCanceled())
                                         return;
             nacho
                         if( iValues > 0 ){
                                 mean[iBand] /= (double) iValues;
                                 variance[iBand] = variance[iBand] / (double) iValues - mean[iBand] * mean[iBand];
+                        }
+                }
                 calculated = true;
                 t2 = new Date().getTime();
             System.out.println("Estadisticas " + grf.getFName() + ": " + ((t2 - t1) / 1000D) + ", secs.");
+        }
         /**
          * Carga estadisticas desde el fichero Rmf si las hay.
          * @param fName Nombre del fichero
          */
         public void loadFromRmf(){
                 if(grf != null){
                         File f = new File(grf.getFName());
                         if(!f.exists())
                                 return;
                         //calculated = true;
                         //TODO: FUNCIONALIDAD: Carga de estad?sticas rmf
+                }
+        }
         /**
          * Salva estad?sticas a fichero rmf.
          * @param fName
          */
         public void saveToRmf(){
                 //TODO: FUNCIONALIDAD: Salvado de estad?sticas a rmf
+        }
         /*
          *  (non-Javadoc)
          * @see org.gvsig.fmap.driver.IStatistics#isCalculated()
          */
         public boolean isCalculated() {
                 return calculated;
+        }
         /*
          *  (non-Javadoc)
          * @see org.gvsig.fmap.driver.IStatistics#setTailTrimValue(double, java.lang.Object)
          */
         public void setTailTrimValue(double percent, Object valueByBand){
                 String s = new Double(percent).toString();
                 tailTrim.put(s, valueByBand);
+        }
         /*
          *  (non-Javadoc)
          * @see org.gvsig.fmap.driver.IStatistics#getTailTrimValue(double)
          */
         public Object getTailTrimValue(double percent){
                 String s = new Double(percent).toString();
                 return tailTrim.get(s);
+        }
             bsanchez
         /*
          * (non-Javadoc)
          * @see org.gvsig.raster.util.ICancellable#isCanceled()
          */
         public boolean isCanceled() {
                 return canceled;
+        }
         /*
          * (non-Javadoc)
          * @see org.gvsig.raster.util.ICancellable#setCanceled(boolean)
          */
         public void setCanceled(boolean value) {
                 canceled = value;
+        }
+}