Application: gvSIG desktop

package org.gvsig.raster.lib.buffer.impl.statistics;

import org.gvsig.raster.lib.buffer.api.Band;

import org.gvsig.raster.lib.buffer.api.BufferManager;

import org.gvsig.raster.lib.buffer.api.statistics.HistogramBand;

import org.gvsig.raster.lib.buffer.api.statistics.StatisticsBand;

import org.gvsig.tools.ToolsLocator;

import org.gvsig.tools.dispose.DisposeUtils;

import org.gvsig.tools.dispose.impl.AbstractDisposable;

import org.gvsig.tools.exception.BaseException;

import org.gvsig.tools.task.SimpleTaskStatus;

/**

 * @author fdiaz

 *

 */

public class DefaultStatisticsBand implements StatisticsBand {

    private boolean calculated;

    private double mean;

    private double max;

    private double secondMax;

    private double min;

    private double secondMin;

    private double variance;

    private HistogramBand histogramBand;

    private long countNoDataValues;

    private int bandLenght;

    private int dataType;

    /**

     * @param band

     */

    public DefaultStatisticsBand(Band band) {

        init();

    }

    /**

     *

     */

    private void init() {

        this.calculated = false;

        this.countNoDataValues = 0;

    }

    @Override

    public boolean isCalculated() {

        return calculated;

    }

    @Override

    public void calculate(SimpleTaskStatus status, Band band) {

        boolean isMyStatus = false;

        if (status == null) {

            status = ToolsLocator.getTaskStatusManager().createDefaultSimpleTaskStatus("Statistics");

            status.add();

            status.message("Calculating statistics band");

            isMyStatus = true;

        } else {

            status.push();

        }

        try {

            calculated = false;

            status.setRangeOfValues(0, band.getRows() * 2);

            long cont = 0;

            this.bandLenght = band.getRows() * band.getColumns();

            this.dataType = band.getDataType();

            MeanOperation meanOperation = new MeanOperation();

            MaxOperation maxOperation = new MaxOperation();

            MinOperation minOperation = new MinOperation();

            SecondMaxOperation secondMaxOperation = new SecondMaxOperation();

            SecondMinOperation secondMinOperation = new SecondMinOperation();

            meanOperation.pre();

            maxOperation.pre();

            minOperation.pre();

            secondMaxOperation.pre();

            secondMinOperation.pre();

            for (int row = 0; row < band.getRows(); row++) {

                status.setCurValue(cont++);

                if (status.isCancelled()) {

                    status.abort();

                    return;

                }

                for (int column = 0; column < band.getColumns(); column++) {

                    Double value = band.getAsDouble(row, column);

                    if (value == null) {

                        this.countNoDataValues++;

                    } else {

                        meanOperation.addValue(value);

                        maxOperation.addValue(value);

                        minOperation.addValue(value);

                        secondMaxOperation.addValue(value);

                        secondMinOperation.addValue(value);

                    }

                }

            }

            meanOperation.post();

            maxOperation.post();

            minOperation.post();

            secondMaxOperation.post();

            secondMinOperation.post();

            this.mean = ((Number) meanOperation.getResult()).doubleValue();

            this.max = ((Number) maxOperation.getResult()).doubleValue();

            this.min = ((Number) minOperation.getResult()).doubleValue();

            this.secondMax = ((Number) secondMaxOperation.getResult()).doubleValue();

            this.secondMin = ((Number) secondMinOperation.getResult()).doubleValue();

            VarianceOperation varianceOperation = new VarianceOperation(this);

            varianceOperation.pre();

            HistogramOperation histogramOperation = new HistogramOperation(this, band);

            histogramOperation.pre();

            for (int row = 0; row < band.getRows(); row++) {

                status.setCurValue(cont++);

                if (status.isCancelled()) {

                    status.abort();

                    return;

                }

                for (int column = 0; column < band.getColumns(); column++) {

                    Double value = band.getAsDouble(row, column);

                    if( value != null ) {

                        histogramOperation.addValue(value);

                        varianceOperation.addValue(value);

                    }

                }

            }

            histogramOperation.post();

            varianceOperation.post();

            this.histogramBand = ((HistogramBand) histogramOperation.getResult());

            this.variance = ((Number) varianceOperation.getResult()).doubleValue();

            if (isMyStatus) {

                status.terminate();

            } else {

                status.pop();

            }

            calculated = true;

        } catch (Exception e) {

            status.abort();

            throw e;

        }

    }

    @Override

    public int getDataType() {

        return this.dataType;

    }

    @Override

    public long getBandLength() {

        return this.bandLenght;

    }

    @Override

    public double getSecondMax() {

        return this.secondMax;

    }

    @Override

    public double getSecondMin() {

        return this.secondMin;

    }

    @Override

    public double getMaximum() {

        return this.max;

    }

    @Override

    public double getMinimum() {

        return this.min;

    }

    @Override

    public double getMean() {

        return this.mean;

    }

    @Override

    public double getVariance() {

        return this.variance;

    }

    @Override

    public double[] getTailTrimValue(double percent) {

        return getHistogramBand().getTailTrimValue(percent);

    }

    @Override

    public double[] getTailTrimValue(int pos) {

        return getHistogramBand().getTailTrimValue(pos);

    }

    @Override

    public HistogramBand getHistogramBand() {

        return this.histogramBand;

    }

    @Override

    public double getMedian() {

        return getHistogramBand().getMedian();

    }

    @Override

    public long getNoDataValuesCount() {

        return countNoDataValues;

    }

}