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 org.gvsig.fmap.dal.coverage.datastruct.BufferHistogram;

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

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

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

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

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

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

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

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

import org.gvsig.raster.impl.datastruct.BufferHistogramImpl;

/**

 * Clase para la gesti?n de histogramas de un raster formado por una lista de ficheros.

 * Para devolver un histograma pedir? el histograma a todos los ficheros que

 * componen el multifichero.

 * 

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

 */

public class MosaicDataStoreHistogram implements Histogram {

        /**

         * Histograma de la imagen completa

         */

        private BufferHistogram                       histogram   = null;

        private Histogram[]                           multihistog = null;

        private MosaicRasterStore                     dataset     = null;

        /**

         * Constructor

         * @param dataset

         */

        public MosaicDataStoreHistogram(Histogram[] h, MosaicRasterStore dataset){

                this.multihistog = h;

                this.dataset = dataset;

        }

        /**

         * Gets the histogram. Builds the histogram of each dataset and join these

         * @return histograma 

         */

        public BufferHistogram getBufferHistogram()

                throws HistogramException, ProcessInterruptedException {

                if(histogram == null) {

                        BufferHistogram[] pres = new BufferHistogram[multihistog.length];

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

                                pres[i] = multihistog[i].getBufferHistogram();

                        }

                        Statistics stats = ((RasterDataStore)dataset).getStatistics();

                        try {

                                stats.calculate();

                        } catch (FileNotOpenException e) {

                                throw new HistogramException("");

                        } catch (RasterDriverException e) {

                                throw new HistogramException("");

                        }

                        BufferHistogram res = new BufferHistogramImpl(((RasterDataStore)dataset).getBandCount(), stats.getMin(), stats.getMax(), ((RasterDataStore)dataset).getDataType()[0]);

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

                                res.union(pres[i]);

                }

                return histogram;

        }

        /**

         * Pone a cero el porcentaje de progreso del proceso de calculo de histograma

         */

        public void resetPercent() {

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

                        multihistog[i].resetPercent();

        }

        /**

         * Obtiene el porcentaje de proceso en la construcci?n del histograma,

         * @return porcentaje de progreso

         */

        public int getPercent() {

                //TODO: Calculo del porcentaje

                return 100;

        }

        /*

         * (non-Javadoc)

         * @see org.gvsig.fmap.dal.coverage.dataset.Histogram#getMaximum()

         */

        public double getMaximum() {

                double max = 0;

                try {

                        for (int i = 0; i < getBufferHistogram().getNumBands(); i++) {

                                double m = getBufferHistogram().getMax(i);

                                if(m > max)

                                        max = m;

                        }

                } catch (HistogramException e) {

                        return 0;

                } catch (ProcessInterruptedException e) {

                        return 0;

                }

                return max;

        }

        /*

         * (non-Javadoc)

         * @see org.gvsig.fmap.dal.coverage.dataset.Histogram#getMinimum()

         */

        public double getMinimum() {

                double min = Double.POSITIVE_INFINITY;

                try {

                        for (int i = 0; i < getBufferHistogram().getNumBands(); i++) {

                                double m = getBufferHistogram().getMin(i);

                                if(m < min)

                                        min = m;

                        }

                } catch (HistogramException e) {

                        return 0;

                } catch (ProcessInterruptedException e) {

                        return 0;

                }

                return min;

        }

}