Application: gvSIG desktop

/* gvSIG. Desktop Geographic Information System.

 *

 * Copyright ? 2007-2017 gvSIG Association

 *

 * 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.

 *

 * For any additional information, do not hesitate to contact us

 * at info AT gvsig.com, or visit our website www.gvsig.com.

 */

package org.gvsig.raster.lib.buffer.impl.operations.tailtrim;

import org.gvsig.fmap.geom.exception.CreateEnvelopeException;

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

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

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

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

import org.gvsig.raster.lib.buffer.api.exceptions.BandException;

import org.gvsig.raster.lib.buffer.api.exceptions.BufferException;

import org.gvsig.raster.lib.buffer.api.exceptions.BufferOperationException;

import org.gvsig.raster.lib.buffer.api.operations.OperationFactory;

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

import org.gvsig.raster.lib.buffer.spi.exceptions.ProcessingOperationException;

import org.gvsig.raster.lib.buffer.spi.operations.AbstractSpecifiedBandsOperation;

import org.gvsig.tools.locator.LocatorException;

/**

 * @author fdiaz

 *

 */

public class TailTrimOperation extends AbstractSpecifiedBandsOperation{

    static public String STATISTICS_PARAM = "statistics";

    static public String TAIL_TRIM_PERCENT_PARAM = "tail_trim_percent";

    private Statistics statistics;

    private double tailTrimPercent;

    private RowProcessor[] rowProcessors;

    /**

     * @param factory

     *

     */

    public TailTrimOperation(OperationFactory factory) {

        super(factory);

    }

    @SuppressWarnings("unchecked")

    @Override

    public void preProcess() throws BufferOperationException {

        super.preProcess();

        BufferManager manager = BufferLocator.getBufferManager();

        statistics = (Statistics) this.getParameter(STATISTICS_PARAM, null);

        if (statistics == null) {

            statistics = this.getInputBuffer().getStatistics(null);

        };

        tailTrimPercent = (double) this.getParameter(TAIL_TRIM_PERCENT_PARAM, 0);

        tailTrimPercent = (tailTrimPercent>100)?100:tailTrimPercent;

        tailTrimPercent = (tailTrimPercent<0)?0:tailTrimPercent;

        //FIXME: Falta la gesti?n del par?metro copyUnprocessedBands, de momento se copian sin tenerlo en cuenta

        int bands = this.getInputBuffer().getBandCount();

        rowProcessors = new RowProcessor[bands];

        int [] bandTypes = new int[bands];

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

            bandTypes[i] = this.getInputBuffer().getBandTypes()[i];

        }

        NoData[] noData = this.getInputBuffer().getBandNoData();

        NoData[] resultNoData = new NoData[noData.length];

        for (int band = 0; band < noData.length; band++) {

            int bandType = this.getInputBuffer().getBand(band).getDataType();

            switch (bandType) {

            case BufferManager.TYPE_BYTE:

                rowProcessors[band] = new ByteRowProcessor(band);

                if (isProcessableBand(band)) {

                    if (noData[band].isDefined()) {

                        resultNoData[band] = noData[band];

                    } else {

                        resultNoData[band] = manager.createNoData((byte) 0, (byte) 0);

                    }

                }

                break;

            case BufferManager.TYPE_USHORT:

                rowProcessors[band] = new UShortRowProcessor(band);

                if (isProcessableBand(band)) {

                    if (noData[band].isDefined()) {

                        resultNoData[band] = noData[band];

                    } else {

                        resultNoData[band] = manager.createNoData((int) 0, (int) 0);

                    }

                }

                break;

            case BufferManager.TYPE_SHORT:

                rowProcessors[band] = new ShortRowProcessor(band);

                if (isProcessableBand(band)) {

                    if (noData[band].isDefined()) {

                        resultNoData[band] = noData[band];

                    } else {

                        resultNoData[band] = manager.createNoData(Short.MIN_VALUE, Short.MIN_VALUE);

                    }

                }

                break;

            case BufferManager.TYPE_INT:

                rowProcessors[band] = new IntRowProcessor(band);

                if (isProcessableBand(band)) {

                    if (noData[band].isDefined()) {

                        resultNoData[band] = noData[band];

                    } else {

                        resultNoData[band] = manager.createNoData(Integer.MIN_VALUE, Integer.MIN_VALUE);

                    }

                }

                break;

            case BufferManager.TYPE_FLOAT:

                rowProcessors[band] = new FloatRowProcessor(band);

                if (isProcessableBand(band)) {

                    if (noData[band].isDefined()) {

                        resultNoData[band] = noData[band];

                    } else {

                        resultNoData[band] = manager.createNoData(Float.MIN_VALUE, Float.MIN_VALUE);

                    }

                }

                break;

            case BufferManager.TYPE_DOUBLE:

                rowProcessors[band] = new DoubleRowProcessor(band);

                if (isProcessableBand(band)) {

                    if (noData[band].isDefined()) {

                        resultNoData[band] = noData[band];

                    } else {

                        resultNoData[band] = manager.createNoData(Double.MIN_VALUE, Double.MIN_VALUE);

                    }

                }

                break;

            default:

                throw new IllegalArgumentException("Unknow type of band '"+band+"'");

            }

        }

