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.cache.buffer.impl.histogram;

import org.gvsig.raster.cache.buffer.Buffer;

import org.gvsig.raster.cache.buffer.exception.HistogramException;

import org.gvsig.raster.cache.buffer.exception.ProcessInterruptedException;

import org.gvsig.raster.cache.buffer.histogram.BufferHistogram;

import org.gvsig.raster.cache.buffer.histogram.HistogramComputer;

/**

 * This class computes the histogram for a <code>Buffer</code>

 * 

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

 */

public class BufferHistogramComputer implements HistogramComputer {

        private BufferHistogram             histogram            = null;

        private Buffer                      buffer               = null;

        protected int                                    progressHistogram    = 0;

        private double[]                    limits               = null;

        private boolean                     refresh              = false;

        /**

         * Constructor

         * @param dataset

         */

        public BufferHistogramComputer(Buffer buffer) {

                this.buffer = buffer;

        }

        /*

         * (non-Javadoc)

         * @see org.gvsig.fmap.dal.coverage.store.props.HistogramComputer#getBufferHistogram()

         */

        public BufferHistogram getBufferHistogram(int numberOfClasses) throws ProcessInterruptedException, HistogramException {

                if(histogram == null || refresh) {

                        //RasterTask task = RasterTaskQueue.get(Thread.currentThread().getId() + "");

                        progressHistogram = 0;

                        double[][] limits = buffer.getStatistics().getAllBandsLimits();

                        BufferHistogramImpl hist = new BufferHistogramImpl(buffer.getBandCount(), limits[0], limits[1], buffer.getDataType(), numberOfClasses);

                        for (int iBand = 0; iBand < buffer.getBandCount(); iBand++) {

                                for (int row = 0; row < buffer.getHeight(); row++) {

                                        switch(buffer.getDataType()) {

                                        case Buffer.TYPE_BYTE:

                                                for (int col = 0; col < buffer.getWidth(); col++) 

                                                        hist.incrementPxValue(iBand, (double)(buffer.getElemByte(row, col, iBand)));

                                                break;

                                        case Buffer.TYPE_SHORT:

                                                for (int col = 0; col < buffer.getWidth(); col++) 

                                                        hist.incrementPxValue(iBand, (double)(buffer.getElemShort(row, col, iBand)));

                                                break;

                                        case Buffer.TYPE_INT:

                                                for (int col = 0; col < buffer.getWidth(); col++) 

                                                        hist.incrementPxValue(iBand, (double)(buffer.getElemInt(row, col, iBand)));

                                                break;

                                        case Buffer.TYPE_FLOAT:

                                                for (int col = 0; col < buffer.getWidth(); col++) 

                                                        hist.incrementPxValue(iBand, (double)buffer.getElemFloat(row, col, iBand));

                                                break;

                                        case Buffer.TYPE_DOUBLE:

                                                for (int col = 0; col < buffer.getWidth(); col++) 

                                                        hist.incrementPxValue(iBand, buffer.getElemDouble(row, col, iBand));

                                                break;

                                        }

                                        /*if (task.getEvent() != null)

                                                task.manageEvent(task.getEvent());*/

                                        progressHistogram = ((iBand * buffer.getHeight() + row) * 100) /(buffer.getHeight() * buffer.getBandCount());

                                }

                        }

                        progressHistogram = 100;

                        histogram = hist;

                }

                return histogram;

        }

        /*

         * (non-Javadoc)

         * @see org.gvsig.fmap.dal.coverage.store.props.HistogramComputer#resetPercent()

         */

        public void resetPercent() {

                progressHistogram = 0;

        }

        /*

         * (non-Javadoc)

         * @see org.gvsig.fmap.dal.coverage.store.props.HistogramComputer#getPercent()

         */

        public int getPercent() {

                return progressHistogram;

        }

        /*

         * (non-Javadoc)

         * @see org.gvsig.fmap.dal.coverage.store.props.HistogramComputer#getMinimum()

         */

        public double getMinimum() {

                if(limits == null)

                        try {

                                limits = buffer.getStatistics().getLimits();

                        } catch (ProcessInterruptedException e) {

                                return 0;

                        }

                return limits[0];

        }

        /*

         * (non-Javadoc)

         * @see org.gvsig.fmap.dal.coverage.store.props.HistogramComputer#getMaximum()

         */

        public double getMaximum() {

                if(limits == null)

                        try {

                                limits = buffer.getStatistics().getLimits();

                        } catch (ProcessInterruptedException e) {

                                return 0;

                        }

                return limits[1];

        }

        /*

         * (non-Javadoc)

         * @see org.gvsig.fmap.dal.coverage.store.props.HistogramComputer#refreshHistogram()

         */

        public void refreshHistogram() {

                refresh = true;

        }

}