Application: gvSIG desktop » gvSIG raster

/* gvSIG. Geographic Information System of the Valencian Government

 *

 * Copyright (C) 2007-2008 Infrastructures and Transports Department

 * of the Valencian Government (CIT)

 *

 * 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., 51 Franklin Street, Fifth Floor, Boston,

 * MA  02110-1301, USA.

 *

 */

package org.gvsig.raster.impl.store.properties;

import java.util.ArrayList;

import org.gvsig.fmap.dal.coverage.exception.FileNotOpenException;

import org.gvsig.fmap.dal.coverage.exception.ProcessInterruptedException;

import org.gvsig.fmap.dal.coverage.exception.RasterDriverException;

import org.gvsig.fmap.dal.coverage.exception.RmfSerializerException;

import org.gvsig.fmap.dal.coverage.store.MultiRasterStore;

import org.gvsig.fmap.dal.coverage.store.props.Statistics;

import org.gvsig.raster.impl.provider.RasterProvider;

import org.gvsig.raster.impl.store.DefaultMultiRasterStore;

/**

 * Estadisticas asociadas a un fichero raster formado por multiples ficheros.

 *

 * @author Nacho Brodin (nachobrodin@gmail.com)

 */

public class MultiDatasetStatistics extends AbstractStatistics {

        protected DatasetStatistics[]     statList      = null;

        protected RasterProvider[]        providerList  = null;

        protected int                     nDataset      = 0;

        /**

         * Carga el objeto DatasetListStatistics desde los datos del RMF.

         * @param datasource

         * @return

         * @throws RmfSerializerException

         */

        public static MultiDatasetStatistics loadDatasetListStatistics(MultiRasterStore store) throws RmfSerializerException {

                DatasetStatistics[] stats = new DatasetStatistics[store.getDataStoreCount()];

                DefaultMultiRasterStore storeImpl = (DefaultMultiRasterStore)store;

                for (int i = 0; i < storeImpl.getProviders().size(); i++) {

                        Statistics statFile = (storeImpl.getProvider(i)).getStatistics();

                        stats[i] = (DatasetStatistics)(storeImpl.getProvider(i)).loadObjectFromRmf(DatasetStatistics.class, statFile);

                        //Con que un dataset no tenga la estadistica calculada se pone a no calculado el DatasetListStatistics

                        if(!stats[i].isCalculated())

                                return null;

                }

                MultiDatasetStatistics result = new MultiDatasetStatistics(stats);

                result.setCalculated(true);

                return result;

        }

        /**

         * Constructor

         */

        public MultiDatasetStatistics() {

        }

        /**

         * Constructor

         */

        public MultiDatasetStatistics(RasterProvider provider) {

                this.providerList = new RasterProvider[]{provider};

                statList = new DatasetStatistics[]{provider.getStatistics()};

        }

        /**

         * Constructor

         */

        public MultiDatasetStatistics(RasterProvider[] datasetList) {

                this.providerList = datasetList;

                statList = new DatasetStatistics[datasetList.length];

                for(int i = 0; i < datasetList.length; i ++)

                        statList[i] = this.providerList[i].getStatistics();

        }

        /**

         * Constructor

         */

        public MultiDatasetStatistics(ArrayList<RasterProvider> datasetList) {

                statList = new DatasetStatistics[datasetList.size()];

                this.providerList = new RasterProvider[datasetList.size()];

                for(int i = 0; i < datasetList.size(); i ++){

                        this.providerList[i] = datasetList.get(i);

                        if(this.providerList[i] != null)

                                statList[i] = this.providerList[i].getStatistics();

                }

        }

        /**

         * Constructor

         */

        public MultiDatasetStatistics(DatasetStatistics[] statList) {

                this.statList = statList;

                providerList = new RasterProvider[statList.length];

                for(int i = 0; i < statList.length; i ++)

                        providerList[i] = this.statList[i].getDataProvider();

                int len = 0;

                for(nDataset = 0; nDataset < statList.length; nDataset ++) {

                        try {

                                statList[nDataset].calcFullStatistics();

                        } catch (FileNotOpenException e) {

                        } catch (RasterDriverException e) {

                        } catch (ProcessInterruptedException e) {

                        }

                        len += providerList[nDataset].getBandCount();

                }

                constructStats(len);

        }