        try {

            this.setOutputBuffer(manager.createBuffer(

                this.getInputBuffer().getRows(),

                this.getInputBuffer().getColumns(),

                bandTypes,

                resultNoData,

                this.getInputBuffer().getProjection(),

                this.getInputBuffer().getEnvelope()));

        } catch (LocatorException | BufferException | CreateEnvelopeException e) {

            throw new BufferOperationException(e);

        }

    }

    @Override

    public void process() throws ProcessingOperationException {

        super.process();

        for (int band=0; band<this.getInputBuffer().getBandCount(); band++){

            if (getBandsToProcess().contains(band)) {

                Band bufferBand = this.getInputBuffer().getBand(band);

                Band outputBufferBand = this.getOutputBuffer().getBand(band);

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

                    Object rowBuffer = bufferBand.createRowBuffer();

                    bufferBand.fetchRow(row, rowBuffer);

                    Object outputRowBuffer = outputBufferBand.createRowBuffer();

                    outputBufferBand.fetchRow(row, outputRowBuffer);

                    rowProcessors[band].processRow(rowBuffer, outputRowBuffer);

                    outputBufferBand.putRow(row, outputRowBuffer);

                }

            } else {

                try {

                    this.getOutputBuffer().getBand(band).copyFrom(this.getInputBuffer().getBand(band));

                } catch (BandException e) {

                    throw new ProcessingOperationException(e);

                }

            }

        }

    }

    @Override

    public void postProcess()  throws BufferOperationException {

        super.postProcess();

    }

    interface RowProcessor {

        void processRow(Object inputRow, Object outputRow);

        Number processValue(Number value);

    };

    private abstract class AbstractRowProcessor implements RowProcessor {

        int band;

        double minValue;

        double maxValue;

        NoData noData;

//        NoData outputNoData;

        public AbstractRowProcessor(int band) {

            this.band = band;

            noData = getInputBuffer().getBand(band).getNoData();

//            outputNoData = outputBuffer.getBand(band).getNoData();

            double[][] tailTrim = statistics.getTailTrimValue(tailTrimPercent);

            minValue = tailTrim[this.band][0];

            maxValue = tailTrim[this.band][1];

        }

    }

    private class ByteRowProcessor extends AbstractRowProcessor {

        public ByteRowProcessor(int band) {

            super(band);

        }

        @Override

        public void processRow(Object inputRow, Object outputRow) {

            byte[] inputByteRow = (byte[])inputRow;

            byte[] outputByteRow = (byte[])outputRow;

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

                outputByteRow[i] = processValue(inputByteRow[i]).byteValue();

            }

        }

        @Override

        public Number processValue(Number value) {

            if(noData.isDefined() && noData.getValue().equals(value)){

                return getOutputBuffer().getBand(band).getNoData().getValue();

            }

            int iValue = 0xFF & ((Byte) value).byteValue();

            if(iValue < minValue || iValue > maxValue){

                return getOutputBuffer().getBand(band).getNoData().getValue();

            } else {

                return value;

            }

        }

    }

    private class ShortRowProcessor extends AbstractRowProcessor {

        public ShortRowProcessor(int band) {

            super(band);

        }

        @Override

        public void processRow(Object inputRow, Object outputRow) {

            short[] inputShortRow = (short[])inputRow;

            short[] outputShortRow = (short[])outputRow;

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

                outputShortRow[i] = processValue(inputShortRow[i]).shortValue();

            }

        }

        @Override

        public Number processValue(Number value) {

            if(noData.isDefined() && noData.getValue().equals(value)){

                return getOutputBuffer().getBand(band).getNoData().getValue();

            }

            int iValue = ((Short) value).shortValue();

            if(iValue < minValue || iValue > maxValue){

                return getOutputBuffer().getBand(band).getNoData().getValue();

            } else {

                return value;

            }

        }

    }

    private class UShortRowProcessor extends AbstractRowProcessor {

        public UShortRowProcessor(int band) {

            super(band);

        }

        @Override

        public void processRow(Object inputRow, Object outputRow) {

            short[] inputShortRow = (short[])inputRow;

            short[] outputShortRow = (short[])outputRow;

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

                outputShortRow[i] = processValue(inputShortRow[i]).shortValue();

            }

        }

        @Override

        public Number processValue(Number value) {

            if(noData.isDefined() && noData.getValue().equals(value)){

                return getOutputBuffer().getBand(band).getNoData().getValue();

            }

            int iValue = 0xFFFF & ((Short) value).shortValue();

            if(iValue < minValue || iValue > maxValue){

                return getOutputBuffer().getBand(band).getNoData().getValue();

            } else {

                return value;

            }

        }

    }

    private class IntRowProcessor extends AbstractRowProcessor {

        public IntRowProcessor(int band) {

            super(band);

        }

        @Override

        public void processRow(Object inputRow, Object outputRow) {

            int[] inputByteRow = (int[])inputRow;

            int[] outputByteRow = (int[])outputRow;

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

                outputByteRow[i] = processValue(inputByteRow[i]).intValue();

            }

        }

        @Override

        public Number processValue(Number value) {

            if(noData.isDefined() && noData.getValue().equals(value)){

                return getOutputBuffer().getBand(band).getNoData().getValue();

            }

            int iValue = value.intValue();

            if(iValue < minValue || iValue > maxValue){

                return getOutputBuffer().getBand(band).getNoData().getValue();

            } else {

                return value;

            }

        }

    }

    private class FloatRowProcessor extends AbstractRowProcessor {

        public FloatRowProcessor(int band) {

            super(band);

        }

        @Override

        public void processRow(Object inputRow, Object outputRow) {

            float[] inputFloatRow = (float[])inputRow;

            float[] outputFloatRow = (float[])outputRow;

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

                outputFloatRow[i] = processValue(inputFloatRow[i]).floatValue();

            }

        }

        @Override

        public Number processValue(Number value) {

            if(noData.isDefined() && noData.getValue().equals(value)){

                return getOutputBuffer().getBand(band).getNoData().getValue();

            }

            float fValue = value.floatValue();

            if(fValue < minValue || fValue > maxValue){

                return getOutputBuffer().getBand(band).getNoData().getValue();

            } else {

                return value;

            }

        }

    }

    private class DoubleRowProcessor extends AbstractRowProcessor {

        public DoubleRowProcessor(int band) {

            super(band);

        }

        @Override

        public void processRow(Object inputRow, Object outputRow) {

            double[] inputByteRow = (double[])inputRow;

            double[] outputByteRow = (double[])outputRow;

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

                outputByteRow[i] = processValue(inputByteRow[i]).doubleValue();

            }

        }

        @Override

        public Number processValue(Number value) {

            if(noData.isDefined() && noData.getValue().equals(value)){

                return getOutputBuffer().getBand(band).getNoData().getValue();

            }

            double fValue = value.doubleValue();

            if(fValue < minValue || fValue > maxValue){

                return getOutputBuffer().getBand(band).getNoData().getValue();

            } else {

                return value;

            }

        }

    }

}