        /**

         * Adds a new data store to the statistics list

         * @param provider

         */

        public void addProvider(RasterProvider provider) {

                DatasetStatistics[] statListAux = new DatasetStatistics[providerList.length + 1];

                RasterProvider[] datasetListAux = new RasterProvider[providerList.length + 1];

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

                        datasetListAux[i] = providerList[i];

                        statListAux[i] = statList[i];

                }

                datasetListAux[providerList.length] = provider;

                statListAux[providerList.length] = provider.getStatistics();

                providerList = datasetListAux;

                statList = statListAux;

        }

        /**

         * Calcula las estadisticas recorriendo todo el fichero.

         */

        public void calcFullStatistics()throws FileNotOpenException, RasterDriverException, ProcessInterruptedException {

//                long t2;

//                long t1 = new Date().getTime();

                int len = 0;

                for(nDataset = 0; nDataset < statList.length; nDataset ++) {

                        statList[nDataset].calcFullStatistics();

                        len += providerList[nDataset].getBandCount();

                }

                try {

                        constructStats(len);

                } catch (ArrayIndexOutOfBoundsException e) {

                        throw new RasterDriverException("Error en el acceso al array de m?ximos y m?nimos");

                }

                calculated = true;

//                t2 = new Date().getTime();

//                System.out.println("Estadisticas MultiFile: " + ((t2 - t1) / 1000D) + ", secs.");

        }

        /*

         * (non-Javadoc)

         * @see org.gvsig.raster.impl.dataset.properties.AbstractStatistics#resetPercent()

         */

        public void resetPercent() {

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

                        statList[i].resetPercent();

                }

        }

        /*

         * (non-Javadoc)

         * @see org.gvsig.raster.impl.dataset.properties.AbstractStatistics#getPercent()

         */

        public int getPercent() {

                try {

                        return (statList[nDataset].getPercent() / statList.length);

                } catch (ArrayIndexOutOfBoundsException e) {

                        return 0;

                }

        }

        /**

         * Cuando se llama a este m?todo fuerza que la siguiente petici?n de estad?sticas

         * no sea le?da de RMF y sean recalculadas por completo.

         * @param forceToRecalc

         */

        public void forceToRecalc() {

                for(int i = 0; i < statList.length; i ++)

                        statList[i].forceToRecalc();

        }

        /**

         * Construye todos los datos del DatasetStadistics actual a partir de toda la

         * lista de DatasetStatistics que tiene este objeto.

         *

         * @param len

         * @throws ArrayIndexOutOfBoundsException

         */

        private void constructStats(int len) throws ArrayIndexOutOfBoundsException {

                max = new double[len];

                min = new double[len];

                secondMax = new double[len];

                secondMin = new double[len];

                maxByteUnsigned = new double[len];

                minByteUnsigned = new double[len];

                secondMaxByteUnsigned = new double[len];

                secondMinByteUnsigned = new double[len];

                mean = new double[len];

                variance = new double[len];

                int count = 0;

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

                        if(statList[i].getMax() == null)

                                return;

                        for(int j = 0; j < statList[i].getMax().length; j ++) {

                                max[count] = statList[i].getMax()[j];

                                min[count] = statList[i].getMin()[j];

                                secondMax[count] = statList[i].getSecondMax()[j];

                                secondMin[count] = statList[i].getSecondMin()[j];

                                maxByteUnsigned[count] = statList[i].getMaxByteUnsigned()[j];

                                minByteUnsigned[count] = statList[i].getMinByteUnsigned()[j];

                                secondMaxByteUnsigned[count] = statList[i].getSecondMaxByteUnsigned()[j];

                                secondMinByteUnsigned[count] = statList[i].getSecondMinByteUnsigned()[j];

                                mean[count] = statList[i].getMean()[j];

                                variance[count] = statList[i].getVariance()[j];

                                count ++;

                        }

                }

        }

        /**

         * Obtiene el flag que informa de si las estad?sticas est?n calculadas o no.

         * @return true indica que est?n calculadas y false que no lo est?n

         */

        public boolean isCalculated() {

                return calculated;

        }

        /**

         * N?mero de bandas

         * @return

         */

        public int getBandCount() {

                int bandCount = 0;

                for(int i = 0; i < statList.length; i ++)

                        bandCount += statList[i].getBandCount();

                return bandCount;

        }

}